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CHAPTER  I 


INTRODUCTION 

The  purpose  of  this  investigation  was  to  make  a  detailed  study  of 
the  fundamental  properties  of  PS  converted  waves  to  ascertain  whether 
they  can  be  detected  consistently  on  seismograms  from  both  underground 
nuclear  explosions  and  earthquakes  and  whetner  they  can  be  used  to  help 
distinguish  underground  nuclear  explosions  from  earthquakes.  The  study 
included : 

a.  Theoretical  considerations  of  the  amplitudes  of  PS  and  other  con 
verted  waves  as  compared  to  observed  amplitudes. 

b.  Harmonic  analysis  of  seismograms  to  aid  in  the  interpretation. 

c.  Analysis  cf  s,  el  sinograms  of  underground  nuclear  explosions  and 
large,  distant  earthquakes  recorded  at  the  University  of  Utah  seismograph 
stations  and  other  selected  stations  to  develop  dependable  criteria  for 
recognition  of  PS  converted  waves  on  seismograms, 

d.  The  analysis  of  selected  seismograms  from  underground  nuclear  ex 
plosions  for  recognition  of  FS  converted  waves  recorded  by  vertical-  and 
horizontal-component  seismometers  at  the  Geotech  field  stations. 
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General  background 


A  PS  converted  wave  ts  a  seismic  body  wave  which  results  from  the 
conversion  of  an  incident  parent  P  wave  at  a  boundary  within  the  earth's 
crust  or  upper  mantle  to  a  refracted  vertically  polarized  shear  wave,  SV. 
Other  types  of  conversion  are  also  possible,  as  for  example,  the  conver¬ 
sion  of  an  incident  SV  wave  to  a  P  wave. 

Although  this  study  was  originally  intended  to  concentrate  on  the  P 
to  S  conversion,  it  was  recognized  during  the  course  of  the  investigation 
that  other  types  of  conversion  x^ithin  the  earth's  crust  and  upper  mantle 
are  also  of  importance  in  seismology  and  therefore  also  in  the  probieS* of 
distinguishing  earthquakes  and  underground  nuclear  explosions.  Consequently 
this  study  was  expanded  to  include  some  of  these  other  types  of  conversions. 
It  should  be  emphasized,  however,  that  our  study  is  concerned  principally 
with  those  types  oi  conversions  that  would  arrive  before  the  conventional 
body  shear  wave  and  for  the  most  part  within  approximately  one  minute  after 
the  arrival  of  the  initial  P  wave. 

Most  of  the  previous  work  by  Soviet  and  American  investigators  on  P 

to  S  conversions  has  been  summarized  and  evaluated  by  Cook,  Algermissen, 

and  Costain  (1962,  p,  4775).  As  of  1962,  the  tentative  results  of  the 
research  on  PS  converted  waves,  principally  for  epicentral  distan’os  less 
than  400  km,  can  be  summarized  as  follows, 

1.  The  amplitude  of  the  PS  converted  waves  were  usually  several  times 

greater  than  that  of  the  first  arrival  of  the  P  wave,  but  not  without  ex¬ 

ception. 

2.  The  successively  later  arrivals  of  PS  converted  waves  from 
successively  deeper  horizons  showed  progressively  greater  ground-velocity 
amplitudes . 
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3.  For  blasts,  the  periods  of  the  various  PS  converted  waves  were 
approximately  the  same  as  those  of  the  parent  P  wave  (Schwind  and  others, 
1960);  whereas  for  earthquakes,  the  periods  of  the  PS  converted  waves 
were  1.3  to  1.7  times  larger  than  those  of  the  parent  P  wave  (Kuz'mina, 
1959). 

It  should  be  emphasized  that  it  is  desirable  to  compile  particle- 
motion  diagrams  or  similar  aids  to  the  identification  of  the  PS  conver¬ 
ted  waves.  As  of  1962,  particle  motion  diagrams  in  PS  converted  waves 
studies  had  been  published  only  by  the  American  workers  (Schwind  and 
others,  1960).  Of  the  Soviet  workers,  Kuz'mina  (1959)  was  the  only  one 
who  reported  compiling  particle-motion  diagrams,  but  his  were  not  pub¬ 
lished. 

The  tentative  interpretations  that  have  been  made  in  the  limited 
time  available  on  this  research  project  indicate  that  converted  body 
waves  may  constitute  a  recognizable  part  of  the  early  arriving  seismic 
energy  within  about  30  sec  after  the  first  arrival  of  the  energy  of  the 
compressional  P  wave  and  before  the  arrival  of  the  conventional  body  shear 
wave  at  epicentral  distances  up  to  about  10,000  km.  We  include  in  conver¬ 
ted  waves  the  energy  that  is  both  refracted  and  reflected  at  boundaries 
within  the  earth's  crust  and  upper  mantle. 
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Data  necessary  in  the  study  of  PS  converted  waves 

To  facilitate  the  identification  of  PS  converted  waves  in  particu¬ 
lar  and  converted  waves  in  general,  it  is  desirable  to  know  the  crustal 
section  and  velocities  in  the  region  of  study.  Even  with  this  informa¬ 
tion  available,  considerable  difficulty  may  be  experienced  with  multiple 
reflections.  This  problem  has  been  discussed  briefly  by  Cook,  Alger- 
missen,  and  Costain  (1962,  p.  4776). 

Particle-motion  diagrams  greatly  facilitate  the  identification  of 
PS  converter  waves  Often,  the  beginning  of  a  phase  on  a  seismogram  is 
difficult  to  determine  without  the  aid  of  a  particle-motion  diagram.  In 
a  comparison  of  observed  data  with  theory,  such  factors  as  near-surface 


layering  and  instrument  response  must  be  considered. 


CHAPTER  II 


EQUIPMENT  AND  INSTRUMENTATION 
Permanent  seismograph  stations 

Most  of  the  seismograms  from  permanent  seismograph  stations  used 
for  analysis  in  this  study  were  recorded  at  the  three  seismograph 
stations  of  the  University  of  Utah.  Table  2-2  gives  the  principal  facts 
of  these  stations  and  pertinent  information  concerning  the  instrumenta¬ 
tion,  including  the  magnifications.  The  Portable  Benioff  seismographs 
at  all  three  stations,  the  Wood-anderson  seismographs  at  Dugway  and 
Price,  and  allied  instrumentation  were  funded  largely  ±.rom  this  project; 
and  this  instrumentation  was  installed  prior  to  October  1,  1963. 

Some  seismograms  used  for  analysis,  especially  during  the  earlier 
part  of  the  study,  were  recorded  at  other  permanent  seismograph  stations. 
Table  2-2  gives  the  available  principal  facts  and  pertinent  information 
concerning  the  instrumentation  of  these  stations. 

For  the  analyses  of  seismograms  from  the  permanent  seismograph 
stations,  the  seismograms  were  usually  used  from  vert'c  il  and  horizontal 
seismometers  which  were  matched  and  which  also  had  essentially  the  same 
magnification.  In  those  cases  in  which  the  magnifications  of  the  matched 
instruments  differed,  the  digitized  values  as  measured  on  the  seismo¬ 
grams  were  normalized. 
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Table  2-1 --  Seismograph  stations  of  the  University  of  Utah 
Dugway,  Utah*  (DUG) 


Latitude:  40°  11.7’  N  Longitude:  112°  49.0’  W 

Elevation:  1^81  meters  Bedrock:  Rhyolite-quartz  dacite  flow  (Ter¬ 

tiary) 

Locality:  The  vault  is  located  5  miles  southwest  of  Dugway,  Utah  and 
is  underground. 

Time:  Standard  timing  system  of  the  World-Wide  Seismograph  System, 
Model  10700,  The  Geotechnical  Corporation,  Garland,  Texas. 

Time  signal  impressed  every  6  hours. 

Instruments:  Standard  Benioff  three-component  seismographs  (Z,N.E)**; 
To  *  1  sec,  Tg  ■  0.75  sec,  Mag.  *  400,000#,  trace 
speed  60  mm  per  minute;  Portable  Benioff  three-compon¬ 
ent  seismographs  (Z,N,E)***, 

To  =  1  sec,  Tg  *  1.5  sec,  Mag.  *  100,000#,  trace  speed 
60  or  120  mm/min;  Wcod-Anderson  horizontal  torsion 
seismographs  (N,E)***,  To  -  0.8  sec,  Mag.  =  2,800, 
trace  speed  60  mm  per  minute. 

Price,  Utah  (PCU) 

Latitude:  39°  36.4'  N  Longitude:  110°  48.3’  W 
Elevation:  1714  meters  Bedrock:  Mancos  shale  (Cretaceous) 

Locality:  The  vault  is  located  on  the  carpus  of  Carbon  College, 

Price,  Utah 

Time:  Simplex  Master  Clock  calibrated  against  WWV.  Time  signals 
impressed  manually  twice  daily. 

Instruments:  Portable  Benioff  three-component  seismographs  (Z,N,E)****; 

To  »  1  sec,  Tg  *  1.5  sec,  Mag.  **  36,100  (Z)  and  14,000 
(N,E),  trace  speed  60  or  120  mm/min;  Wood-Anderson 
horizontal  torsion  seismographs  (N,E);  To  *  0.8  sec, 

Mag.  ■  2,800,  trace  speed  60  mm/min. 


*  One  of  the  stations  in  the  World-Wide  Network  of  Standardized 
Seismographs  operated  in  conjunction  with  the  U.  S.  Coast  & 
Geodetic  Survey. 

**  Instruments  placed  in  operation  May  6,  1962. 

***  Instruments  placed  in  operation  on  September  30,  1963. 

****  Instruments  placed  in  operation  on  May  16,  1962. 

#  All  three  components  (Z,N,E). 
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Table  2-1 (continued)--  Seismograph  stations  of  the  University  of  Utah 

Salt  Lake  City,  Utah  (SLC) 

Latitude:  40°  45.9'  N  Longitude*  1110  50.9’  W 
Elevation:  1425  meters  Bedrock:  Quaternary  alluvium 
Locality:  The  vault  is  located  on  the  campus  of  the  University  of 
Utah,  Salt  Lake  City,  Utah 

Time:  IBM  Master  Clock,  calibrated  against  WWV.  WWV  receiver 

equipped  with  radio  time-signal  converter  and  amplifier  for 
impression  of  time  signals  directly  on  each  seismogram. 

Time  signal  impressed  twice  daily. 

Instruments:  Portable  Ben.ioff  three-component  seismographs 

(Z,N,E) *****;  To  «  1  sec,  Tg  *  1.5  sec,  Mag.  -  30,000 
(Z)  and  25,000  (N,E),  trace  speed  60  or  120  mm/min. 


*****  Instruments  placed  in  operation  on  April  8,  1962. 
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Table  2-2--  Other  permanent  seismograph  stations,  from  which  seismo¬ 
grams  were  used  in  this  study. 

Boulder  City,  Nevada  (BCN) 

Latitude:  35°  58'  51"  N 
Longitude:  114°  50'  02"  W 
Elevation:  776  meters 


Instrumentation  January  15,  1962  (written  communication  Jan.  15, 
from  Capt.  R.  A.  Earle,  U.  S.  Coast  and  Geodetic  Survey) 


3  Benioff  Moving  Coil  Short  Period 

N-S 

0.45  sec 
1.4  sec 
Critical 
50,000  ca 


Flaming  Gorge,  Utah  (FGU) 

Latitude:  40°  55*  35.5"  N 
Longitude:  109°  23*  10.3"  W 
Elevation:  1982.5  meters 


Tg 

To 

Seismometer  damping 
Magnification 
(Operating) 


0.45  sec 
1.4  sec 
Near  critical 
100,000  ca 


E-W 
0.45  sec 
1.4  sec 
Critical 
50,000  ca 


Note:  The  above  applies  to  film  magnified  8  times. 


Instrumentation  January  8,  1962  (written  communication  January  8, 

1962  from  Capt.  R.  A.  Earle,  U.  S.  Coast  and  Geodetic  Survey) 


3  Benioff  Moving  Coil  Short  Period 

N-S 
0.35  sec 
1.41  sec 
3:1 

Near  critical 
65,000  ca 

60 

Note:  Avove  applies  to  film  magnified  8  times. 


Tg 
To 

Galv.  Damping 
Seismo.  Damping 
Magnification 
(Operating) 

Relative  dial  setting  Mag. 


0.33  sec 
1.06  sec 
6:1 

Near  critical 
150,000  ca 


60 


E-W 
0.58  sec 
1,48  sec 
Near  critical 
3:1 

65,000  ca 
60 
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Table  2- 2 (continued) . --  Other  permanent  seismograph  stations,  from 

which  seismograms  were  used  in  this  study. 


La  Paz,  Bolivia  iLPB) 

Latitude:  16°  31'  57. 6"  S 
ongitude:  68®  05*  54.1  W 
Elevation:  3292  meters 

Instrumentation:  Standardized  seismographs,  World-wide  network  of 

U.  S.  Coast  and  Geodetic  Survey 

Magnification  of  Z,  E-W,  and  N-S:  approximately  50,000  at  1  cps. 

La  Palma,  El  Salvador  (LPS) 

Latitude:  14°  17'  32.0"  N 
Longitude:  89°  09'  43.0  W 
Elevation:  1000  meters 

Instrumentation:  Standardized  seismographs,  World-wide  network  of 

U.  S.  Coast  and  Geodetic  Survey 

Magnification  of  E-W,  and  N-S:  approximately  100,000  at  1  cps. 
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Stations  of  the  Long-Range  Seismic  Measurements  Program 

(Geotech  stations) 


Sites  in  the  western  part  of  the  United  States  recording  under¬ 
ground  nuclear  explosions  at  the  Nevada  test  site  are  shown  in  Fig¬ 
ure  5-1  ,  These  sites  are  occupied  as  part  of  the  Long-Range  Seismic 
Measurements  (LRSM)  program  and  are  operated  by  the  Geotechnical  Corpor¬ 
ation  as  part  of  the  VELA  UNIFORM  program. 

Data  are  recorded  in  the  field  on  an  Ampex  14-channel  magnetic- 
tape  system.  Model  314,  and  are  frequency-modulated  with  a  270-cps 
center-frequency  carrier.  Field  recordings  are  made  at  a  tape  speed  of 
0.3  ips. 

Selected  segments  of  data  from  each  of  these  LRSM  stations  for  each 
event  are  re-Tecorded  onto  a  single  reel  of  tape  to  form  a  magnetic  tape 
composite  of  a  single  event.  A  typical  composite  contains  data  from  about 
forty  stations,  recorded  in  the  order  of  increasing  distance  from  the  sig¬ 
nal  epicenter.  Approximately  30  minutes  of  data  are  recorded  from  each 
station.  Both  short-period  and  long-period  data  are  recorded. 

Three  short-period  seismometers  at  each  LRSM  site,  either  large  or 
small  Benioff  seismometers,  are  oriented  in  vertical,  radial  and  transverse 
directions  with  respect  to  the  azimuth  to  the  epicenter.  The  response  of 
the  short-period  system  peaks  at  about  2.6  cps.  Magnifications  varied 
from  site  to  site. 

We  have  further  edited  the  magnetic  tape  composites  of  certain  under¬ 
ground  nuclear  explosions,  and  put  the  well-recorded  short-period  data  on 
new  composites,  which  were  later  digitized.  This  is  discussed  more  fully 
in  Chapter  IV. 
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CHAPTER  III 
THEORETICAL  STUDIES 


Energy  and  amplitude  ratios  of  converted  waves 

A  study  was  made  of  amplitude  ratios,  energy  ratios,  and  phase  angles 
of  reflected  and  refracted  1*  and  SV  waves  derived  from  SV  waves  or  P  waves 
incident  on  a  plane  interface.  All  computations  were  programmed  in  FORTRAN 
for  the  IBM  7040  and  1620  systems.  These  programs  are  listed  in  the 
Appendix;  Program  18  deals  with  incident  P  waves  and  Program  19  with  inci¬ 
dent  SV  waves.  Incident  angles  were  varied  from  0°  to  88®.  Both  real 
and  imaginary  coefficients  were  considered  in  the  calculations.  Reflec¬ 
tion  and  transmission  coefficients  were  calculated  for  P  and  SV  waves 
for  various  velocity  and  density  ratios. 

Results  of  a  plane  SV  wave  incident  on  a  plane  interface  are  des- 
scribed  in  Costain,  Cook  and  Algermissen  (1963;  1965)  and  by  Costain  and 
Cook  (1965).  In  these  publications,  energy  ratios  from  Knott's  energy 
equation  and  amplitude  ratios  and  phase  angles  from  the  Zoeppritz  equa¬ 
tions  were  calculated  for  a  plane  SV  wave  incident  on  a  plane  elastic  dis¬ 
continuity  using  Program  19.  Incident  angles  were  varied  from  0®  to  88° 
in  increments  of  2®  except  near  the  critical  angles,  and  near  phase  changes 
of  180°  in  the  real  domain,  where  the  ratios  were  calculated  in  incre¬ 
ments  of  0.25®.  Compress ional  velocity  ratios  (transmitted/incident)  of 
0.7,  0.8  and  0.9,  and  density  ratios  (transmitted/incident)  of  0.7,  0.8, 

0.9  and  1.0  were  used.  Poisson's  ratio  was  0.25  for  each  medium.  Some 
calculations  for  a  Poissons  ratio  of  0.400  are  given  in  Costain  and  Cook 
(1965). 

A  layered  crustal  section  was  investigated  and  a  FORTRAN  program 
written  (Program  9A)  to  calculate  the  theoretical  amplitudes  to  be  expec- 
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ted  from  converted  and  unconverted  P  and  S  waves  on  the  basis  of  reflec¬ 
tion  coefficients  from  the  Zoeppritz  equations.  The  program  assumes  layer 
thicknesses  large  compared  with  the  wavelength  involved.  Some  results  are 
presented  in  Costain,  Cook  and  Algermissen  (1963). 

Theoretical  travel-time  curves 

Theoretical  travel-time  curves  were  computed  (using  Program  DR-13 
listed  in  the  Appendix)  for  many  crustal  models  in  the  western  part  of 
the  United  States.  These  curves  include  all  P  and  S  head  waves  from  each 
interface,  and  converted  head  waves  such  as  PPPSS,  SSPSS,  SSPPP,  etc.  The 
models  are  based  on  seismic  refraction  studies  by  the  United  States  Geolo¬ 
gical  Survey,  University  of  Utah,  and  others.  The  travel-time  curves  were 
used  to  aid  in  the  interpretation  of  the  events  on  the  seismograms  out 
to  the  arrival  of  the  body  S-wave.  In  addition,  certain  reflected  and 
multiply  reflected  phases  were  added  to  the  travel-time  curves  (Appendix: 
Programs  5  and  21) . 

Other  theoretical  work  was  primarily  of  a  nature  to  aid  in  the  inter¬ 
pretation  of  the  early  part  of  the  seismogram,  and  included: 

a)  Computation  of  the  power  density  spectrum  from  the  cosine 
transform  uf  the  autocorrelation  function  (Program 

No.  DR  12-B) . 

b)  Computation  of  head  wave  coefficients  for  a  PPP  head  wave 
using  the  method  of  Zvolinskii  (Program  No.  1). 

A  complete  list  of  the  theoretical,  utility,  and  data  analysis  pro¬ 
grams  is  contained  in  the  Appendix. 
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Theoret  *r.  1  ParticX  e*  D  lag  r  aiyfio  tion  l 

Particle-motion  in  an  infinite,  homogeneous,  isotropic  medium  was 
examined  theoretically  for  certain  simple  combinations  of  P  and  S  motion. 

The  following  diagrams  were  constructed  according  to  the  diagram  of 
Figure  3-la. 

Fig.  3-lb:  Two  P  waves  of  amplitudes  A^  and  A2 

(Ai  -  10;  A2  =  12)  arriving  at  a  point  at  different  angles  of 
incidence,  i  (Ai  *  10°).  Both  wave  trains  have  the  same 
frequency  (1  cps), but  differ  in  phase  by  an  amount  ^  l(o.  The 
x-component  of  motion  is  given  by: 

x  =  10  sin  (2  TT  t)  sin  10°  +  12  sin  (2  TT  t  +  7r/fe)  sin  70°. 
The  vertical,  or  y-component  of  motion  is  given  by: 

y  =  10  sin  (2  77  t)  cos  10°  +  12  sin  (2  77  t  +  '%)  cos  70° 
The  resultant  motion  is  shown  in  Figure  3- lb  . 

Fig.  3-lo:  Pure  P  and  pure  S  arriving  simultaneously  with  the  S-amplitude 

(As  =  15)  greater  than  the  P-amplitude  (Ap  =  10).  The  two  wave 
trains  arrive  at  the  same  angle  of  incidence  (i  -  20°)  but 
differ  in  phase  by  an  amount  Both  have  a  frequency  of 

1  cps.  The  equations  that  describe  the  x-  and  y-  components 
of  motion  are 

x  =  10  sin  (2 TT' t)  cos  20°  -  15  sin  (2  7T t  +  %)  sin  20° 
y  =  10  sin  (2  7Tt)  sin  20°  +  15  sin  (2*^t  +  cos  20® 
The  resultant  motion  is  shown  in  Figure  3-lc. 

Fig*  3-ld:  pure  P  and  pure  S  arriving  simultaneously  with  continuously 
decreasing  amplitudes,  the  amplitude  of  S  decreasing  faster 
than  the  amplitude  of  P.  Initially  As>Ap.  Both  wave  trains 
arrive  at  the  same  angle  of  incidence  (i  s  30°)  and  have  the 

-tf 

same  frequency  (1  cps),  but  differ  in  phase  by  an  amount  . 
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Figure  3-1.  Theoretical  vertical  particle-motion  diagrams, 


0.4t 

x  =  (20/  <2  )  sin  (2  t)  cos  30® 

-  (50/e°'6t)  sin  (2-^t  +  %)  sin  30® 
y  =  (20/e°‘4t)  sin  ( 2 ^  t)  sin  30® 

+  (50/e°’6t)  sin  2^t  +  "%)  cos  30® 

The  resultant  motion  is  shown  in  Figure  3-ld. 

Fig,  3-lir:  P  and  S  both  decreasing  in  amplitude  at  different  rates 

(As/'Ap),  arriving  at  the  same  angle  of  incidence,  differing 
in  phase  by  an  amount”/^  ,  and  with  the  frequency  of  the  P-wave 
equal  to  twice  the  frequency  of  the  S-wave. 
x  =  (20/e°‘4t)  sin  (4“^t)  cos  30®  - 

(50 /e°*6t)  sin  (2^1  +~%)  sin  30® 
y  =  (20/ e.0'^)  sin  (4"7/*t)  sin  30®  + 

(50/<^°-6t)  sin  ( 2  ~TT t  +  *|)  cos  30® 

6 

The  resultant  motion  is  shown  in  Figure  3-1*". 
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Fig.  3-lf:  and  S  arriving  simultaneously  followed  by  P2  after  a 


time  t  (t  s  0.45  sec) . 

As  decreases  faster  than  Ap2  which  decreases  faster  than  Ap^. 

?1  and  S  arrive  at  the  same  angle  of  incidence;  ?2  arrives  at 
a  slightly  different  angle.  All  arrivals  have  the  same  fre¬ 
quency  but  differ  in  phase.  This  example  is  shown  in  Figure  3-lf. 
Fig.  3-lg:  P  and  S  arriving  at  same  angle  of  incidence.  Both  waves  have 
the  same  constant  amplitude  (As  =  Ap)  but  differ  in  phase  and 
frequency  (frequency  of  S  s  6/5  cps;  frequency  of  P  =  4/5  cps). 
x  =  10  sin  (8  7Tt/5)  cos  30°-10  sin  (12  ^t/5  +  ^)  sin  30° 
y  =  10  sin  (8  7?t/5)  sin  30e-10  sin  (12  "77 1/5  +  Tf)  cos  30° 
This  example  is  shown  in  Figure  3-lg. 

Fig.  3-lh:  Pi*  ?2  anc*  S  all  with  the  same  constant  amplitude. 

Pi  and  S  arrive  at  the  same  angle  of  incidence. 

?1  and  P2  are  in  phase  with  each  other,  but  out  of  phase  with  S. 
All  have  different  frequencies  (Pi  =  4/5  cps;  P2=8/5  cps; 

S  =  6/5  cps). 

x  -  10  sin  (S'*- 1/5)  cos  30°-10  sin  (12^  t/5  +  J)  sin  30° 

4-  10  sin  (16  ^ t /5)  cos  (  77"/4.5) 
y  =  10  sin  (8 ^ t/5)  sin  30°  +  10  sin  (12^t/5)  +  ”? )  cos  30° 

+  10  sin  (16  ^t/5)  sin  ("^/4.5) 
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The  plots  in  Figure  3-1  are  representative  of  at  least  a  few  of  the 
types  of  particle-motion  diagrams  that  can  be  expected  in  the  simplest 
cases.  The  examples  described  above  do  not  include  the  effects  of  in¬ 
strument  response,  the  effects  of  a  free  surface,  the  effects  of  layer¬ 
ing,  inhomogeneities,  or  anisotropy.  Diagrams  similar  to  those  in  Figure 
3-1  are  frequently  observed  on  actual  particle-motion  diagrams  of  ground 
motion.  However,  this  should  not  imply  that  the  theoretical  and  observed 
particle-motion  diagrams  that  resemble  each  other  have  a  common  origin. 
The  simple  theoretical  s.udy  she  -s  that  differences  in  wave  type,  fre¬ 
quency,  phase,  and  rates  of  change  of  amplitude  of  phases  arriving 
simultaneously  at  a  point  does  give  a  combined  affect  which  resembles 
actual  particle -motion  diagrams  in  many  instances.  Further  study  is 
necessary  to  determine  the  relative  importance  of  each  of  these  effects, 
as  well  as  the  effects  of  a  layered  half-space. 

Figure  3-2  shows  four  representative  particle-motion  diagrams  that 
in  sonic  respects  closely  resemble  the  theoretical  diagrams  of  Figure  3-1. 
All  four  diagrams  show  motion  in  a  vertical-radial  plane.  Diagrams  (a), 
(b),  (c),  and  (d)  of  Figure  3-2  are  to  be  compared  with  diagrams  (e) , 

(g),  (d)  and  (h) ,  respectively,  of  Figure  3-1. 
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Figure  3-2.  Observed  vertical  particle-motion  diagrams. 

(a)  Particle-motion  diagram  of  the  AARDVARK  event  as  recorded  at 
station  MNNV.  Time:  Pn  +  16.95  sec. 

(b)  Particle-motion  diagram  of  the  STILLWATER  event  as  recorded  at 
station  WINV.  Time:  Pn  +  15.96  sec. 

(c)  Particle-motion  diagram  of  the  STILLWATER  event  as  recorded  at 
station  WINV.  Time:  Pn  +  24.24  sec. 

(d)  Particle-motion  diagram  of  the  STILLWATER  event  as  recorded  at 
station  WINV.  Time:  Pn  +  18.04  sec. 


CHAPTER  TV 


METHOD  OF  ANALYSIS  OF  DATA  RECORDED  BY  THE  GEOTECH  FIELD  STATIONS 

Data  available  for  analysis 

Data  recorded  by  the  Geotech  field  teams  were  received  from  several 
sources.  Data  in  digital  form  on  1/2"  IBM-compatible  computer  tapes 
were  received  from  the  Seismic  Data  Laboratory  in  Alexandria,  Virginia, 
These  data  included 


Explos ion 

Station 

AARDVARK 

MVCL 

AARDVARK 

FMUT 

AARDVARK 

LCNM 

AARDVARK 

DRCO 

AARDVARK 

MNNV 

AARDVARK 

PMWY 

FISHER 

WINV 

FISHER 

MVCL 

FISHER 

FMUT 

FISHER 

MNNV 

FISHER 

LCNM 

GNOME 

LCNM 

GNOME 

TPNV 

GNOME 

SVAZ 

GNOME 

FSAZ 

GNOME 

WMAZ 

In  addition  to  the  above,  a  total  of  seventeen  analog  1"  FM  magnetic 
tape  composites,  eacn  containing  20  to  40  analog  FM  recordings  of  a  single 
explosion  recorded  at  different  stations, were  edited  and  spliced,  and  a 
new  composite  made.  Editing  was  done  by  the  Geophysics  Department,  Univer¬ 
sity  of  Utah;  splicing  and  playback  to  get  the  final  new  composite  was 
done  by  Space  Technology  Laboratories,  Redo'ido  Beach,  California,  The 
original  seventeen  tape  composites  were  copies  of  recordings  made  by  the 
Geotech  field  teams  across  the  United  States  for  the  following  seventeen 
explosions : 

STOAT  GNOME 

FISHER  AARDVARK 
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MAD 

CHINCHILLA 

FEATHER 

PLATYPUS 

C0DSAW 

HARD HAT 

AGOUTI 

CIMMARON 

CHENA 

ARMADILLO 

MINK 

STILLWATER 

PAMPAS 

Time  segments  containing  the  best  short-period  data  from  the  above 
seventeen  explosions  were  recombined  on  a  new  composite.  The  new  com¬ 
posite  reduced  the  number  of  1"  reels  from  seventeen  to  three.  These 
three  reel,  contain  all  of  the  best  short-period  data  for  each  of  the 
above  seventeen  explosions.  Short-period  calibration  data  were  put  on 
two  additional  composite  reels.  The  signal  composite  also  contains 
long-period  data  since  the  original  1"  analog  FK  magnetic  tapes  were 
fourteen- track  recordings  of  both  short-period  and  long-period  data; 
however,  when  editing,  emphasis  was  placed  on  the  short-period  data. 

Digitizing  intervals  used 

The  three  1"  analog  FM  composite  reels  were  then  digitized  and  the 
digitized  data  placed  on  1/2"  computer  tape  for  subsequent  processing 
on  the  IBM  7040  digital  computer.  The  digitizing  was  done  by 

a)  Hercules  Powder  Co. ,  Salt  Lake  City,  Utah 
(digitizing  interval  0.02  sec) 

b)  Space  Technology  Laboratories,  Redondo  Beach,  Calif, 
(digitizing  interval  0.03  sec) 

A  description  of  the  contents  of  the  ti  r-ee  composite  analog  FM  mag¬ 
netic  tapes  made  by  the  University  of  Utah  and  Space  Technology  Labora¬ 
tories  is  on  the  following  pages. 
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Composite  Analog  Tapes 
of  Data  Recorded  by 
The  Geotech  Field  Teams 


Composite  Analog  Reel  No.  1 


GNOME 


LCNM, 

POTX, 

SSTX, 

TCNM, 

SMTX, 

FSAZ, 

WMAZ, 

KNOT, 

WNSD, 

SFAZ, 

LMNV, 

cwar, 

WRAR, 

JSTN, 

MMTN, 

MVCL, 

NGWS, 

BLWV, 

DHNY, 

BGME. 

AARDVARK 

LMNV, 

MNNV, 

ATNV, 

KNOT, 

TFCL, 

VTOR, 

HLID, 

DRCO, 

LCNM, 

PTOR, 

MFAR, 

CNWS, 

MMTN, 

BLWV, 

GLTX, 

FISHER 

LMNV, 

MNNV, 

KNOT, 

WMAZ, 

FMUT, 

TCNM, 

VNUT, 

MLNM, 

LCNM, 

PMWY, 

MAD 

TPNV, 

KNOT, 

WMAZ, 

WINV, 

LMNV, 

HLID, 

PMWY, 

SSTX, 

WNSD, 

SMTX, 

FEATHER 

MNNV, 

WNSD, 

KNOT, 

HBOK. 

WMAZ, 

FMOT, 

WINV, 

MINK 

WINV, 

VNUT, 

MNNV, 

BFCL, 

FMOT, 

CHENA 

MNNV, 

CPCL, 

FSAZ, 

MVCL, 

VNUT, 

MLNM, 

LPTX, 

SVAZ, 

HBOK, 

AMOK, 

TOOK, 

DRCO, 

SJTX, 

PMWY, 

MPAR, 

VNUT, 

FMUT, 

CPCL, 

TPNV, 

MNNV, 

HLID, 

CVTN, 

WINV, 

FMUT, 

WINV, 

CPCL, 

FSAZ, 

MVCL, 

VNUT, 

PMWY, 

HSNB, 

WNSD, 

LPTX, 

TOOK, 

SEMN, 

SJTX, 

CWAR, 

JSTN, 

NGWS, 

DANY, 

BGME. 

WINV, 

CPCL, 

FSAZ, 

MVCL, 

SFAZ, 

HLID, 

SSTX, 

AMOK, 

LPTX. 

MNNV, 

FMOT, 

CPCL, 

FSAZ, 

MVCL, 

VNUT, 

MPAR, 

BLWV. 

CPCL, 

MVCL, 

VNUT, 

DRCO, 

PMWY, 

SSTX, 

FSAZ. 


HLID. 
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Composite  Analog  Reel  No.  2 


/ 


HARD HAT 

DVCL,  MNNV,  ATNV,  KNOT,  TNCL,  WMAZ,  FMUT,  WINV,  CPCL,  FSAZs  MVCL, 
SFAZ,  SFAZ-CC,  VNUT,  VTOR,  SVAZ,  HLID,  DRCO,  MLNM,  TCNM,  BFCL,  LCNM, 
EPTX,  TCNM,  PTOR,  PMWY,  RTNM,  EFTX,  GNNM,  BMTX,  SSTX,  LPTX. 

PAMPAS 

DVCL,  MNNV,  ATNV,  ATNV-CC,  BFCL,  KNUT,  TNCL,  TNCL-CC,  WMAZ,  FMUT, 
WINV,  CPCL,  FSAZ,  SFAZ,  VTOR,  HLID,  TCNM,  MVCL,  VNUT,  SVAZ,  DRCO, 
LCNM,  KNUT-CC,  WMAZ-CC,  SVAZ-CC,  PTOR,  EPTX,  EFTX,  SSTX. 

CIMARRON 

DVCL,  ATNV,  BFCL,  KNUT,  TNCL,  WMAZ,  FMUT,  WINV,  FSaZ,  SFAZ,  VTOR, 

HLID,  MNNV,  CPCL,  MVCL,  VNUT,  SVAZ ,  DRCO,  TCNM,  PTOR,  LCNM,  EPTX, 

EFTX. 

AGOUTI 

DVCL,  MNNV,  MNNV-CC,  ATNV,  ATNV-CC,  KNUT,  TNCL,  TNCL-CC,  WMAZ,  FMUT, 

WINV,  FSAZ,  SFAZ,  DVCL-CC,  BFCL,  KNUT-CC,  WINV-CC,  CPCL,  MVCL,  VNUT, 

SVAZ,  HLID,  TCNM,  BFCL-CC,  FMUT-CC,  CPCL-CC,  FSAZ-OC,  MVCL-CC, 
HLID-CC,  DRCO,  DRCO-CC,  MLNM,  KLNM-CC,  TCNM-CC,  LCNM,  LCNM-CC,  GNNM, 


J 

CC  -  denotes  cavity  collapse 
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ARMADILLO 


DVCL ,  DVCL-CC ,  MNNV,  ATNV,  ATNV-CC,  KNOT,  KNUT-CC,  TNCL ,  TNCL-CC, 
WMAZ ,  WMAZ-CC,  FMUT,  FSAZ,  VTOR,  WINV,  MNNV-CC,  BFCL ,  FMUT-CC, 
CPCL,  MVCL ,  SFAZ,  VTOR-CC,  HI ID,  TCNM,  WINV-CC,  CPCL-CC,  FSAZ-CC, 
MVCL-CC,  SFAZ-CC,  HLID-CC,  MLNM,  MLNM-CC ,  TCNM-CC,  PTOR,  PTOR-CC, 
LCNM,  LCNM-CC- ,  PMWY,  PtfWY-CC,  EPTX,  EPTX-CC,  GNNM,  GNNM-CC,  BMTX, 

VNUT  VNUT-CC,  SVAZ,  SVAZ-CC,  DRCO,  DRCO-CC,  BFCL-CC, 

CODSAW 

DVCL,  MNNV,  KNUT,  TNCL,  WMAZ,  WINV,  VTOR,  BFCL,  FSAZ,  SFAZ,  HLID, 
TCNM,  DVCL-CC,  MNNV-CC,  KNUT-CC,  TNCL-CC,  WMAZ-CC,  FMUT,  WINV-CC, 
CA-CPCL,  CPCL,  MVCL,  VNUT,  VTOR-CC,  SVAZ,  DRCO,  LPTX. 

CHINCHILLA 

DVCL,  ATNV,  KNUT,  TNCL,  WMAZ  WINV,  VTOR,  MNNV,  CW-BFCL,  BFCL, 
FMUT,  CW-FSAZ.  FSAZ,  CW-SFAZ,  SFAZ,  CW-HLID,  CS-TCNM,  TCNM,  CPCL, 

CW-MVGL,  MVCL.  CW-VNUT,  VNUT,  CW-SVAZ,  SVAZ,  HLID,  CW-DRCO,  DRCO, 

LCNM,  EPTX,  SSTX. 

STILLWATER 

DVCL,  MNNV,  ATNV,  ATNV-CC,  KNUT,  TNCL,  WMAZ,  FMUT,  WINV,  VTOR, 
MNNV-CC,  BFCL,  KNUT-CC,  TNCL-CC,  WINV-CC,  FSAZ,  MVCL,  SFAZ,  HLID 
BFCL-CC,  FMUT-CC,  CPCL,  CPCL-CC,  FSAZ-CC,  MVCL-CC,  SFAZ-CC, 

VNUT,  VNUT-CC ,  VTOR-CC,  SVAZ,  HLID-CC,  DRCO,  MLNM,  MLNM-CC,  TCNM, 
TCNM-CC ,  LCNM,  EFTX,  SSTX,  DHNT. 

STOAT 

DVCL,  KNUT,  WMAZ,  FMUT,  WINV,  MNNV,  CPCL,  FSAZ,  MVCL,  SFAZ,  HLID, 
TCNM,  VNUT,  SVAZ,  KLNM,  DRCO,  LCNM, 

PLATYPUS 

DVCL,  MNNV,  ATNV,  WMAZ ,  FSAZ,  SSTX. 


CC  -  denotes  cavity  collapse 
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The  new  analog  composite  was  not  completed  in  time  to  process  more 
than  a  small  fraction  of  the  data  on  the  three  reels.  Most  of  the  data 


that  were  processed  were  studied  using  the  alphanumeric  and  CalComp 
particle-motion  plots  described  below. 

Other  data  recorded  by  the  Ceotech  field  teams  which  were  received 
but  not  processed  by  us  in  any  way  because  of  insufficient  time  include: 

a)  Eighteen  underground  explosions  on  reels  of  1"  FM  magnetic 
tape  received  in  March  1964  from  the  Seismic  Data  labora¬ 
tory. 

b)  Sixteen  earthquakes  and  underground  explosions  recorded 
by  the  Geotech  recording  teams  and  digitized  by  the 
Seismic  Data  Laboratory.  These  data  are  in  digital 
form  on  1/2"  IBM-campatible  computer  tape. 

Only  underground  explosions  were  studied  from  the  recordings  made 
by  the  Geotech  field  teams. 
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Particle-motion  diagrams  were  made  directly  from  digitized  data.  No 


corrections  for  instrument  lespons.  other  than  magnification  were  made. 

The  particle-motion  diagrams  were  displayed  by  one  of  the  following  methods; 
a)  Hand  plots.  The  diagrams  were  plotted  ny  hand  during  the  early 

stages  of  the  project  since  machine  methods  were  not  readily  avail¬ 
able  to  us.  The  following  particle-motion  diagrams  were  construct¬ 


ed  using  this  method; 


Explos ion 
AARDVARK 


FISHER 

GNOME 


Station 


MNNV, 

ATNV, 

KNUT, 

TFCL, 

FMUT, 

WINV, 

CPCL, 

F5AZ, 

MVCL, 

VNUT, 

FTOR, 

HLID, 

DRCO, 

LCNM, 

PTOR, 

PMWY, 

MNNV, 

LC.  .i, 

MVCL, 

VTOR, 

NCID 

KNUT, 

Mv'CL, 

LCNM 

MMTN, 

LPTX, 

WNSD, 

SJTX, 

POTX, 

SSTX, 

TCNM, 

SMTX, 

MLNM, 

LPTX, 

SVAZ, 

AMOK, 

FSAZ, 

WMAZ, 

KNUT, 

PMWY, 

MPAR, 

VNUT, 

FMUT, 

CWAR, 

WRAR, 

b)  Alphanumeric  plots  -  These  are  machine  plots  done  on  the  IBM  7040 
system  with  an  IBM  1403  printer  as  output.  This  method  greatly 
increased  plotting  efficiency  over  the  hand  plots,  but  the  com¬ 
puter  facilities  required  for  obtaining  this  type  of  plot  were 
not  available  to  us  until  March  1964.  The  program  itself  was 
not  operational  until  September  1964.  The  method  prints  out  a 
known,  preselected  sequence  of  letters  and  numbers  on  the  IBM 
1403  printer  at  positions  on  the  paper  corresponding  to  the 
digitized  data.  The  known  sequence  of  symbols  is  then  connect¬ 
ed  by  hand  to  get  the  final  particle-motion  diagram.  The  data 
for  the  particle-motion  diagrams  are  read  directly  from  digi¬ 
tized  data  on  magnetic  tape  or  from  cards,  using  Program  DR-18FP1 
or  Program  DR-18FP2  respectively,  described  in  the  Appendix.  An 
example  of  an  alphanumeric  plot  is  shown  in  Figure  4-1 
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Figure  4-1.  Particle-motion  diagram  illustrating  the  alphanumeric  type  of  plot. 


Table  4-1  Summary  of  particle-motion  plots  of  data  recorded  at  the  ueotech  stations 
Epicentral  Plot 

Distance  Yield  Interval  Magnification  (K) 

Explosion  Station  (KM)  Plot*  Gain^  (KT . )  Seic~*  (Sec. )  Vertical _ Radial _ Tran 
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L  denotes  large  lenioff  seismometer;  S  denotes  small  Benioff  seismometer . 


The  following  particle-motion  diagrams  were  constructed  using  this  method 


Exp los ion  Station 


CIMMARON,  HARDHAT,  PAMPAS,  STOAT  FMUT 
AARDVARK,  CHINCHILLA,  CODSAW,  STOAT  MNNV 
AGOUTI,  CHINCHILLA,  STILLWATER,  STOAT  WINV 
AGOUTI  (Cavity  collapse)  ATNV 
AGOUTI  CPCL 
AGOUTI  SFAZ 


c)  CalComp  plots .  The  University  of  Utah  Computer  Center  acquired  an 
IBM  1627  (CalComp)  plotter  in  October  1964,  and  plotting  techniques 
using  data  on  magnetic  tape  were  perfected  in  December  1964.  This 
method  is,  of  course,  preferable  to  either  of  the  above  for  dis¬ 
playing  particle-motion  diagrams. 

The  following  particle-motion  diagrams  were  constructed  using  this  method 
Explosion  Station 

ARMADILLO  FSAZ,  SFAZ,  WMAZ 

In  general,  particle-motion  diagrams  were  constructed  only  for  the 
first  20  to  30  seconds  of  motion  after  the  arrival  of  the  initial  P  wave. 
In  a  few  cases  particle-motion  diagrams  were  made  for  the  time  interval 
around  the  expected  arrival  of  the  body  S-wave  to  help  identify  this  S 
phase. 

Particle-motior  diagrams  were  always  made  for  both  a  horizontal 
plane  and  a  vertical-radial  plane. 

A  typical  sequence  of  particle-motion  diagrams  plotted  on  a  CalComp 
plotter  is  shown  in  Chapter  5,  Figure  5-10 

A  summary  of  all  particle-motion  plots  for  the  Geotech  stations  for 
which  we  have  at  least  two  events  per  station  is  given  in  Table  4-1  . 
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CHAPTER  V 


INTERPRETATION  OF  DATA  RECORDED  AT  THE  CEOTECH  FIELD  STATIONS 

Introduction 

A  limited  study  of  the  theoretical  amplitudes  of  P  to  S  and  S  to  P 
conversions  was  made  during  the  early  part  of  this  project.  This  is  dis¬ 
cussed  in  Chapter  III.  The  theory  de^lt  with  plane  waves,  but  more  sig¬ 
nificantly,  with  wave  lengths  that  were  small  compared  with  layer  thick¬ 
nesses.  More  recent  results  of  the  Haskel 1-Thompson  matrix  method  have 
shown  the  importance  of  thin,  low-velocity  surface  layers  on  the  motion 
at  the  surface  arising  from  a  plane  compressional  wave  incident  at  the 
base  of  a  layered  half-space.  In  particular,  Hannon  (1964)  has  con¬ 
structed  a  number  of  theoretical  seismograms  using  frequency -dependent 
transmission  coefficients  which  may  be  computed  by  the  Haskell-Thompson 
matrix  method.  The  method  determines  the  ratio  of  the  vertical  and 
horizontal  motions  at  the  surface  to  the  total  amplitude  (or  the  ampli¬ 
tudes  of  the  components)  of  the  plane  wave  incident  at  the  bottom  inter¬ 
face  of  the  layered  system  (Haskell,  1953;  Hannon,  1964). 

The  amplitude  ratios  are  functions  of  frequency,  angle  of  incidence, 
layer  thickness,  compressional  and  shear  velocities  in  each  layer,  and 
densities  of  the  layers.  The  steady-state  transmission  coefficients  give 
the  amplitude  and  phase  of  the  vertical  and  horizontal  surface  motion  to 
the  motion  of  the  steady-state  plane  wave  incident  at  the  Dottora  inter¬ 
face.  They  include  the  important  effects  of  frequency  and  layer  thick¬ 
ness  on  the  transmission  coefficients,  af  well  as  the  angle  of  incidence, 
velocities,  and  densities  of  the  layers  that  affect  the  amplitude  ratios 
as  computed  from  the  Zoeppritz  equations  (Chapter  III,  this  report). 
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Hannon  (1964)  used  four  crustal  models  with  vary iu /  thicknesses  of 
alluvial  layers  ranging  from  zero  to  30,  to  730  to  1430  meters  in  thick¬ 
ness.  It  is  worthwhile  here  to  summarize  his  conclusions  (Hannon,  1964, 
p.  2072): 

1)  The  presence  of  less  than  100  feet  of  sediments  causes  the  surface 
amplitudes  to  increase  by  a  factor  of  two  for  the  source  functions 
used  by  him.  Ringing  was  a  predominant  feature. 

2)  For  thicker  alluvial  sections,  the  surface  motion  amplitudes  're 
even  greater,  and  reached  values  three  times  those  for  the  model 
with  no  alluvium. 

3)  "The  near  surface  layeiing  exerts  a  greater  influence  on  the 
horizontal  component  than  on  the  vertical,  causing  the  horizon¬ 
tal  component  to  have  a  more  irregular  shape.  This  emphasizes 
the  importance  of  converted  SV  waves.  This  effect  is  further 
illustrated  by  the  fact  that  the  arrival  times  predicted  by  ray 
theory... do  not  agree... as  well  for  the  horizontal  as  they  do  for 
the  vertical.  Thus  the  presence  of  layering  may  produce  ellip¬ 
tical  particle  motion." 

4)  "The  different  character  introduced  by  the  multiple  and  conver¬ 
ted  waves  tends  co  introduce  a  longer  period  oscillation  in  the 
record  which  is  dependent  on  the  structure  rather  than  on  the 
source." 

It  is  apparent  from  Hannon's  work  that  care  must  be  used  in  inter¬ 
preting  shear  motion  on  a  particle-motion  diagram,  and  in  particular,  in 
identifying  the  shear  on  the  basis  of  its  amplitude  relative  to  other  ob¬ 
served  phases  on  the  seismogram.  The  effect  of  layered  structure  may 
drastically  affect  both  the  amplitude  of  the  shear,  and  the  apparent 
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arrival-time  on  a  component  seismogram.  As  Hannon  showed,  for  the  pulse 
widths  he  used,  the  near-surface  layering  exerted  a  greater  influence  on 
the  horizontal  component  than  on  the  vertical,  thus  emphasizing  the  im¬ 
portance  of  converted  SV  waves . 

In  the  interpretations  discussed  in  this  Chapter,  the  amplitudes 
of  certain  of  the  converted  waves,  for  example,  PPS,  seem  larger  than 
would  be  expected  for  the  assumed  velocity  contrasts  in  the  simplified 
crustal  sections  used.  However,  the  crustal  models  assumed  are  a  first 
approximation  only,  based  on  refraction  data  recorded, at  relatively  widely 
separated  stations  in  many  cases,  and  more  resolution  may  be  required  to 
detect  and  explain  the  effects  of  thin,  near  surface  layering  on  the 
amplitudes  of  the  arrivals,  although  the  models  used  do  agree  with  the 
arrival  times  of  the  phases  interpreted  on  the  seismograms,  in  most 
cases.  Also,  it  is  apparent  that  "ringing"  is  present  on  many  of  the 
seismograms,  which  may  oveiride  a  weaker  PS  converted  wave  following  a 
stronger  P  phase. 

In  view  of  the  theoretical  work  discussed  above  which  takes  into 
account  the  frequency  content  of  the  pulse  as  well  as  the  layer  velocities, 
densities  and  incident  angles,  the  interpretations  presented  in  this 
Chapter  must  be  regarded  as  provisional,  especially  where  the  interpreta¬ 
tion  of  shear  arriving  before  the  body  S  wave  is  concerned.  True  shear 
arrivals  during  this  time  would  have  to  be  a  result  of  P  to  S  conversion 
somewhere  along  the  transmission  path.  The  effect  of  the  recording  sta¬ 
tion  location  and  the  layering  beneath  the  station  on  the  amplitudes  of 
the  arrivals  is  now  recognized  as  having  an  important  effect  on  the  char¬ 
acter  of  the  seismogram.  Quantitative  measurements  of  amplitude  ratios  of 
S  to  P  or  P  to  S  or  even  P  to  P  from  particle-motion  diagrams  must  be  made 
with  care.  Further  study  of  the  enhancement  of  the  amplitudes  of  early 


shear  arrivals  as  a  function  of  frequency  content  of  the  pulse  is  warranted. 
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General  Discussion 


Figure  5-1  shows  the  geographic  locations  of  the  Geotech  recording 
teams  in  the  western  part  of  the  United  States  with  respect  to  the  Nevada 
test  site  (NTS).  We  were  not  able  to  interpret  all  of  the  data  we  have 
on  magnetic  tape  in  digital  form  recorded  at  these  stations.  A  few  events 
recorded  at  selected  stations  are  discussed  in  this  chapter. 

Figure  5-2  shews  the  locations  of  all  of  the  underground  explosions 
in  the  Nevada  test  site  for  which  particle-motion  diagrams  were  made  at 
at  least  one  Geotech  field  recording  station  for  this  study.  Table  5-1 
lists  some  of  the  pertinent  data  for  each  of  these  eleven  explosions. 
Emphasis  was  placed  on  Geotech  stations  at  which  two  or  more  underground 
explosions  were  recorded  at  the  same  site  from  the  same  approximate  source, 
that  is,  i'ne  Nevada  test  site.  Particle-motion  diagrams  were  also  made 
for  the  GNOME  event  in  New  Mexico,  but  no  other  data  are  available  for 
the  same  source  area  and  epicentral  distances,  and  so  these  particle- 
motion  diagrams  are  not  presented  herein. 

Figure  5-3  shows  CalComp  plots  of  three  underground  nuclear  explo¬ 
sions  recorded  at  Winnemucca,  Nevada  (WINV)  by  the  Geotech  field  record¬ 
ing  teams.  The  signals  fom  the  underground  nuclear  explosions  STOAT  and 
AGOUTI  are  very  similar.  The  digital- to-analog  playbacks  of  each  of  these 
two  recordings  are  shown  in  Figure  5-4  on  an  expanded  scale.  The  traces 
are  almost  identical  phase  for  phase  for  the  entire  fifteen  seconds  of 
the  signal.  Beyond  that,  the  signal  as  shown  is  poor  because  of  clipping 
in  the  digitizing  electronics  were  the  FM  analog  tapes  were  digitized. 

The  third  explosion,  STILLWATER,  in  Figure  5-3,  does  not  have  the  same 
degree  of  similarity  in  the  first  fifteen  seconds,  but  the  epicenter  for 
this  explosion  was  approximately  nine  kilometers  from  STOAT  and  AGOUTI. 
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Figure  5-1.  Location  map  of  LRSM  recording  stations  in  the  western  part 
of  the  United  States. 
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Figure  5-2.  Location  map  of  underground  nuclear  explosions  at  the 
Nevada  test  site. 
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Figure  5-3.  CalComp  digital -to-analog  playback  of  the  events  AGOUTI,  STILLWATER  and  STOAT  as  recorded  at 
station  WTNV, 
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First  15  sec  only 


See  Fig.  5pTT  Other  considerations  such  as  yield,  iepth  of  uriai  , 
geological  environment  at  the  source,  etc.  would  also  be  expected  to 
affect  the  source  function. 

Figure  5-5  shows  an  annotated  version  of  Figure  5-3  for  AGOUTI, 

STIL1  /ATER,  and  STOAT  as  recorded  at  WINV.  Several  things  are  apparent 
on  the  seismograms: 

1)  Patterns  resembling  beat  frequency  patterns  or  possibly 
"ringing"  are  well  developed  during  the  first  15  seconds 
of  the  seismogram  for  STILLWATER.  The  patterns  are  still 
evident  after  15  seconds  in  the  stronger  arrivals.  En¬ 
velopes  have  been  drawn  in  to  accentuate  the  patterns. 

Such  patterns  are  also  observed  for  AGOUTI  but  are  less 
well  developed.  Also,  clipping  in  the  digitizing  elec¬ 
tronics  has  obscured  the  patterns  after  the  first  15 
seconds  of  the  AGOUTI  seismogram. 

2)  For  the  seismogram  of  STILLWATER,  the  wave  train  labeled 
ABCDEFG  at  approximately  17.5  seconds  after  the  arrival  of 
Pn  is  very  similar  in  appearance  to  the  wave  train  labeled 
A'B'C'D' E'F'G'  at  approximately  24.5  seconds  after  the  arrival 
of  Pn>  Examination  of  the  particle-motion  diagrams  for 
these  two  wave  trains  show  the  ground  motion  for  cycle  A 

to  be  S-like  and  prograde,  whereas  the  ground  motion  for 
A1  is  P-like  and  prograde.  Cycles  C,  D,  C,  F,  and  G  show 
prograde  P-raotion,  but  C*  ,  D' ,  E',  and  F'  show  prograde 
S-motion,  for  the  most  part.  The  first  wave  train,  then, 
is  predominantly  P-motion,  the  later  wave  train  is  predomi¬ 
nantly  S-motion.  Both  wave  tra'.ns  are  of  approximately  the 
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Figure  5-5.  CalComp  digital-to-analog  payback  of  the  events  AGOUTI,  STILLWATER  and  STOAT  as  recorded  at 
station  WINV.  Annotated. 


same  amplitude,  as  shown  on  the  seismogram  recording 
the  vertical  component  of  ground  motion. 

3)  Wave  trains  A. ..G  and  A'...G'  on  the  STILLWATER  seismogram 
are  parts  of  short-period  events  which  modulate  another 
wave  whose  period  is  approximately  four  seconds.  The  ampli¬ 
tude  of  the  "carrier"  wave  seems  to  decrease  slightly  with 
increasing  time  and  the  wave  modu'ated  by  the  A1 ...G'  wave 
train  does  not  appear  to  be  a  s  pie  continuation  of  the 
one  modulated  by  wave  train  A. ..G.  The  time  of  onset  of 
the  wave  modulated  by  wave  train  A. ,.G  appears  to  be  at 
approximately  15.5  seconds  after  the  arrival  of  Pn  at 
WJ.NV. 

4)  A  spectral  analysis  plot  of  STILLWATER-WINV  shows  two 
dominant  peaks  in  the  frequency  spectrum,  one  at  1.8  cps 
and  one  at  2.3  cps.  These  two  peaks  seem  to  be  associated 
primarily  with  approximately  the  first  14  seconds,  and 

with  the  last  14  seconds,  respectively,  of  the  seismogram. 

Many  of  the  main  features  discussed  above  appear  in  other  seismograms 
obtained  by  the  Geotech  stations  that  we  have  studied  to  date.  The  modu¬ 
lated  low  frequency  "carrier"  waves  can  be  seen  for  both  STOAT  and  AGOUTI 
(Figure  5-5)  but  clipping  in  the  digitizing  electronics  for  AGOUTI  makes 
the  wave  difficult  to  see.  These  low-frequency  waves  are  recorded  over 
the  interval  of  time  during  which  multiple  reflections  are  arriving  from 
the  base  of  the  crust.  The  wave  as  sketched  in  for  theSTOAT  seismogram 
appears  to  be  of  higher  frequency  than  that  f  >r  either  STILLWATER  or  AGOUTI. 
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Figure  5-6  shows  an  annotated  series  of  three  underground  nuclear 
explosions  recorded  at  Fillmore,  Utah  (FMUT) .  The  "beat"  frequency 
patterns  and  low  frequency  waves  are  not  as  well  developed  as  they  are 
for  STILLWATER -WINV.  The  sites  for  the  three  explosions  were  widely 
separated  at  the  Nevada  test  site  (see  Figure  5-2),  and  the  seismograms 
bear  little  resemblance  to  one  another. 

Wave  trains  after  Pj^P  (or  P?)  on  the  vertical  seismograms  for 
CIMARRON  and  HARDHAT  on  Figure  5-6  for  station  FMUT  are  similar  to  wave 
trains  recorded  for  STILLWATER  at  WINV  (Figure  5-5).  However,  the 
CIMARRON -HARDHAT  group  is  not  as  well  developed  later  on  the  same  seismo¬ 
gram.  The  crustal  tranmission  paths  are  probably  quite  different  in  de¬ 
tail  for  STILLWATER -WINV  and  CIMARRON -FMUT.  There  is  a  decrease  in 
amplitude  for  first-motion  and  later  arrivals  arriving  at  FMUT  from  the 
Nevada  test  site  as  compared  with  corresponding  epicentral  distances  in 
Arizona  and  Nevada. 

The  beat-frequency  patterns  may  be  a  resonance  effect  due  to  the 
presence  of  a  low-velocity  sedimentary  surface  layer  for  which  the  re- 
flecti— i  coefficient  is  effectively  unity  at  the  top  of  the  layer.  Van 
Nostrand  (1964)  has  shown  that  "resonance  is  the  most  important  factor 
in  signal  amplitude  variation  resulting  from  changes  in  near  surface 
geology".  The  maximum  signal  amplitude  occurs  when  the  thickness  of 
the  low-velocity  surface  layer  approaches  a  quarter-wavelength.  For 
frequencies  of  1  cps  in  a  low-velocity  sedimentary  layer  having  a  P- 
wave  velocity  of  3.0  km  per  second,  maximum  signal  amplitudes  will 
be  obtained  for  layer  thickne  ses  of  775  meters,  or  about  2,500  feet. 

These  quarter-wavele  gt'n  resonant  thicknesses  will  considerably  distort 
the  signal  shape,  which  then  becomes  a  function  of  the  thickness  of  the 
thickness  of  the  surface  layer.  Van  Nostrand  (1964)  shows  the  original 
input  signal  considerably  increased  in  amplitude  for  a  surface  layer  whose 

thickness  is  A/4,  and  considerably  lengthened 
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(16-30:59.4) 


in  cine  for  a  surface  layer  whose  thickness  is  3 

The  resonance  effect  will  seriously  interfere  with  the  identifica¬ 
tion  of  weak  arrivals  closely  following  stronger  ones,  since  an  original 
two-cycle  event  may  be  drawn  out  into  as  many  as  eight  cycles  by  the 
resonance  effect. 
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Discussion  of  STILLWATER -WINV 


From  the  time  of  arrival  of  Pn  on  the  seismogram  of  STILLWATER-WlNV  out 
to  about  four  seconds,  the  motion  is  predominantly  P-like.  The  amplitude 
increased  gradually  almost  in  a  linear  fashion  out  to  approximately  3.5 
seconds.  Some  of  this  may  be  a  resonance  effect  due  to  "ringing"  asso¬ 
ciated  with  a  low-velocity  sedimentary  layer  beneath  the  station. 

Table  5-2  lists  theoretical  arrival  times  for  some  "phases"  which 
were  searched  for  on  the  seismogram  of  STILLWATER-WINV. 

The  interval  from  4.5  sec  to  6.5  sec  after  the  arrival  of  Pn  is 
characterized  by  about  four  cycles  of  predominantly  S-like  motion  with 
the  attitude  of  the  major  axis  of  the  motion  in  a  vertical  plane  gradually 
increasing  until  the  motion  is  nearly  vertical.  The  amplitude  increases 
in  a  linear  fashion  during  the  interval  of  A.  t  =  4.5  sec  to^t  =  approxi¬ 
mately  7.5  sec. 

The  motion  for  the  first  five  cycles  in  this  interval  is  alternately 
prograde  and  retrograde  for  each  cycle.  This  behavior  may  indicate  the 
simultaneous  arrival  of  wave  trains  of  different  type  and  frequency.  The 
alternating  prograde  and  retrograde  behaviour  in  the  theoretical  particle- 
motion  diagrams  was  a  characteristic  of  the  simultaneous  arrival  of  wave 
trains  of  slightly  different  frequency  content. 

At  about  6.5  sec  after  the  first  motion,  the  motion  becomes  large  in 
amplitude  and  nearly  vertical  and  remains  prograde  for  the  next  11  cycles 
until  approximately  12  sec  after  the  arrival  of  Pn,  The  amplitudes  reach 
a  maximum  at  about  7.5  sec  after  the  arrival  of  P„. 
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Table  5-2.-- Theory -al  arrival  times  of  certain  phases  for  the 
crusta.  i  action  of  Figure  5-7a  (STILLWATER-WINV) 


Phase 

pn 


Description 

First  arrival, 
head  wave 


Theoretical 
Arrival  Time 
(Sec) 


67.7 


Phase-Pn 
Arrival  Time, 
t ,  sec 


SPP 

and 

PPS 


Herd  waves, 
converted  at 
base  of  crust 


71.9 


A. 2 


SPS 


Head  wave,  up  and 
down  paths  shear 

Direct  wave 


76.1 

80.2 


8. A 
12.6 


PM1P_ 
(or  P?) 


P-wave  reflected 
once  from  base  of 
crust  (or  a  guided 
wave  in  the  crust) 


80.9 


13.2 


PMlppPp 


P-wave  reflected  once 
from  base  of  crust  at 
critical  angle,  then 
reflected  from  the  surface 
of  earth  back  to  base 
of  crust  where  it  con¬ 
tinues  as  a  P  head  wave 


73.6 


5.9 


PM2P 

?M3P 


P-wave  reflected  twice 

from  base  of  crust  82.6 

P-wave  reflected  three 

times  from  base  of  crust  85.2 


1A.9 

17.5 


PMAP 


P-wave  reflected  four 
times  from  base  of  crust 


88.7 


21.0 


PM2PpPp 


Smip 

and 

PM1S 


P-wave  reflected  twice 
from  base  of  crust  and 
continuing  as  a  head  wave  79.6 

S-wave  (or  P-wave) 

converted  on  reflection 

at  base  of  crust  to  a 

P-wave  (or  S-wave)  90.27 


11.9 


22.6 


Table  5-2 
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The  phase  SPP  theoretically  should  arrive  about  4.2  sec  after  the 
arrival  of  Pn.  See  cable  5-2.  This  is  also  the  approximate  arrival  time 
of  PPS.  The  PPS  head  wave  is  probably  of  lower  amplitude  than  SPP  but  it 
may  take  several  cycles  before  the  SPP  becomes  dominant  on  the  particle- 
motion  diagrams.  The  vertical  trace  of  STILLWATER  shows  a  linear  in¬ 
crease  in  amplitude  for  the  wave  group  in  the  interval  At  =  4.5  sec 
to  A  r  =  7.5  sec.  It  has  not  yet  been  determined  whether  the  linear 
increase  in  amplitude  is  associated  with  a  "beat"  or  is  a  buildup  of  a 
particular  phase,  or  is  the  combined  effect  of  two  or  more  arrivals,  or 
is  a  resonance  effect.  Some  evidence  in  favor  of  a  resonance  effect  is 
the  repetition  of  the  phenomenon  on  the  seismogram  several  times  in  the 
first  20  seconds.  The  linear  increase  in  amplitude  seems  to  be  common 
to  the  Pn  group,  the  SPP  group,  and  the  P^P  group,  suggesting  a  common 
station  effect  which  promotes  resonance  for  all  significant  arrivals, 
regardless  of  the  transmission  path. 

The  theoretical  travel-time  of  an  SPS  phase  based  on  the  model  of 
Figure  5-7a  is  76.1  sec,  or  8.4  seconds  after  the  arrival  of  Pn.  There 
is  some  indication  of  shear  on  the  particle-motion  diagram  at  the  ex¬ 
pected  arrival  time  for  SPS  and  it  is  possi1  le  that  this  phase  is  present. 
The  motion  is  predominantly  S-like  for  several  seconds  after  the  expected 
arrival  time  of  SPS.  However,  the  motion  is  not  clear-cut  and  it  is  not 
possible  to  estimate  the  time  of  onset  of  the  phase. 

The  dominant  prograde  motion  changes  to  a  single  cycle  of  retrograde 
motion  at  about  12  seconds  after  Pn.  Following  this  change  are  several 
cycles  of  well-developed  shear  motion  just  preceding  the  arrival  of  Pj^P. 


STILLWATER-WINV  (A=484  Km.) 
AGOUTI  — W  I NV  (A=475  Km.) 

STOAT  -WINV  (,A=475  Km.) 

a  CRUSTAL  SECTION  -  N.  CENTRAL  NEVADA 


ARMADILLO-FSAZ  (A  =479.2  Km.) 
ARMADILLO “S  FAZ  (A  =  577.3  Km.) 

ARM  ADI  LLO-WMAZ  (A  =  387.9  Km) 

b.  CRUSTAL  SECTION  —  NORTHERN  ARIZONA 


Figure  5-7.  Simplified  crustal  sections 

(a) 

(b) 


North  central  Nevada 
Northern  Arizona 


PjjjP  appears  to  become  well  developed  at  approximately  14  seconds  after 
the  arrival  of  Pn.  This  agrees  well  with  the  theoretical  arrival  time 


for  ?Ml 


P  but  it  would  probably  also  agree  with  the  approximate  travel 


time  of  a  guided  wave  in  the  crust,  usually  called  P.  This  phase  is  the 


most  distinctive  arrival  on  most  of  the  seismograms. 


The  S-like  motion  proceding  Pj^P  (or  P)  for  STILLWATER-W1NV  is  very 

well  developed  and  is  prominent  on  the  radial  seismometer  trace.  However, 

the  observed  arrival  time  of  12.0  sec  after  Pn  does  not  coincide  with  any 

arrival  which  could  be  expected  to  have  a  large  amplitude  for  the  crust* 1 

section  being  used.  The  S-like  motion  may  be  the  result  of  the  complex 

interaction  oi  crustal  layering,  the  arrival  of  low  amplitude  portions 

of  the  PjfjP  phase,  and  other  phases  such  as  SPS(?).  The  S-like  motion 

preceding  PM  P  for  STILLWATER-MNNV  is  not  as  well  developed  on  other 
1;1 

seismograms  for  other  explosions. 

The  portion  of  the  seismogram  after  the  arrival  of  PM^P  (or  P?)  out 
to  approximately  27  seconds  is  believed  to  be  due  mainly  to  multiple  re¬ 
flections  from  the  base  of  the  crust,  or  multiple  reflections  of  P.  This 
energy  may  also  include  Sj^P,  that  is  shear  energy  converted  to  P  on  re¬ 
flection  at  the  Mohorovicic  discontinuity.  No  attempt  has  been  made  to 
identify  the  onset  time  of  each  multiple  on  the  seismograms.  The  approxi¬ 
mate  expected  arrival  times  for  9m^P  (i=2,  3,  4)  have  been  noted  in  paren¬ 
theses  on  STILLWATER-WINV. 


There  is  w<  .1  developed  S-like  motion  appearing  on  the  particle-motion 
diagrams  at  approximately  22.5  sec  which  coincides  with  the  expected  arrival 


time  of  Pm^S. 
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Discussion  of  the  ARMADILLO  event  as  recorded 
at  WMAZ,  FSAZ  and  SFAZ 

Stations  WMAZ,  FSAZ  and  SFAZ  are  in  northern  Arizona  and  in  line 
with  the  Nevada  test  site.  Epicentral  distances  from  NTS  for  WMAZ, 

FSAZ  and  SFAZ  are  387.9,  479.2  and  577.3  km,  respectively.  Inspection 
of  Figure  5-8  shows  some  of  the  events  which  make  up  the  early  part  of 
the  seismogram  at  each  station.  Arrival  times  are  shown  as  measured 
from  the  arrival  or  Pn. 

At  WMAZ,  using  the  crur  al  section  of  Diment,  Stewart  and  Roller 
(1961),  the  first  good  indication  of  shear  on  the  particle-motion  diagrams 
is  at  4.4  sec  after  the  first  motion.  This  coincides  wel1  th  t.o  ex¬ 
pected  arrival  of  PK  a  converter;  head  •  ?e  from  the  base  ov  >  he  crust. 

It  is  also  apparent  from  the  particle-motion  diagrams  that  the  observed 
shear  motion  is  "riding"  on  P  energy  of  much  larger  amplitude  which  could 
be  the  SPP  phase. 

For  this  crustal  model,  SPS  and  P^P  should  both  arrive  at  approxi¬ 
mately  8.5  sec  after  Pn.  Inspection  of  the  particle-motion  diagrams 
around  this  time  shows  a  strong  P-arrival  at  about  7.7  sec.  Several 
cj'des  of  large  amplitude  P-motion  follow.  The  P-motion  is  prograde. 

This  P-motion  is  believed  to  be  due  to  the  arrival  of  P^P,  a  reflection 
from  the  base  of  the  crust,  or  P,  a  guided  wave  in  the  crust.  There  is 
also  some  indication  of  shear  at  about  8.0  sec  lasting  for  about  a  cycle 
which  may  indicate  the  SPS  arrival.  At  approximately  10  sec  there  are 
several  cycles  of  good  S-like  moi  ion.  This  motion  does  not  coincide 
with  the  expected  arrival  time  of  any  of  the  phases  in  Figure  5-8.  The 
motion  is  retrograde. 
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EVENT  P£P 


After  the  arrival  of  P) ,  the  seismogram  out  to  'pproximately 

25  sec,  at  least,  is  probably  due  to  multiple  reflections  from  the  base 
of  the  crust,  or  multiple  reflections  of  P. 

It  is  possible  that  some  interned iate  layering  would  provide  a 
velocity  structure  giving  a  travel-time  of  about  10.0  sec  for  a  phase 
that  could  be  called  SPS,  as  in  the  Berg,  Cook,  Narans,  and  Dolan  (1960) 
model  where  SPS  is  generated  from  the  base  of  the  9  km  layer. 

It  is  difficult  to  draw  positive  conclusions  from  a  study  of  a 
single  event  such  as  ARMADILLO  recorded  at  a  single  station  such  as  WMAZ. 
However,  shear  motion  i^  appart  it  on  the  particle-motion  diagrams  at  4,4 
sec  and  8  sec  which  does  coincide  approximately  with  the  expected  arrival 
time  for  PPS  and  SPS,  respectively.  Further,  except  for  the  shear  at 
10  sec,  these  are  the  only  intervals  of  shear  apparent  on  the  particle- 
motion  diagrams  out  to  about  10  sec.  The  SPS  motion  is  considerably  larger 
than  the  PPS  motion. 

There  is  S-like  motion  recorded  at  FSAZ  at  about  A. 2  sec  after  the 
arrival  of  Pn.  This  coincides  very  well  with  the  expected  arrival  time 
for  the  PFS  head  wave  from  the  base  of  the  crust  (see  Figure  5-7 b  and 
Figure  5-8.  There  are  no  intervals  of  shear  preceding  this  on  the  seismo¬ 
gram,  The  S-like  motion  lasts  for  approximately  one  second  and  is  followed 
(at  5.5  sec  after  Pn)  by  3  cycles  of  larger  amplitude  P-motion  lasting 
to  almost  8  sec  after  Pn  on  the  seismogram.  Most  of  this  P-motion  is 
probably  associated  with  the  SPP  phase.  At  approximately  8  sec  after  Pn, 
several  cycles  of  higher  frequency  S-like  motion  are  apparent  on  "he 
particle-motion  diagrams,  coinciding  very  well  with  the  theoretical 
arrival  time  of  SPS  from  the  base  of  the  crust.  The  theoretical  arrival 
time  of  Pj^P  (or  p)  at  FSAZ  from  the  Nevada  test  site  using  the  crustal 
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section  of  Diment,  Stewart,  and  Roller  (1961)  is  approximately  11.2  sec. 

This  coincides  almost  exactly  with  the  arrival  of  high- amplitude  prograde 
P-like  motion  (see  Figure  5-10), and  it  is  almost  certainly  due  to  the 
arrival  of  PM!**  (or  P) .  Following  this  arrival  are  the  usual  large  ampl- 
tudes  due  to  multiple  reflections  from  the  base  of  the  crust. 

Station  SFAZ  recording  the  ARMADILLO  event  shows  frequencies  consider¬ 
ably  higher  than  those  recorded  at  FSAZ  for  the  same  event  (see  Figure  5-9). 
This  is  probably  a  station  effect.  The  signal  appears  to  show  the  effects 
of  reverberations.  The  motion  is  P-like  with  a  small  minor  axis.  No  S-like 
motion  is  evident  until  approximately  6  seconds  after  Pn  when  several 
cycles  of  S-like  motion  appear  lasting  approximately  two  seconds.  This 
S-like  motion  does  not  coincide  with  any  expected  S  arrival-time  on 
Figure  5-8.  Ringing  may  have  a  severe  effect  on  body-wave  arrivals  at 
this  station  for  thir  event.  P-motion  resumes  at  approximately  8  seconds 
after  Pn.  Shear-motion  is  again  apparent  beginning  approximately  10 
seconds  after  Pn  and  continuing  approximately  2  seconds.  This  interval 
of  shear  coincides  well  in  arrival-time  with  a  similar  interval  recorded 
at  station  WMAZ  and  noted  above.  The  S-like  motion  is  not  as  well  develop¬ 
ed  on  the  seismogram  recorded  at  FSAZ  for  the  ARMADILLO  event,  for  this 
time  interval,  but  the  particle-motion  does  indicate  that  some  shear  is 
also  present  in  the  dominant  P-like  motion  associated  with  the  arr.  val  of 

PMiP- 

The  theoretical  arrival  time  of  Pj^P  at  SFAZ  for  the  ARMADILLO  event 
based  on  the  crustal  section  of  Diment,  Stewart,  and  Roller  is  approximate¬ 
ly  15  seconds  after  Pn.  This  time  on  the  record  (see  Figure  5-9)  does 
not  coincide  with  :he  arrival  of  a  large  amplitude  P-phase.  The  large 
amplitudes  usually  associated  with  Pj^P  (or  P)  are  not  seen  until  about 
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(’6:30:57.5)  S?S 


Figure  5-9.  CalComp  digltal-to-analog  playback  of  the  AKJiADILLO  event  as  recorded  at  WMAZ,  FSAZ,  and  SFAZ 


16  seconds  after  Pn.  The  arrival  of  Pj^P  (or  P)  would  then  seem  to  be 
approximately  one  second  late  at  SFAZ  for  the  ARMADILLO  event.  It  is  to 


be  noted,  however,  that  the  usual  high  amplitude  following  P^  P  for  at 
least  10  or  15  seconds  which  we  associate  with  multiple  reflections  from 
the  base  of  the  crust  are  not  present  either,  on  this  seismogram  (see 
Figure  5-9)..  Therefore,  evidently  less  P-energy  has  traveled  as  a  guided 
wave  for  this  event  as  compared  with  FSAZ  and  WMAZ,  and  the  actual  onset 
of  Pj^P  at  SFAZ  probably  occurs  about  a  second  before  the  larger  ampli¬ 
tudes  appear  at  16  seconds  after  the  ari ival  of  Pn  at  SFAZ.  The  station 
is  about  100  km  from  FSAZ  and  attenuation  in  the  crustal  wave  guide  will 

made  the  onset  of  P..  P  more  difficult  to  pick  for  SFAZ.  It  can  be  seen 

M1 

on  the  seismogram  for  SFAZ  in  Figure  5-9  that  the  signal  appears  to  be  modu¬ 
lated  by  a  low- frequency  v<>  'e  that  begins  at  the  expected  time  of  arrival 


CHAPTER  VI 


METHOD  OF  ANALYSIS  OF  DATA  RECORDED  AT  PERMANENT  SEISMOGRAPH  STATIONS 

Data  available  for  analysis  from  permanent  seismograph  stations 
Table  6-1  summarizes  the  available  seismic  data  from  permanent  seismo¬ 
graph  stations  that  were  used  in  the  present  investigation.  A  total  of 
159  events  are  listed.  The  methods  of  analysis  of  these  data  will  be  given 
in  a  later  section. 
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Table  6-1  Data  from  permanent  seismograph  stations  that  were 
digitized  in  the  present  investigation 

STATION  CODE: 


BCN 

-  Boulder  City,  Nevada 

FGU 

-  Flaming  Gorge,  Utah 

PCU 

-  Price,  Utah 

SLC 

-  Salt  Lake  City,  Utah 

DUG 

-  Dugway,  Utah  (world-wide  station) 
Heavy  Mass  Benioffs 

DUG  (Port) 

-  Dugway,  Utah 

Portable  Benioffs 

SPB 

-  Bolivia  -  (world-wide  station) 

LPS 

-  El  Salvador  -  (world-wide  station) 

SYMBOL  KEY 

*  -  Seismic  event. 

**  -  Event  digitized  at  95  minute  angle 

from  recorded  base  line. 

+  -  Event  originally  *  and  later  confirmed 

in  the  February,  1964  edition  of  "The 
Effects  of  Nuclear  Weapons". 

#  -  Event  has  IBM  computer  display. 


DATE 

STATION 

SECONDS 

DIGITIZED 

REMARKS 

9-19-57 

BCN 

10 

RAINIER 

10-8-58 

BCN 

10 

Confirmed  Blast 

1  16-59 

BCN 

10 

Central  California 

2-27-59 

BCN 

6 

Central  Utah 

3-16-59 

BCN 

11 

Central  Nevada 

3-23-59 

BCN 

10 

Central  Nevada 

4-1-59 

BCN 

10 

West  Nevada 

7-24-59 

BCN 

10 

N.  California 

10-13-59 

BCN 

10 

Central  Arizona 

10-17-61 

FGU 

10 

Central  Utah 

-58- 


SECONDS 


DATE 

STATION 

DIGITIZED 

5-12-62 

DUG 

15 

6-1-62 

DUG 

18 

6-6-62 

DUG 

15 

6- 13-62 

DUG 

13 

6-21-62 

DUG 

20 

6-30-62 

DUG 

15 

7-4-62 

SLC 

15 

7-6-62 

SLC 

10 

7-6-62 

DUG 

14 

**  7-6-62 

DUG 

14 

7-13-62 

DUG 

7 

**  7-13-62 

DUG 

15 

7-24-62 

DUG 

15 

7-27-62 

DUG 

16 

8-3-62 

PCU 

10 

8-3-62 

SLC 

10 

8-24-62 

DUG 

18 

9-6-62 

SLC 

11 

9-14-62 

DUG 

20 

9-16-62 

DUG 

20 

9-24-62 

DUG 

14 

9-28-62 

DUG 

1j 

9-29-62 

DUG 

18 

**  10-5-62 

DUG 

9 

10-12-62 

DUG 

20 

10-19-62 

DUG 

16 

10-19-62 

DUG 

12 

**  10-19-62 

DUG 

31 

10-24-62 

DUG 

12 

10-27-62 

DUG 

15 

11-12-62 

DUG 

12 

12-7-62 

DUG 

14 

12-11-62 

DUG 

15 

12-12-62 

DUG 

20 

12-29-62 

DUG 

15 

1-6-63 

DUG 

15 

1-27-63 

SLC 

15 

**  2-22-63 

DUG 

25 

**  2-22-63 

DUG 

28 

**  2-24-63 

DUG 

25 
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REMARKS 

AARDVARK 

#  RACOON  (twice) 

#  PACKRAT 

#  DESMOINES 

#  DAMAN  I 

#  SACRAMENTO 
Bingham  Blast 
Bingham  Blast 
SEDAN 

#  SEDAN 

MERRIMAC 

#  MERRIMAC 
Mexico 

#  WICHITA 
N.  Chile 

N.  Chile 

#  BOBAC 
Magna,  Utah 

#  HYRAX 
Mexico 

Hokkaido,  Japan 
W.  Columbia 

#  Confirmed  Blast 
MISSISSIPPI 

#  ROANOKE 

#  Confirmed  Blast 
Mexico 
Mexico 
Mexico 

#  Confirmed  Blast 

Andreanof  Island 

#  Confirmed  Blast 
Central  Utah 

#  NUMBAT 
N.  Chile 

Kurile  Islands 
Central  Utah 
North  Pole 
Dominican  Republic 
Guatemala 


SECONDS 


DATE 

STATION 

DIGITIZED 

REMARKS 

**  3-1-63 

DUG 

14 

#  * 

3-7-63 

DUG 

15 

South  Chile 

**  3-15-63 

DUG 

17 

#  * 

4-24-63 

PC. 

5 

Price,  Utah 

4-24-6.' 

SLC 

15 

Price,  Utah 

5-4-63 

DUG 

15 

Andreanof  Islands 

5-6-63 

PCU 

15 

Price,  Utah 

5-6-63 

DUG 

15 

Price,  Utah 

5-16-63 

PCU 

15 

Price,  Utah 

5-17-63 

DUG 

15 

Kurile  Islands 

5-22-63 

DUG 

24 

Mexico 

5-25-63 

DUG 

15 

Denver 

5-25-63 

PCU 

15 

Denver  (bad) 

5-27-63 

SLC 

15 

Fox  Island 

6-1-63 

DUG 

15 

Samoa 

6-3-63 

DUG 

18 

Honshu 

6-3-63 

DUG 

15 

Columbia 

**  6-5-63 

DUG 

12 

+  YUBA 

**  6-6-63 

DUG 

20 

#+  HUTIA 

6-10-63 

DUG 

15 

Kamchatka 

**  6-14-63 

DUG 

14 

+  MATACO 

6-24-63 

DUG 

15 

Fox  Island 

6-25-63 

SLC 

15 

Yellowstone 

6-25-63 

DUG 

15 

Yellowstone 

**  6-25-63 

DUG 

13 

+  KENNEBEC 

7-2-63 

DUG 

15 

Denver 

7-4-63 

SLC 

20 

Tonga  Islands 

7-4-63 

DUG 

20 

Tonga  Islands 

7-8-63 

DUG 

15 

N.  California 

7-15-63 

DUG 

15 

Kamchatka 

7-16-63 

DUG 

15 

Georgia,  USSR 

7-17-63 

DUG 

15 

West  Nevada 

**  7-20-63 

DUG 

12 

S .  Nevada 

8-8-63 

DUG 

15 

Fox  Island 

**  8-' 2-63 

DUG 

15 

PEKAN 

8-13-63 

DUG 

15 

Kodiak  Island 

9-8-63 

SLC 

20 

Fiji  Islands 

**  9-13-63 

LPS 

35 

#  BILBY 

**  9-13-63 

LPB 

37 

BILBY 

9-13-63 

SLC 

20 

BILBY 
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SECONDS 


DATE 

STATION 

’ TGITIZED 

REMARKS 

**  10-11-63 

DUG 

15 

+  GRUNION 

**  10-16-63 

PCU 

20 

+  CLEARWATER 

**  10-16-63 

SLC 

17 

+  CLEARWATER 

**  10-16-63 

DUG  (Port.) 

17 

+  CLEARWATER 

**  10-26-63 

DUG 

20 

#  SHOAL 

10-26-63 

DUG 

20 

SHOAL 

**  11-14-63 

DUG 

14 

ANCHOVY 

11-19-63 

DUG 

14 

Bingham  Blast 

11-19-63 

SLC 

20 

Bingham  Blast 

**  11-22-63 

DUG 

12 

+  GREYS 

**  12-6-63 

DUG 

18 

Mono  County,  Calif. 

**  12-10-63 

DUG 

35 

West  Bolivia 

**  12-12-63 

DUG 

13 

+  EAGLE 

**  12-14-63 

DUG 

37 

Central  Alaska 

**  12-16-63 

DUG 

26 

Ecuador 

**  12-19-63 

DUG 

32 

Argentina 

**  12-20-63 

DUG 

32 

N.  Peru 

**  12-21-63 

DUG 

45 

Tonga  Islands 

**  12-29-63 

DUG 

40 

N.  Chile 

**  12-30-63 

DUG 

34 

Kurile  Islands 

**  1-10-64 

DUG 

25 

Kurile  Islands 

**  1-16-64 

SLC 

20 

* 

**  1-19-64 

DUG 

25 

Solomon  Islands 

**  1-23-64 

DUG 

12 

* 

**  1-30-64 

DUG 

20 

* 

**  2-1-64 

DUG 

31 

Fox  Island 

**  2-6-64 

DUG 

30 

Columbia 

**  2-9-64 

DUG 

25 

Fiji  Islands 

**  2-13-64 

DUG 

25 

Bolivia-Brazil 

**  2-13-64 

DUG 

12 

★ 

**  2-19-64 

DUG 

29 

N.  Chile 

**  2-20-64 

DUG 

24 

Kurile  Islands 

**  2-22-64 

DUG 

30 

Kurile  Islands 

**  3-12-64 

DUG 

20 

* 

**  3-15-64 

DUG 

20 

#* 

**  3-23-64 

DUG 

20 

#  Walker  Lake,  Nevada 

**  3-24-64 

DUG 

20 

#  Walker  Lake,  Nevada 

**  3-31-64 

DUG 

25 

#  Kurile  Islands 

**  4-7-64 

DUG 

20 

#  Walker  Lake,  Nevada 

**  4-8-64 

DUG 

40 

Kurile  Islands 
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DATE 

**  4-9-64 
kk  4-14-64 
**  4-14-64 
**  5-14-64 
**  5-15-64 

**  6-11-64 
**  6-12-64 
**  6-18-64 
**  6-20-64 
**  6-25-64 

**  6-30-64 
**  7-7-64 
**  7-16-64 
**  7-19-64 
**  8-8-64 

**  8-10-64 
**  8-18-64 
**  8-24-64 
**  9-23-64 
**  1C-22-64 


STATION 

SECONDS 

DIGITIZED 

REMARKS 

DUG 

35 

El  Salvador 

DUG 

1A 

DUG 

33 

Kurile  Islands 

DUG 

13 

DUG 

15 

#* 

DUG 

20 

#* 

DUG 

20 

* 

DUG 

20 

★ 

DUG 

20 

Mexico 

DUG 

12 

★ 

DUG 

14 

★ 

DUG 

33 

Oregon  Coast 

DUG 

12 

DUG 

25 

#* 

DUG 

30 

Nicaragua 

DUG 

40 

S.  Peru 

DUG 

40 

Peru-Brazil 

DUG 

DUG 

30 

Montana 

19 

Las  Vegas,  Nevada 

DUG 

35 

SALMON 
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Data  processing 


Each  of  the  traces  of  the  vertical-component  and  two  horizontal- 
component  paper  seismograms  was  hand-digitized  for  a  total  length  of 
15  to  40  sec;  tne  length  of  the  digitizing  of  each  event  usually  de¬ 
pended  on  the  clarity  of  the  record  and  the  quality  of  amplitudes  of 
the  events.  The  digitizing  interval  was  0.05  sec.  The  digitizing 
was  done  with  a  Gaertner  microscope  xy-coordinate  comparator 
(Fig. 6-1  )•  The  readings  were  corrected  for  the  tilt  of  the  zero 
reference  line  of  the  seismogram  trace  by  drawing  a  straight  construc¬ 
tion  line  at  an  angle  of  9.5  minutes  of  arc  counter-clockwise  from  the 
recorded  base  line  at  the  beginning  of  the  time  interval  to  be  digi¬ 
tized.  For  a  recorder-drum  speed  of  60  mm/min,  the  error  in  amplitude 
involved  in  using  this  construction  line  as  the  corrected  zero  refer¬ 
ence  line  is  only  0.055  mm  and  0.110  mm  at  times  of  20  and  40  sec,  re¬ 
spectively,  after  the  start  of  digitizing;  and  this  error  in  amplitude 
was  considered  negligible  in  this  study. 

The  events  that  were  digitized  in  this  manner  are  given  in  table 
6-1  .  In  the  early  part  of  our  work  some  seismograms,  as  for  example 
those  from  the  Boulder  City  station,  were  digitized  at  0.1-sec  inter¬ 
vals;  and  also  some  of  these  were  not  corrected  for  tilt.  However, 
these  were  not  used  in  the  later  detailed  analysis  discussed  later  in 
this  section. 

All  paper  seismograms  that  were  digitized  were  carefully  checked 
for  differential  shrinkage  of  the  paper  by  measuring,  .with  the  Gaertr.er 
microscope  xy-coordinate  comparator,  the  distance  between  the  corxes- 
ponding  successive  minute  marks  on  the  Z,  E-W,  and  N-£  seismograms  in 
that  portion  of  each  seismogram  analyzed.  For  essentially  all  of  the 


-63- 


Figure  6-1.  Gaertner  microscope  xy-coordinate  comparator  used  for  hand 
digitizing  seismograms. 


seismograms  digitized  out  to  a  total  time  interval  of  about  20  sec 
(20  mm  for  a  recorder-drum  speed  of  60  mm/min) ,  the  differential 
paper  shrinkage  between  the  thiee-component  records  was  less  than  0.03 
sec  and  hence  considerably  less  than  one  digitizing  interval  (0.05  sec) 
t»ud  therefore  considered  negligible.  For  some  of  the  seismograms 
digitized  out  to  40  sec  (40  mm),  however,  the  differential  paper  shrink¬ 
age  slightly  exceeded  0.05  mm;  and  proper  adjustments  in  the  reading 
procedures  were  made  so  that  corresponding  points  on  the  three-component 
traces  could  be  read  and  therefore  compared  with  an  accuracy  of  one 
digitizing  interval  (0.05  sec). 

It  had  been  determined  from  graphical  syrtiiesis  that  an  error  in 
the  digitizing  interval  of  0.05  sec  for  seismic  waves  of  1.0-sec  period 
would  not  appreciably  change  the  character  of  the  particle-motion  dia¬ 
gram. 

The  digitized  data  for  the  Z,  E-W,  and  N-S  components  were  punched 
on  IBM  cards  for  subsequent  data  processing  on  either  the  IBM  7040  or 
IBM  1401,  or  both,  which  were  available  at  the  University  of  Utah  Com¬ 
puter  Center.  Various  computer  programs  were  used  in  the  presentation 
of  the  data  and  the  subsequent  analyses  of  the  data  that  will  now  be 
described. 

The  computer  also  computed  and  printed  out  the  values  of  epi- 
central  distances  (using  the  Richter  methoc  of  calculations  for  short 
distances  and  Turner's  method  (Macelwane  and  Sohon,  1932,  p.  285-287, 
for  teleseismic  distances)  and  printed  out  various  forms  of  data  pre¬ 
sentation  that  were  used  in  the  analysis  of  the  data. 


-65- 


Particle-motion  diagrams 


To  facilitate  the  compilation  of  the  particle-motion  diagrams,  the 
computer  computed  and  printed  out  (from  the  E-W  and  N-S  component  values) 
the  resolved  radial-horizontal  component  value  as  measured  along  a  hori¬ 
zontal  azimuth  toward  the  epicenter.  The  particle-motion  diagrams  were 
at  first  hand-plotted.  Later,  the  data  were,  processed  with  the  computer 
so  that  the  particle-motion  diagrams  were  plotted  automatically  on  a  Cal- 
Comp  plotter.  For  each  event,  particle-motion  diagrams  were  plotted  in 

(1)  a  vertical  plane  along  the  azimuth  toward  the  epicenter  and  (2)  a 
horizontal  plane;  in  this  report  these  will  be  designated  as  vertical  and 
horizontal  particle-motion  diagrams,  respectively. 

Other  presentations  of  data 

The  computer  printed  out,  in  the  form  of  a  closely  spaced  dot  pre¬ 
sentation  (at  0.05-sec  intervals)  on  the  abscissa  the  following  four  traces 
(Fig.  6-2): 

(1)  The  values  of  the  resolved  radial-horizontal  component  (Fig.  6-2D). 

(2)  The  values  of  the  vertical  component  (Z)  (Fig.  6-2E).. 

(3)  The  values  of  the  horizontal  north- south  component  (NS) 

(Fig.  6-2 F) . 

(4)  The  values  of  the  horizontal  east-west  component  (EW)  (Fig.  6-2G)- 
Figure  7-8  also  shows  an  example  of  this  type  of  presentation  for  the 
underground  nuclear  explosion  DAMAN  I,  in  which  trace  D  is  at  the  top  c 

the  diagram  and  traces  E,  F,  and  G  are  successively  below  it. 

To  facilitate  the  recognition  of  the  various  phases,  three  types  of 
closely  spaced  dot  presentations  were  printed  out  by  the  computer: 
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Figure  6-2.  Various  presentations  of  data  f or <  DAMAN  I  underground 
nuclear  explo  ion,  Nevada  test  site,  June  21,  1962:  dot  presen 
tions  to  help  iistinguish  character  of  ground  motion — (A)  S-  an 
P-like,  (B)  prograde  and  retrograde,  (C)  linear  and  elliptical; 
(D  to  E)  digital-to-analog  seismograms. 


A.  A  presentation  to  distinguish  between  P-like  and  S-like  motion. 

F igureb-3  shows  hew  this  presentation  is  obtained.  On  the  vertical 
particle-motion  diagram  (Fig.  6-3  A),  "P-like"  (points  a  and  b)  and 
"S-like"  (points  c  and  d)  motion,  respectively,  are  shown  between  two 
successive  points  digitized  at  0.05-sec  intervals.  The  absolute  magni¬ 
tudes  of  these  particle-motion  vectors  are  then  drawn  to  scale  (in  mm) 

as  ordinate  (Fig.  6-3  C)  so  that  P-like  and  S-like  motion  will  be  plotted 
above  and  below  the  reference  line,  respectively,  at  intervals  of  time 
corresponding  to  the  digitizing  interval.  A  predominance  of  the  length 
and  number  of  dots  above  the  abscissa  during  a  given  time  interval  indi¬ 
cates  P-Mke  motion  and  a  predominance  of  the  length  and  number  of  dots 
below  the  abscissa  indicates  S-like  motion.  Figure  6-2  A  gives  an  ex¬ 
ample  of  this  type  of  presentation  for  DAMAN  I. 

B.  A  presentation  to  distinguish  between  prograde  and  retrograde 
motion.  As  shown  on  the  vertical  particle-motion  diagram  (Fig.  6-3  A), 
the  value  of  the  angle  of  orientation  of  the  line  between  two  successive 
points,  as  measured  counterclockwise  from  the  horizontal  line  in  the 
direction  toward  the  epicenter,  is  either  0  to  Tf/2  or  ff  tc  3  7T/2 
radians  for  S-like  motion  and  either  7772  to  ff  or  3^/2  to  2Jf~ 

for  P-like  motion.  The  value  of  rf\.  is  drawn  to  scale  (in  radians)  as 
ordinate  (Fig.  6-3  D)  above  the  abscissa  at  intervals  of  time  correspond¬ 
ing  to  the  digitizing  interval.  This  presentation  indicates  the  direction 
of  apparent  rotation  of  the  motion  on  the  vertical  particle-motion  diagram: 
(1)  if  the  values  of  the  loci  of  the  top  of  the  dots  increase  as  time  in¬ 
creases,  the  motion  is  prograde  (counterclockwise);  and  (2)  if  the  values 
of  the  loci  of  the  top  of  the  dots  decrease  as  time  increases,  the  motion 
is  retrograde  (clockwise).  It  should  be  noted  that  forc(  =  0,  the 
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A.  VERTICAL  ARTICLE-MOTION  DIAGRAM 
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Figure  6-3.  Portion  of  (A)  vertical  particle-motion  diagram  and  (B)  horizontal  particle-motion  diagram  to 
show  method  of  compiling  dot  presentations  to  help  distinguish  character  of  ground  motion--(C)  S-  and 
P-like;  (D)  prograde  and  retrograde  (measurement  of  angle  ) ,  and  (E)  linear  and  elliptical. 


direction  of  motion  is  toward  the  epicenter.  Figure  6-2B  gives  an 
example  of  this  type  of  presentation  for  DAMAN  I. 

C.  A  presentation  to  distinguish  between  linear  and  elliptical 
motion.  On  the  horizontal  particle-motion  diagram  (Fig.  6-3  B) ,  the 
length  of  the  resultant  line  of  two  successive  points  (e  and  f)  is  pro¬ 
jected  on  both  (a)  the  horizontal  azimuthal  direction  away  from  the  epi¬ 
center,  for  which  the  component  value  is  plotted  to  scale  (in  mm)  above 
the  reference  line  (Fig.  6-3  E)  (this  component  is  the  resolved  radial- 
horizontal  component  value)  and  also  (b)  the  direction  perpendicular  to 
the  azimuthal  direc  ion,  for  which  the  component  value  is  plotted  to 
scale  (in  mm)  below  the  reference  line  (Fig.  6-3  E) .  In  this  presenta¬ 
tion,  a  purely  compressional  (P)  wave  or  shear  (SV)  wave  polarized  in 
a  vertical  plane  along  an  azimuthal  direction  from  the  epicenter  would 
have  values  in  the  direction  of  the  azimuth  only  (dots  above  the  abscissa 
only).  Thus  a  preponderance  of  dots  above  the  abscissa  indicate  linear- 
type  motion;  and  a  number  of  dots  below  the  abscissa  indicate  a  corres¬ 
ponding  amount  of  elliptical-type  motion  in  the  horizontal  plane.  Figure 
6-2c  gives  an  example  of  this  type  of  presentation  for  DAMAN  1. 

It  should  be  emphasized  that  the  motions  indicated  by  these  presenta- 
surface 

tions  are  ground /motions  only,  and  that  caution  must  be  used  when  relating 
this  motion  to  body  waves. 

In  Figure  6-2,  for  DAMAN  I  underground  nuclear  explosion,  these  dot 
presentations  are  displayed  together  with  the  digital-to-analog  seismo¬ 
grams  and  vertically  above  them  in  the  same  time  sequence.  This  type  of 
composite  presentation  affords  a  powerful  method  of  analysis  of  the  phases. 


-70- 


CHAPTER  VII 


INTERPRETATION  of  data  recorded  at  permanent  seismograph  stations . 

Consistency  of  the  character  of  body-wave  phases  of  underground  nu¬ 
clear  explosions  from  Nevada  test  site  and  recorded  at  Dugway, 

Utah 

Investigators  of  refraction  seismic  studies  of  earth  crustal  struc¬ 
ture  have  long  recognized  the  similarity  of  the  seismograms  obtained  at 
a  given  station  when  shots  are  repeated  at  the  same  or  near-by  shotpoints. 
For  example,  for  seismograms  observed  from  the  underground  nuclear  explo¬ 
sions  LOGAN  and  BLANCA  at  an  epicentral  distance  of  1610  km,  Romney  (1959, 
Fig.  2,  p.  1492)  recognized  the  similarity  of  the  waveforms  from  both 
shots;  he  reported  that  "each  of  the  HARDTACK  underground  nuclear  explo¬ 
sions  produced  almost  identical  waveforms  at  a  given  station".  To  our 
knowledge,  however,  little  quantitative  analysis  has  ever  been  made  of 
this  similarity. 

We  will  first  discuss  the  striking  consistency  of  the  body-wave 
phases  after  the  arrival  of  the  first  P  energy,  and  later  show  how  this 
consistency  can  be  used  for  the  positive  detection  and  identification  of 
certain  underground  nuclear  explosions  from  the  Nevada  test  site  as  re¬ 
corded  at  the  Dugway,  Utah  station. 

Available  data 

The  seismograms  used  in  this  study  were  obtained  at  the  Dugway 
seismograph  station,  Utah  (Fig.  7-1),  which  is  one  of  the  stations  in  the 
University  of  Utah  network.  The  instrumentation  is  the  type  used  in  the 
standardized  world-wide  network  operated  in  cooperation  with  the  U.S. 

Coast  and  Geodetic  Survey.  The  three-component  system  is  matched,  with 
a  magnification  of  400,000  for  each  of  the  three  seismometers,  namely 
the  vertical  (Z) ,  and  the  two  horizontal  (EW  and  NS)  seismometers.  The 
drum  speed  of  the  recorder  is  60  mm/min. 
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Figure  7-1.  Map  showing  Nevada  test  site  (NTS)  and  Dugway,  Utah  seismograph  station.  =  437  to  448  km. 
Bearing  of  NTS  from  Dugway  =  S.  40°  W. 


The  seismograms  for  a  total  of  22  announced  underground  nuclear  ex¬ 
plosions  at  the  Nevada  test  site  (NTS)  were  used  for  this  study.  Figure  7-2 
shows  the  locations  of  15  of  these  explosions,  for  which  the  latitudes 
and  longitudes  have  been  declassified.  The  announced  Richter  magnitudes 
of  these  events  were  3.9  to  5.8,  In  addition,  several  other  events  which 
have  not  yet  been  declassified  but  which  are  believed  to  be  underground 
nuclear  explosions  from  NTS,  were  analyzed;  however,  these  are  not  in¬ 
cluded  in  this  report. 

Tables  7-1  and  7-2  give  unclassified  pertinent  data  incident  to  the 
announced  nuclear  explosions  studies:  table  7-1  gives  principal  facts  that 
facilitate  the  intepretation;  and  table  7-2  includes  data  either  observed 
on  the  seismograms  or  computed  therefrom. 

Comparison  of  underground  nuclear  explosions 

r.n  the  present  study,  an  analysis  was  made  of  seismograms  obtained 
at  a  single  permanent  seismograph  station,  namely  Dugway,  Utah,  from  (1) 

22  announced  underground  nuclear  explosions  at  the  Nevada  test  site  (NTS) 
and  (2)  various  other  events  which,  on  the  basis  of  a  brief  examination 
of  the  original  seismograms  only,  were  considered  tentatively  to  be 
possible  explosions  from  the  NTS  and  therefore  worthy  of  detailed  study. 

During  this  analysis,  a  striking  consistency  in  the  character  of 
the  phases  of  the  body  waves  during  the  first  15  to  20  sec  after  the  first 
P  arrival  was  discovered  for  several  of  the  confirmed  underground  nuclear 
explosions  and  for  several  of  the  events  that  we  now  consider  to  be  posi¬ 
tively,  but  as  yet  unclassified,  underground  nuclear  explosions.  The  con- 
istency  of  the  c  -ter  of  the  phases  is  so  striking  that  no  only  can 
some  of  the  body  phases,  including  certain  provisionally  interpreted  con- 
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Figure  ?-2.  Topographic  map  (U.S.  Geological  Survey)  of  Yucca  Flat 
at  Nevada  test  site  and  the  names  and  locations  of  15  of  the  22 
announced  underground  nuclear  explosions  included  in  this  study. 


December  it,  19A1,  from  Lieut.  Robert  E  Rlecker,  t'.S  Air  Force,  who  obtained  this  information  from  published  data 


Date 

(CCT) 

Name 

Magnitude 

(Richter 

scale) 

1/ 

Yield 

(kt) 

Epi central 
distance 
to  Dcgway, 
Utah 
(km) 

Time  of 
delonat  ton 
(CCT) 

hr:mln: sec 

Latitude 

(North) 

Longitude 

(West) 

Type  of 
material 

Shot 

Depth 

(feet) 

5-12-62 

AARDVARK 

4.9 

33 

446.1 

19:00:00.103# 

37*03"55" 

116*01*49" 

Tuff 

1,434 

(1,424)** 

6-1-62 

racoon 

- 

Low 

448.0 

17:00:00.137# 

3  7  *02' 44" 

116*02*04" 

Alluvium 

550 

(539)** 

S-6-62 

PACES. 

4.4  t  0.5 

Low 

448.2 

17:00:00.122# 

37  *02 '45" 

116*02*21” 

Alluvium 

869 

(860)** 

6-13-62 

DES  MOINES 

4.4 

Low 

440.1 

21:00:00.116#*** 

37*13' 20" 

116*09'44" 

Tuff 

655 

(1,424)*J 
(1.205)  ~ 

6-21-62 

DAMAn  I 

Low 

448.0 

17:00:00.131# 

37*02  35" 

116*01*50" 

Alluvium 

860 

(854)** 

6-30-62 

Sacramento 

4.1  t  0.4 

Low 

5 

21:30:00.155# 

37  *07  ’03" 

116*02*51" 

Alluvium 

500 

7-6-62 

SEDAN 

4.9  !  0.4 

100 

437.3 

17:00:00.147 

37*10‘37" 

116*02*43*’ 

Alluvium 

600 

(635)** 

7-13-62 

MERR1MAC 

4.6  t  0.3 

Low 

447.1 

16:00:00.2# 

37*03 ’18" 

116*02*00" 

Alluvium 

1,370 

(1,356)** 

7-27-62 

WICHITA 

4.1  t  0.6 

Low 

442.0 

21:00:00.2# 

37*07  *47" 

116*03*23" 

Alluvium 

750 

(493)** 

8-24-62 
(2nd  expl.) 

BCBAC 

4.0  t  0.5 

Low 

447.4 

17:00:00.13# 

37  *02  ".0" 

116*01 *26" 

Alluvium 

680 

(675)** 

9-14-62 

HYRAX 

4.1  t  0.7 

Low 

441  .  . 

17:10:00.1# 

37*02*38.8 

T  16*01 *16.0" 

A! luvlua 

718 

(720)** 

9-29-62 

ALLEGHENY* 

- 

Low 

- 

17:00:00.2## 

- 

- 

Tuff** 

692** 

10-5-62 

MISSISSIPPI 

- 

lntermed** 

440.8 

17:00:00.15# 

37*08 ' 22" 

116*03*01" 

Tuff 

1,620** 

10-12-62 

ROANOKE 

3.9  t  0.3 

Low 

442.2 

15:00:00.155# 

37*07*22" 

116*03*03" 

Alluvium 

(Tuff) 

514 

(510)** 

10-19-62 

BANDICOOT* 

- 

Low 

- 

18:00:00.08#'* 

- 

- 

Alluvium** 

800** 

10-27-62 

SANTEE' 

- 

Low 

- 

15:00:00.14## 

- 

- 

Alluvium** 

1,050** 

12-7-62 

TENDRAC* 

- 

Low 

.  - 

19:00:00.09## 

- 

- 

- 

- 

12-12-62 

NUMBAT 

4.3  t  0.5 

Low 

446.9 

18:45:00.1*** 

37*02*46.0" 

116*00'56 . 2" 

Alluvium 

775 

6-6-63 

HUTIA 

- 

Low** 

- 

14:00** 

- 

- 

Al luvlim 

442** 

8-12-63 

PEKAN* 

- 

Low 

- 

23:45** 

- 

- 

Al luvlura** 

997** 

9-13-63 

BILBY 

5.8 

200 

446.0 

17:00:00.13 

37*03*38'* 

116*01*18" 

Tuff 

2,314 

(2,413)** 

10-26-63 

SHOAL 

- 

About  12** 

- 

17:00:00.100  -/ 

4/ 

39*12*01"  - 

118*22*49" 

Granite** 

1,205** 

11-14-63 

ANCHOVY* 

_ 

Low 

- 

16:00** 

- 

- 

Alluvium 

853** 

;  no 

fob 


i 


h 


) 


3 


5 


tow  yield  means  less  than  20  kllo'ona;  intermediate  means  20  to  999  kilotons,  inclusive.  All  yield  data  taken  from  Glass* 
if  Depression  is  subsidence  of  earth  into  underground  cavity,  as  distinguished  from  crater  by  throw-out  of  earth.  All  depress 


if  Excavation  experiment;  formed  crater  about  1,280  ft.  diameter,  320  ft.  maximum  depth  (Glasstone,  1964). 

4/  Written  comunicatlon  from  C.  C.  Bates  toT.  W,  Caless,  Oct.  14,  1964;  also  Caless,  1965,  p  1-1.  "ed 

#(1)  Written  communication  December  11,  1963  from  Lieut.  Robert  E.  Rlecker,  U.S.  Air  Force  and  (2)  written  communication,  July  i  C. 

March  29,  1963  from  Major  General  A.  W.  Betts,  U.S.  Atomic  Energy  Commission. 

nbe 

##  Written  communication,  July  9,  1962,  October  8,  196c  and  March  29,  1963  from  Major  General  A.  W.  Betts,  U.S.  Atomic  Energy  ‘*en 

*  Written  cosmunicatlon  from  Dixon  Stewart,  U.S.  Atomic  Energy  Consul  ssion,  received  October  21,  1964.  ■  9 

**  Glrsstone,  Samuel  (editor),  1964,  The  effects  of  nuclear  weapons:  published  by  U.S.  Atomic  Energy  Committee,  April  1962,  A 
February  1964.  Some  shot -depth  data  differ  somewhat  from  those  or. ginally  given  u»  bv  Lieut.  R.  E.  Rlecker. 

,  1 

***  The  detonation  times  of  2200  and  1745  for  DES  KOIKES  and  NUM8AT,  respectively,  on  p.  6776  of  Glasstore  are  incorrect  (writ:  -*P 

from  Brigadier  General  D.  L.  Crowson,  U.S.  Air  Fotce). 


:nd  nuclear  explosions  at  Nevada  test  site  (NTSl .  All  data,  unless  otherwise  nottd,  were  obtained  in  written  . umnunlcat ion, 
iobert  E.  Riccker,  U.S.  Air  Force,  who  obtained  this  information  from  published  data. 
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i/ 

Yield 

(kt) 

Eplcent  ral 
distance 
to  hugway, 
Utah 
(km) 

Time  of 
detention 
(GCT) 

hr :mln : sec 

Latitude 

(North) 

Longitude 

(West) 

Type  of 
material 

Shot 

Depth 

(feet) 

Surface 

Elevation 

(feet) 

Shot 

Elevation 

(feet) 

2/ 

D*pi  cation- 
(Dlaa  6 

Dept ' 

(fee.) 

38 

446.1 

19:00:00. 103d 

3 7*03 ’55" 

116*01 '49" 

Tuff 

1,434 

(1,424)** 

4.072 

2,638 

924;  !2 

Low 

448,0 

17:00:00.137# 

3 7 *02’ 44“ 

116*02 ’04" 

Alluvium 

5  SO 

(539)** 

4,021 

3,47! 

314;  24 

5 

Low 

448.2 

17:00:00.122# 

3 7 *02 ’45" 

116*02'21" 

Al ) uvium 

869 

(860)** 

4,022 

3,153 

598 ;4» 

Low 

440.1 

21:00:00.116#*** 

37*13' 20" 

116*09'44" 

Tuff 

655 

0.624)*#, 
(1,205)  ” 

5,652 

4,997 

None 

Low 

448.0 

17:00:00.131# 

37*02 ' 35" 

116*01 ' 50“ 

Alluvium 

860 

(854)** 

4,016 

3,156 

5 56; 96 

4 

Low 

442.5 

21:30:00.155# 

37*07 '03" 

11  J2'51" 

Alluvium 

500 

4,240 

3,740 

356; 41 

4 

100 

437.3 

17:00:0. 147 

37*10' 37" 

116 *02 '43" 

Alluvium 

600 

(635)** 

4,317 

3,657 

3/ 

) 

Low 

447.1 

16:00:00.2# 

37*03' 18" 

116*02 '00" 

Alluvium 

1,370 

(1,356)** 

4.341 

3,671 

662; 50 

Low 

442.0 

21:00:00.2# 

37*0, "47" 

116*03 '23" 

Alluvium 

750 

(493)** 

4,239 

3,489 

390:36 

> 

Low 

447.4 

17:00:00.13# 

37*02'46" 

116 *01 '26" 

Alluvium 

680 

(675)** 

4,013 

3,333 

425:45 

7 

Low 

447.4 

17:10:00. 1# 

37*02 “38.8 

116*01'16.0" 

Alluvium 

716 

(720)** 

4,016 

3,298 

474; 100 

Low 

- 

17:00:00.2## 

- 

Tuff** 

692** 

- 

100; 10 

Interned** 

440.8 

17:00:00.15# 

J 7*08 ' 22" 

116*03*01" 

Tuff 

1,620** 

2,609 

425 ; 125 

3 

Low 

442.2 

15:00:00.155# 

37*07 ’22" 

116  *03 '03” 

Alluvium 

(Tuff) 

514 

(510)** 

4,197 

3,683 

80 ;  5 

Low 

- 

18:00:00, 08*1 

- 

- 

Alluvium** 

800** 

- 

- 

300; 100 

Low 

- 

15:00:00. 14## 

• 

- 

Alluvium** 

1,050** 

- 

- 

400; 20 

Low 

.  - 

19:00:00.09## 

- 

- 

• 

- 

- 

- 

- 

5 

Low 

446.9 

18:45:00.1*** 

3 7 *02 ’46.0" 

116*00'56 . 2" 

Alluvium 

775 

4,030 

3,255 

500;  (>0 

Low** 

- 

14:00** 

- 

- 

Alluvium 

442** 

- 

- 

300; 20 

Low 

- 

23:45** 

- 

- 

Al luvlum** 

997** 

- 

- 

550; 60 

200 

446.0 

17:00:00. 13 

37*03'38" 

116*01' 18" 

Tuff 

2,314 

(2,413)** 

4,074 

1,760 

1800; CO 

About  1 2** 

- 

17:00:00.100  - 

.  .  v 

39*12  01"  “ 

4/ 

118*22'49"  ~ 

Granite** 

1,205** 

5,248  ~ 

- 

- 

Low 

_ 

16:00** 

. 

• 

A. luvlum 

853** 

. 

. 

600;  6  5 

20  kllotons;  intermediate  means  20  to  999  kilotons,  inclusive.  All  yield  data  taken  from  Glasstone,  1964. 

f  eavth  Into  underground  cavity,  as  distinguished  from  crater  by  throw-out  of  earth.  All  depression  data  taken  from  Glasstone,  1964. 


ied  crater  about  1,280  ft.  diameter,  320  ft.  maximum  depth  (Glasstone,  1964). 


'1.  C.  C.  Bates  to  T.  W  Caless,  Oct.  14,  1964;  also  Caless,  1968,  p  J-l. 

•nber  11,  1963  from  Lieut.  Robert  E,  Rlecker,  U.S.  Air  Force  and  (2)  written  communication,  July  9,  1962,  October  8,  1962,  and 
General  A.  W.  Betts,  U.S.  Atomic  Energy  Commission. 


j  9,  1962,  October  8,  1962  and  March  29,  1963  from  Major  General  A.  W.  Betts,  U.S.  Atomic  Energy  Commission. 
Dixon  Stewart,  J.S.  Atomic  Energy  Conrolsslon,  received  October  21,  1964, 


1964  The  effects  of  nuclear  weapons:  published  by  U.S.  Atomic  Energy  Committee,  April  1962,  Appendix  B; 
depth  data  differ  somewhat  from  thoseoriginally  given  us  by  Lieut  R.  E.  Pieckcr. 

“0  and  1745  for  DFS  MOINES  and  NUMBAT,  respectively,  on  p,  6776  of  Glass  re  are  incorrect  (written  eommu 
.  Crowson,  U.S.  Air  Force). 


revised  edition  reprinted 


icitlon,  January 


Table  7i2i —  Obaerved  travel  tlaea  and  other  pertinent  data  of  announced  underground  nuclear  exploaiona,  Nevada  teat  aitc  (NTS) 


Date 

(GCT) 

Nirae 

Epiceotral 

Distance 

to 

Dugway , 

Utah 

(taa) 

Time  of 
Detonat ion 
(GCT) 
h  m  5 

Timed  of 
Start  of 
Digitising 
(COT) 
h  m  s 

Observed 
(or  Estimated) 
Time^  of 
Arrive!  of  Pn 
(GCT) 
h  m  s 

Obaerved  (or 
Estimated) 
Travel  Time 
of  Pn 
(sec) 

Title  of  Pn(?) 

arrow  to 

Left  (L)  or 
Right  (*) 
of  G.00  on 
Digltal-to- 
Analog 
Seimnogram 
(aec) 

Length 
of  Tioe 
Digitized 
_ I»ec\_ 

10-26-63 

SHOAL 

489.8 

17  00  OO. 100 

17  01  08.8 

17  01  08.8 

68.  7 

0.0 

6-6-62 

PACKRAT 

448.2 

17  00  00.  122 

17  01  02.1 

(17  01  02.5) 

(62.4) 

0.4  to  R 

15 

6-21-62 

DAMAN  1 

448.0 

17  00  00.131 

17  01  02.1 

( 17  01  02.0) 

(61.9) 

0.1  toL 

20 

6-1-62 

RACCOON 

448.0 

17  00  00  137 

17  01  03.8 

* 

* 

* 

18 

8-24-62 

BOBAC 

447.4 

17  00  00.13 

17  01  02.1 

(17  01  02.3) 

(62.2) 

0.2  to  t 

18 

9-14-62 

’.mux 

447.4 

17  10  00.1 

17  11  02.1 

(17  11  02.5) 

(62.4) 

0.4  to  R 

20 

7-13-62 

MERR1MAC 

447.1 

16  00  00  2 

16  01  01.0 

(16  01  02.2) 

(62.0) 

1.2  to  R 

15 

12-12-62 

NUMBAT 

446.9 

18  45  00.1 

18  '<6  02.1 

(18  46  02.4) 

(62.3) 

0.3  to  R 

20 

5-12-62 

AARDVARK 

446.1 

19  00  00.103 

19  01  02.5 

19  01  02.4** 

62.3(7) 

0.1  to  L 

14 

9-1 <-63 

B1LBY 

446.0 

17  00  00.13 

17  01  02.2 

17  01  02.3 

62.2 

0.1  to  R 

6-30-62 

Sacramento 

442.5 

21  30  00.155 

21  31  02.0 

(21  31  01.8) 

(61.6) 

0.2  to  L 

15 

10-12-62 

ROANOKE 

442.2 

15  00  00.155 

15  01  02.1 

(15  01  01.8) 

(61.6) 

0.3  to  L 

20 

7-27-62 

WT  CHITA 

442.0 

21  00  00.2 

21  01  02.1 

(21  01  01.7) 

(61.5) 

0.4  to  L 

16 

10-5-62 

MISSISSIPPI 

440.  8 

It  00  00.15 

17  01  02.1 

(17  01  02.0) 

(61.9) 

0.1  to  L 

9 

7-6-62 

SEDAN 

437.3 

17  00  00.147 

17  01  01.1 

(17  01  01.3) 

61.2 

0.2  to  R 

14 

6-13-62 

DES  KOINES 

440.1 

21  00  00. U6 

21  01  02.1 

(21  01  01.2) 

(61.1) 

0.9  to  L 

13 

9-29-62 

ALLEGHENY 

17  00  00.2 

17  01  02.1 

(17  01  01.6) 

(61.4) 

0.5  to  L 

18 

10-19-62 

BANDICOOT 

18  00  00.08 

18  01  02.1 

(18  Oi  02.6) 

(62.5) 

0.5  to  R 

16 

10-27-62 

SANTEE 

15  00  00.14 

15  01  02.1 

* 

* 

a 

15 

12-7-62 

TENDJUC 

19  00  00.09 

19  01  03.6 

(19  01  02.6) 

(62.5) 

1.0  to  L 

14 

6-6-83 

HUTIA 

14  00 

14  01  03.8 

* 

* 

* 

20 

8-12-63 

PEKAN 

23  45 

23  46  02.3 

(23  46  02.3) 

1/ 

(62.2(7))- 

* 

15 

11-14-63 

ANCHOVY 

16  00 

16  01  05.3 

(16  01  02.4) 

<62. 30))-7 

2.9  to  L 

14 

*  Aaplitude  too  ame'l  for  accurate  determination. 

**  Cm  EV  seismogram  only  (vertical  seismogram  haa  a  considerable  amount  of  local  obiae.  may  have  arrlyad 
•lightly  earlier  than  indicated  on  this  aeisaogran,  which  alao  had  a  conalderable  amount  of  local  noise 

#  Corrected  for  any  neceasaiy  clock  correction. 

— ^  Aaauaea  that  the  time  of  detonation  vaa  about  90.1  aec  after  the  minute  indicated  in  tie  "Time  of  detonation”  column. 
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verted  wave?,  be  recognized,  but  also  certain  underground  nuclear  explo¬ 
sions  can  be  positively  detected  and  identified  by  this  similarity. 

In  the  present  report,  this  consistency  of  the  character  of  the  phases 
of  certain  underground  nuclear  explosions  will  be  described;  and  an  inter¬ 
pretation  of  the  recognizable  bcdy-wave  phases  and  a  comparison  of  their 
relative  amplitudes  will  be  given. 

Crustal  models  and  possible  body-wave  phases  for  area  of  investigation 

The  crustal  model  of  Berg  and  others  (1960)  for  the  eastern  part  of 
the  Basin  and  Range  province,  as  determined  from  chemical  explosions  at 
Promontory  and  Lakeside,  Utah  and  from  underground  nuclear  explosions  at 
the  Nevada  test  site,  comprises  (1)  a  top  9-km-thtck  layer  with  a  velocity 
of  5.73  km/sec;  (2)  a  16-km-thick  layer  with  a  velocity  of  6.33  km/sec; 
and  (3)  a  47-km- thick  layer  with  a  velocity  of  7.59  km/sec  overlying  a 
mantle  at  a  depth  of  72  km  with  a  velocity  of  7.97  km/sec.  The  crustal 
model  of  Pakiser  and  Hill  (1963),  which  was  obtained  with  an  unreversed 
northward- trending  seismic  refraction  profile  from  underground  nuclear 
explosions  at  the  Nevada  test  site  that  extended  northward  to  Elko, 

Nevada  and  beyond,  comprises  essentially  a  one-layer  crust  with  a  velo¬ 
city  of  6.03  km/sec  overlying  the  mantle  at  a  depth  of  28  km  and  with  a 
velocity  of  7.84  km/stc.  Pakiser  and  Hill  (1963,  p.  5763)  cite  definite 
evidence  for  an  intermediate  layer  of  velocity  6.7  km/sec  along  reversed 
profiles  extending  from  Elko,  Nevada. 


-77‘ 


The  epicentre!  distance  of  the  events  for  most  of  the  seismograms 

investigated  in  this  phase  of  the  study  was  437  to  448  km  (table  7-2). 

For  these  epicentral  distances  and  using  the  earth  crustal  model  of  Berg 

and  others  (i960)  for  the  eastern  part  of  the  Basin  and  Range  province, 

Figure  7-3  shows  the  ray  paths  and  some  of  the  various  phases  hat  could 

theoretically  be  expected  to  arrive  at  Dugway  from  NTS.  Listing  them  in 

the  order  of  expected  arrival  at  Dugway,  we  have: 

Pn  The  refraction  from  the  25-km-depth  boundary. 

(l. e. ,  PPPPP) 

PPPPPPP  The  refraction  from  the  72-km-boundary. 

and/or 

PjXjP  The  reflection  off  th°  72-km-depth  boundary. 

SSPPP  A  P  phase,  which  is  a  head  waves  which  originated 

at  the  explosion  as  shear  and  which  was  converted 
to  P  at  the  25-km-depth  boundary,  and  continued  on  as 
P. 

SSPSS  An  S  wave  which  is  converted  co  P  at  the  25-km-depth 

boundary  and  travels  back  up  through  the  crust  as  S. 

PPP  The  refraction  from  the  9-km-depth  boundary, 

and/or 

PjjP  The  reflection  off  the  25-km-depth  boundary. 

PPS  Head-wave  conversions  from  the  9-km-depth  boundary. 

and 

SPP 

SPS  A  shear  phase,  which  originated  at  the  explosion  as 

shear  and  which  was  converted  to  P  at  the  9-km-depth 
boundary,  and  then  converted  back  to  S. 

P  '  A  channel  wave  in  the  upper  part  of  the  crust. 
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MODEL  OF  BERG  AND  OTHERS  (I960) 


Figure  7-3.  Ray  paths  of  some  possible  phases  traveling  from  Nevada  test  site  to  Dugway,  ass 
layers  and  earth  crustal  model  of  Berg  and  others  (1960). 


Using  the  model  of  Berg  and  others  (I960),  Figure  7-4  shows  the  true 
angles  at  which  some  of  these  phases  should  arrive  at  Dugway  from  NTS. 

For  this  same  model,  Figure  7-5  shows  ti.e  travel-time  curves  for  these 
phases.  It  should  be  noted  that  for  an  epicentral  distance  of  about 
448  km,  for  example,  the  various  phases  arrive  at  intervals  of  about 
2  to  4  sec  after  each  other;  at  a  couple  of  times,  however,  there  is  a 
possible  overlapping  of  the  phases. 

Using  the  digitized  data,  both  vertical  particle-motion  diagrams 
(that  is,  in  a  vertical  plane  along  the  great-circle  azimuth  away  from 
the  explosion)  and  horizontal  particle-motion  diagrams  (that  is,  in  a 
horizontal  plane)  were  plotted--at  first  by  hand -plotting  and  later  in 
the  project  with  a  CalComp  plotter. 

Figure  7-6  shows  the  vertical  particle-motion  diagrams  for  DAMAN  I 
underground  nuclear  explosion.  The  scale,  in  mm  on  the  original  seismo¬ 
gram,  are  shown;  the  scale  may  change  from  one  diagram  to  the  next.  Also 
note  that  "toward"  and  "away"  from  the  explosion  are  to  the  right  and  left, 
respectively--which  is  opposite  to  the  usual  convention.  "Up"  and  "down" 
motion  are  as  indicated.  P-like  motion  would  be  dominantly  in  quadrants 
II  and  IV,  whereas  S-like  motion  would  be  dominantly  in  quadrants  I  and 
III.  The  dots  are  at  the  0.05-sec  digitizing  intervals,  and  the  numbers 
i re  in  seconds  after  the  start  of  digitizing. 
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Horizontal  and  Vertical  Scale 


Figure  7-4.  Ray  paths  of  some  possible  phases  arriving  at  Dugway  from  BILBY  underground  nuclear  explosion, 
Nevada  test  site,  September  13,  1963  assuming  horizontal  layers  and  earth  crustal  model  of  Berg  and  others 
(1960). 
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Figure  7-5.  Theoretical  travel-time  curves  for  various  refracted, 

reflected,  and  converted  waves  using  the  earth  crustal  model  of  Berg 
and  others  (1960)  for  the  eastern  part  of  the  Basin  and  Range 
province.  Horizontal  layering  is  assumed. 


Figure  7-6.  Vertical  particle-motion  diagrams,  DAMAN  I  underground  nuclear  explosion,  Nevada  test  site, 
June  21,  1962,  time  interval  0.00  to  15.35  sec;  0.00  =  17h  01m  02.1s  GCT. 


Using  the  same  crustal  model  of  Berg  and  others  (1960)  just  described, 
we  provisionally  interpret  the  phases  indicated  on  this  particle-motion 
diagram  as  follows: 

Approximate  time  Phase 

(sec) 

(Not  seen  on  this  record--ampiitude  Pn 

too  small) 


2.5 

4 

7  1/2 

8 

10.1 

11  1/2 

13 


PmP  and/  ?r  PPPPPPP 

SSPPP 

SSPSS 

PIIP  and/or  PPP 
PPS  and  SPP 
SPS 

P;  or  a  Multiple  of  P 


The  phases  correspond  well  in  arrival  time  with  those  shown  in 
table  7-3  for  this  model.  However,  there  are  still  some  problems  in 
identifying  positively  several  of  the  phases  that  need  to  be  resolved. 

Table  7.3  shows  the  theoretical  travel  times  of  early-arriving  body- 
wave  phases,  including  certain  converted  head  waves  and  certain  reflected 
waves,  for  an  epicentral  distance  of  448  km  (Nevada  test  site  to  Dugway, 
Utah)  for  the  crustal  models  of  Pakiser  and  Hill  (1963)  and  Berg  and 
others  (1960).  The  time  interval  /\t .  is  the  calculated  time  of  arrival 
of  each  phase  after  the  arrival  of  Pn  as  >  >ed  on  each  model. 

Although  the  travel  times  for  some  of  the  phases  differ  markedly  be¬ 
tween  the  model,  the  travel  times  for  other  of  the  phases  are  essentially 
the  same. 
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Table  7-3.  Theoretical  travel  times  of  early-arriving  body-wave 

phases,  including  certain  converted  head  waves  and  cer¬ 
tain  reflected  waves,  for  an  epicentral  distance  of  448 
km  (Nevada  test  site  to  Dugway,  Utah)  for  the  crustal 
models  of  Pakiser  and  Hill  (1963)  and  Berg  and  others 
(1960).  Time  is  the  calculated  time  oi  arrival  of 
each  phase  after  the  arrival  of  Pn,  as  based  on  the  model. 


Crustal  model 


Phase 

Pakiser  and  Hill  (1963) 

Berg  and  others  (1960) 

Travel 

Time 

(sec) 

At 

(sec) 

Travel 

Time 

(sec) 

At 

(sec) 

63.0 

0.0 

63.9  “ 

0.0 

PPPPPPP 

-- 

65.2 

1.3 

PIIIP 

-- 

65.3 

1.4 

SSPPP 

67.3 

4.3 

67.8 

3.9 

SSPSS 

71.6 

8.6 

71.7 

7.8 

PPP**  and/or  PI];P 

75.2 

12.4 

72.3 

8.4 

?PS  and/or  SPP 

-- 

74.0 

10.1 

SPS 

— 

-- 

75.7 

11.8 

Pg 

74.6 

11.6 

78.3 

14.4 

pns  or  SjjP 

81.3 

18.3 

77.3 

13.4 

Sn  (SSSSS) 

109.3 

46.5 

110.7 

46.8 

128.5 

65.7 

135.5 

71.6 

—  The  travel  time  of  the  first  arrival  shown  on  Fig.  7  of  Berg  and 
others  (1960,  p.  523)  is  62.9  sec,  which  is  in  excellent  agreement 
with  that  of  Pakiser  and  Hill  (1963). 

**  PPP  is  for  model  of  Berg  and  others  (1960)  only. 
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Figures  7-7  and  7-8  show  the  digital- to- analog  seismograms  of 
BT.LBY  and  DAMAN  I,  respectively.  In  each  figure,  from  top  to  bottom, 
are: 

D.  The  resolved  radial-horizontal  component 

E.  The  Z  component 

F.  The  north-south  component 

G.  The  east-west  component 
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Figure  7-8.  Digital-to-annlMg  seismogram,  DAMAN  I  underground  nuclear  explosion,  Nevada  test  site,  June  21 
1962;  0.00  *  17h  01m  02.1s  GCT. 


Observed  phases 


Arrival  of  Pn- -  Pn  (the  first  arriv. 1  of  P,  as  a  refraction  off  the 
25-km-depth  layer,  also  designated  as  PPPPP)  was  observed  at  Dugvay  with 
certainty  for  BILBY  only.  BILBY  gave  an  excellent  first  P-motion,  and 
its  observed  travel  time  of  62.2  sec  (table  7-2)  is  believed  accurate 
to  within  0.2  sec.  For  an  epicentral  distance  of  446.0  km,  this  value 
of  travel  time  compares  favorably  with  the  value  of  (1)  62.7  sec  taken 
from  the  travel-time  curve  of  Pakiser  ind  Hill  (1963)  for  eas ‘^rn  Nevada 
and  (2)  32.9  sec  taker,  from  the  travel-time  curve  of  Berg  and  others  (1960, 
p.  525),  provided  the  first-arrival  curve  on  Figure  7  of  Berg  and  others 
is  used.  It  should  be  noted  that  although  at  an  epicentral  distance  of 
448  km,  the  observed  travel-time  curve  of  Berg  and  others  (1960,  p.  524; 
Fig.  7,  p.  525)  for  the  underground  nuclear  explosion  BLANCA  indicates 
a  value  of  62.9  sec^  the  computed  travel  time  from  the  formula  for  their 
crustal  model  is  63.9  sec  because  their  formula  was  obtained  by  averaging 
data  other  than  that  from  BLANCA  alone.  This  1-sec  discrepancy,  however, 
is  reflected  in  the  travel-time  curves  of  the  crustal  model  of  Berg  and 
others  (1960)  in  Figure  7.5  by  a  1-sec  delay  in  the  absolute  time  of 
arrival  of  Pn.  Other  minor  variations  can  therefore  also  be  expected 
in  comparing  the  theoretical  (based  on  the  crustal  model)  and  observed 
values  of  the  time  intervals,  £,t,  after  the  first  arrival  of  Pn  (i.e., 
PPPPP  in  Figure  7-5  >. 


One  of  the  portable  stations  (9B)  used  in  the  measurements  from 
BLANCA  was  in  Skull  Valley,  Utah  within  about  20  miles  from  the 
Dugway  station. 
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The  observed  travel  times  of  the  first  arrivals  at  Dugwey  from  the 
underground  nuclear  explosions  SEDAN,  AARDVARK,  and  MISSISSIPPI  (table  7-1) 
give  partial  confirmation  for  the  observed  travel  time  for  BILBY,  provided 
their  different  epicentral  distances  are  taken  into  account.  It  is  found, 
however,  that  the  discrepancy  in  travel  times  between  BILBY  and  these  3 
explosions  may  be  as  large  as  half  a  second;  and  this  discrepancy  is 
attributable  largely  to  the  inaccuracy  of  picking  the  first  arrival 
amidst  the  microseismic  or  local  noise  on  the  seismograms  of  these  3 
explos ions , 

For  essentially  all  of  the  seismograms  used  in  this  study,  the 
arrival  of  Pn  was  not  definitely  observed  on  the  seismograms,  even  with 
our  microscope,  as  it  was  usually  masked  by  microseismic  noise.  However, 
to  study  the  consistency  of  the  body-wave  phases,  it  is  not  necessary  to 
observe  the  first  arrival  of  Pn  on  all  seismograms. 

Because  of  the  similarity  of  the  character  of  the  early  body-wave 
phases,  especially  as  observed  on  the  ditial-to-  ialog  seismograms  (as 
discussed  below),  the  positive  determination  of  the  arrival  of  Pn  on  the 
BILBY  digital -to-analog  seismograms  permitted  an  estimation  to  be  made 
of  the  probable  time  of  the  arrival  of  Pn  on  each  of  the  other  digitnl- 
to-analog  seismograms --even  though  the  Pn  arrival  could  not  actually  be 
observed  on  these  othe’  seismograms.  The  probable  time  of  the  arrival 
of  Pn  on  each  of  the  other  digital-to-analog  seismograms  was  determined 
by  (1)  matching  them,  over  a  light  table,  with  the  digital-to-analog 
seismograms  of  BILBY,  and  marking  on  them  the  position  where  the  arrival 
of  Pn  would  be  expected.  On  the  reproductions  in  this  paper,  this  posi¬ 
tion  is  indicated  with  "Pn  (?)"  and  a  vertical  arrow.  In  using  this 
procedure,  the  assumption  is  made  that  the  time  interval  between  the 
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arrival  of  Pn  and  the  very  early  body-wave  phases  chat  are  observed  (even 
when  Pn  is  not  observed)  does  not  vary  appreciably  :rom  seismogram  to 
seismogram  for  many  (but  not  necessarily  all)  underground  nuclear  explo¬ 
sions  that  occur  within  perhaps  10  km  of  each  other  within  NTS,  This 
assumption  is  believed  valid;  but  its  validity  is  still  being  tested. 

It  should  be  noted  that  the  accuracy  of  any  predicted  time  of  deto¬ 
nation  cf  an  event  that  is  detected  and  identified  as  an  underground 
nuclear  explosion  within  NTS  will  be  predicated  principally  on  the  accuracy 
of  the  determination  of  the  arrival  time  of  Pn  for  BILBY. 

Phases  after  Pn . - -  When  the  digital- to-analog  records,  which  are 
printed  out  by  the  IBM  1403  printer,  are  compared  by  placing  any  two 
records  over  each  other  on  a  light  table,  the  consistency  of  the  character 
of  the  body  phases  during  the  first  15  to  20  sec  after  the  arrival  of  Pn 
is  striking. 

During  approximately  the  first  2  sec  after  the  arrival  of  Pn  in  the 
seismograms  used  in  this  study,  little  consistency  is  generally  observed 
because  of  microseismic  background  noise;  if  any  is  seen,  one  must  inter¬ 
pret  it  with  caution.  Beginning  with  about  21/2  sec  and  thereafter  for 
most  of  the  digital-to-analog  traces  included  in  this  study,  the  amplitudes 
of  the  body  phases  are  generally  sufficiently  above  the  microseismic 
noise  level  so  that  consistencies  of  the  character  of  the  phases  can  be 
recognized  when  they  exist. 

and/  (table  7-4) 

PlIlP(?)/°r  PPPPPPP(-) . --  A  small  but  distinct  P  arrival/that  was 

observed  at  a  time  interval,  Z^t,  of  about  2  1/2  sec  after  Pn  from  essen¬ 
tially  all  of  the  nuclear  explosions  of  Richter  magnitude  4.1  or  greater 
is  provisionally  interpreted  as  a  reflection  or  refraction  from  the  72-km- 
depth  boundary  postulated  by  Berg  and  others  (1960)  in  the  eastern  part  of 
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the  Basin  and  Range  province.  The  distance  for  critical  reflection 
from  the  72-km  layer  of  the  model  of  Berg  and  others  (1960)  is  about 
350  km. 

SSPPPf?) . --  The  phase  tentatively  identified  as  SSPPP  is  a  promi¬ 
nent  phase  that  occurs  on  essentially  all  the  22  seismograms  recorded  at 
Dugway  from  underground  nuclear  explosions  detonated  at  NTS.  The  phase 
arrives  about  4  to  4  1/2  sec  after  the  arrival  of  Pn,  is  usually  charac¬ 
terized  by  a  W-shape  on  the  digital- to-analog  seismogram,  is  much  larger 
than  Pn,  and  is  prominent  even  when  Pn  is  too  small  to  be  identified. 

The  SSPPP  phase  is  designated  for  energy  which  apparently  originates  as 
a  shear  wave  at  NTS,  travels  down  through  the  two  upper  crustal  layers 
as  a  shear  wave  to  the  25-km-depth  boundary,  travels  as  a  P  head-wave 
along  this  boundary,  and  then  travels  up  through  the  two-layer  crust  as 
a  P  wave. 

The  phase  shows  about  equal  amplitudes  of  deflection  on  the  vertical 
(table  7-4) 

and  north-south  component^.  Particle-motion  diagrams  of  this  interval 
show  the  motion  to  be  P-like  with  an  apparent  incident  angle  (from  the 
vertical)  of  about  30°.  For  reasons  as  yet  unexplained,  the  east-west  com¬ 
ponent  shows  little  motion  although  uhe  energy  from  NTS  should  be  arriving 
at  about  equal  amplitude  on  both  horizontal  components. 

On  seismograms  with  energies  appreciably  above  the  microseism  level, 
a  very  symmetrical  envelope  can  be  drawn  around  the  peaks  of  the  SSPPP(?) 
phase  (see  DAMAN  I  and  MERRIMAC)  This  envelope  suggests  that  the  SSPPP(?) 
phase  may  start  at  a  time  interval  A.t  as  eirly  as  4.4  sec  (after  Pn) , 
which  is  close  to  the  theoretical  times  of  arrival  of  3.9  and  4.3  sec  for  the 
models  of  Berg  and  others  (1960)  and  Pakiser  and  Hill  (1963),  respectively. 
The  particle-motion  diagram  at  about  4  1/2  sec  shows  a  very  strong  dis- 
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tinctive  arrival  of  a  P-like  wave;  this  P-like  wave  is  perhaps  the  most 
striking  characteristic  found  when  comparing  all  records  and  is  probably 
the  start  of  the  SSPPP  phase. 

The  evidence  secrr.a  strong  that  the  phase  starting  at  about  2  1/2 
sec  is  not  the  start  of  the  high  amplitude  deflections  between  4  1/2 

and  7  sec  (the  SSPPP  phase).  The  apparent  gradual  build-up  of  the  north- 

south  component  suggests  that  the  opposite  may  be  possible.  In  most  cases 
the  amplitude  on  the  Z  component  decreases  after  2  1/2  sec,  indicating 
an  arrival  of  different  energy  at  4  1/2  sec,  although  a  destructively 
interfering  phase  causing  this  decrease  cannot  be  ruled  out.  However, 
the  particle-motion  diagrams  show  that  the  P  energy  at  2  1/2  sec  has 

a  much  steeper  angle  of  incidence  than  the  P  phase  at  4  1/2  sec;  and  on 

many  particle-motion  diagrams,  there  is  intervening  S-li!ce  motion. 

If  this  phase  is  correctly  identified  as  SSPPP,  the  large  amplitude 
of  this  phase  indicates  that  a  large  amount  of  shear  energy  must  be  gener 
ated  by  some  underground  nuclear  explosions.  The  large  amount  oi  shear 
energy  that  corresponds  approximately  with  the  theoretical  times  of  Sn 
(i,e.  SSSSS--see  below)  on  this  model  also  give  additional  supporting 
evidence.  If  this  presumption  of  a  large  amount  of  shear  energy  genera¬ 
ted  by  some  underground  nuclear  explosions  is  incorrect,  however,  the 
phase  designated  SSPPP  may  be  caused  by  a  channeling  effect  that  has 
not  ye*  been  recognized  to  date  by  the  refraction  work,  or  by  a  ringing 
effect  of  the  Dugway  station. 

SSPSS (?) . --  The  phase  designated  SSPSS(?)  is  an  S  wave  which  is  con 
verted  to  P  at  the  25-km-depth  boundary  and  travels  back  up  through  the 
crust  as  S.  Except  for  this  last  leg,  its  path  is  identical  to  that  of 
SSPPP(?).  The  phase  arrives  about  7  to  7  1/2  sec  after  Pn  and  is  much 
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(table  7-4) 

smaller  than  SSPPP(?y,  This  phase  apparently  corresponds  with  the  phase 

tentatively  identified  as  SPS(?)  by  Pakiser  and  Hill  (1963,  p.  5761),  who 

observed  the  phase  arriving  about  8  sec  after  Pn  and  by  Ryall  and  Stewart 

(1963,  p,  5821),  who  observed  the  phase  arriving  about  10  sec  after  Pn. 

PllP(?)  and/or  PPP(?).--  A  prominent  P  phase  arrives  consistently 

about  8  1/2  sec  after  Pn:  it,  too,  is  sometimes  characterized  by  a  W-shape 

on  the  d igi tal- to- analog  seismograms.  The  amplitude  of  the  deflection 

on  the  N-S  component  is  usually  about  twice  that  on  either  the  vertical 
(table  7-4) 

or  E-W  component/.  On  some  seismograms,  as  for  example,  on  the  N-S  com¬ 
ponent  of  DAMAN  I,  the  shape  of  the  phase  closely  resembles  that  of  the 
SSPPP(?)  phase.  Also  a  striking  comparison  was  founu  in  the  first  de¬ 
flection  of  both  phases,  if  present,  immediately  after  4  1/2  and  8  1/2  sec. 
The  first  peaks  of  these  phases  line  up  perfectly  on  nearly  all  records 
regardless  of  any  other  characteristics  on  the  record  or  the  locataion 
(within  NTS)  or  magnitude  of  the  explosion. 

On  some  traces,  there  is  some  evidence  of  "ringing".  A  consistent 
osciliacion  or  "ringing"  on  the  east-west  component  1  tarts  at  7.6  to  7.8 
sec  after  Pn,  and  gradually  increases  in  amplitude  to  persist  for  approxi¬ 
mately  4  sec.  This  appears  on  13  of  the  17  records  in  Area  I  and  3  of 
the  6  in  Area  II.  At  9  sec,  this  component  is  sometimes  still  large  enough 
to  influence  the  particle-motion  diagram  and  may  explain  the  inconsistent 
analysis  of  the  tvpe  of  motion  (S-like  or  P-like)  due  to  this  phase  as 

observed  on  the  vertical  and  north-south  components. 

—  / 

This  phase  is  apparently  the  same  as  P  (after  Mohorovicic)  observed 
by  Ryall  and  Stuart  (1963,  p.  5821)  and  interpreted  by  them  as  a  phase 
"associated  with  waves  reflected  beyond  the  critical  angle,  from  the  base 
of  the  crust"  In  our  notation,  this  phase  is  designated  PnP. 
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head-wave  conversions  from  the  9-km-depth  layer  tnat  arrive  about  10  to 
10  1/2  sec  after  Pn.  On  some  particle-motion  diagrams,  good  S-like  mo¬ 
tion  is  observed;  on  others,  elliptical-type  S  motion  is  observed;  and 
on  still  others,  inconsistencies  are  observed  that  suggest  that  these  two 
phases  may  be  intermixed.  The  amplitudes  of  these  phases  arc  generally 
much  smaller  than  P££p  and/or  PPP. 

SPS(?) . --  A  prominent  phase,  tentatively  identified  as  SPS,  occurs 

(table  7-4) 

at  about  11  1/2  sec  after  P  /;  it  usually  comprises  2  cycles  on  the  record. 
This  phase  was  observed  in  essentially  all  the  records  studied.  On  the 
particle-motion  diagrams,  the  S-like  motion  is  usually  distinct  and  large.” 

The  validity  of  our  interpretation  of  the  PPP,  PPS,  SPP,  and  SPS 
phases  rests  on  the  existence  of  an  intermediate  layer  at  a  depth  of  about 
9  km,  which  was  not  found  conclusively  by  Pakiser  and  Hill  (1963)  beneath 
the  Nevada  test  site. 

P';or  Multiple  P.  -  -  A  large-amplitude  phase  arrives  consistently  about 
(table  7-4) 

13  sec  after  This  phase  may  be  (1)  a  multiple  reflection  of  P  from 

the  25-km-depth  boundary  (which  goes  back  to  the  surface,  then  back  to  the 
same  boundary  and  finally  back  to  the  surface)  or  (2)  a  guided  or  trapped 
wave  in  the  upper  part  of  the  crust.  The  theoretical  arrival  times  for 
such  a  multiple  reflection  is  10.4  and  14. 5(?)  sec  for  the  crustal  models 
of  Berg  and  others  (1960)  and  Pakiser  and  Hill  (1963^  respectively. 

Sn  'i.e.,  SSSSS).--  At  the  time  interval  ^t  (after  Pn)  of  50  to  53 
sec,  large  amplitudes  (table7-4).  on  the  horizontal  comoonents  are  ob¬ 
served.  This  time  is  approximately  the  theoretical  time  for  a  body  Sn  (i.e 
SSSSS)  wave.  T^ese. large  amplitudes  give  strong  support  to  the  hypothesis 
that  a  large  amount  of  shear  energy  is. generated  in  most  of  tne  underground 
nuclear  explosions  studied.  However,  particle-motion  diagrams  were  not 
made  on  these  seismograms  at  this  time  interval. 

17 

It  is  difficult  to  explain  why  the  amplitude  of  SPS(?)  is  proportionally 
so  much  larger  than  that  of  SPP(?),  -95- 


Consistency  of  patterns 


Figures  7-8,  7-9,  and  7-lDshow  the  digital-to-anaiog  seismograms  for 

8  announced  underground  nuclear  explosions  within  NTS.  On  each  figure, 
the  4  letters  "a,  b,  c,  and  d"  are  placed  on  key  peaks  or  troughs  that 
occur  consistently  on  all  the  seismograms  shown;  and  these  key  peaks  and 
troughs  can  be  used  to  match  the  seismogram  of  one  explosion  against 
another.  Once  they  are  matched  (over  a  light  table,  for  example),  there 
are  literally  tens  or  scores  of  other  consistencies  that  can  be  found 
between  these  underground  explosions. 

The  first  explosion,  DAMAN  I,  is  known  to  be  located  within  1  1/2 
km  of  the  second,  MERRIMAC,  as  their  latitudes  and  longitudes  have  been 
declassified.  The  character  of  the  seismogram  of  the  third  explosion  is 
so  consistent  with  that  of  the  other  two  that  it  was  probably  detonated 
in  the  same  general  area  within  NTJ  as  the  other  two. 

It  should  be  emphasized  that  certain  other  announced  underground 

(area  I) 

nuclear  explosions  detonated  in  this  same  area/within  NTS  as  the  previous 
3  examples  do  not  show  the  excellent  consistency  shown  by  these  3  examples. 
The  main  point  is,  however,  that  if  we  do  observe  a  seismogram  essentially 
identical  to  these  three,  we  believe  we  can  correctly  detect  and  identify 
it  as  positively  an  underground  nuclear  explosion  from  this  area  within 
NTS. 

Similar  comparisons  have  been  made  on  other  seismic  events  which  are 
believed  to  be  underground  nuclear  explosions  from  the  Nevada  test  site 
but  which  have  not  yet  been  declassified. 

Figure  7-11  shows  the  digital- to-analog  seismogram  for  SEDAN,  which 

(area  II) 

was  detonated  in  another  area/within  NTS.  In  detail,  its  character  is 
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Figure  7-10.  Digital-to-analog  seismogram,  announced  underg  nd  nuclear  explosion,  Nevada  test  site 


Figure  7—  LI .  Digital-to-analog  seismogram,  SEDAN  underground  nuclear  explosion,  Nevada  test  site,  July 
1962;  0.00  =  17h  01m  Oi.ls  GCT. 


ever,  some  of  the  corresponding  phases  can  be  recognized. 

Figure  7-12  shows  the  digital- t^-analog  seismogram  for  DES  MOINES, 

iarea  III) 

which  was  detonated  in  yet  a  third  area/within  NTS.  Its  character  is 
markedly  different  from  any  of  the  preceding  seismograms.  For  example, 
it  has  a  ir  ich  higher  frequency  content.  Also  some  ringing  may  be  present, 

Positive  detection  and  identification  of  certain  underground 

nuclear  explosions  from  Nevada  test  site 

From  oui  study  we  conclude  as  follows; 

Provided  the  explosions  occurred  in  the  same  area,  the  various  pre¬ 
sentations  indicate  a  striking  consistency  of  the  body-wave  phases  from 
some  (but  not  universally  from  all)  underground  explosions  within  that 
area.  The  characteristics  of  the  body-wave  patterns  vary  from  three 
different  areas  of  detonation  within  the  Mercury  test  site,  presumably 
because  of  different  source  functions  (depending  on  rock  type,  geologic 
environment,  degree  of  coupling,  yield,  etc.).  This  striking  consistency 
has  been  used  successfully  for  the  positive  detection  and  identification 
of  certain  underground  nuclear  explosions  detonated  at  the  Nevada  test 
site  and  recorded  at  Dugway,  Utah  and  to  correctly  determine,  within  a 
few  kilometers,  the  area  of  the  detonation  within  NTS;  and  it  is  believ<  * 
that  the  times  of  detonation  of  these  certain  explosions  can  be  correctly 
predicted  within  one  second  usually  and  within  a  few  tenths  of  a  second 
sometimes . 

We  doubt  that  an  earthquake  within  the  NTS  area  could  give  the  ex¬ 
actly  identical  character  as  that  from  DAMAN  I  or  MERRIMAC,  for  example, 
but  we  have  no  data  from  earthquaxes  within  NTS  to  prove  it. 
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Figure  7-12.  Digital -to-analog  seismogram,  DES  MOINES  underground  nuclear  explosion,  Nevada  test  site 
June  13,  1962;  0.00  =  21h  01m  02.1s  GCT. 


Semipalatinsk  area,  USSR,  and  recorded  at  Dugway,  Utah 

We  next  will  compare  seismograms  from  4  probable  underground  explo¬ 
sions  in  the  Semipalati nsk  area,  USSR,  again  as  recorded  at  Dugway,  Utah. 
The  Rich!  r  magnitudes  of  the  events,  as  reported  by  the  U,S.  Coast  and 
Geodetic  Survey  preliminary  epicenter  cards,  were  5.6  or  greater:  and  the 
epicentral  distance  was  about  9,800  km  (about  90°). 

The  4  Soviet  probable  underground  explosions  have  many  striking  con¬ 
sistencies.  In  particular,  the  events  of  May  16,  1964  and  July  19,  1964 
are  almost  identical  in  many  respects.  As  given  by  the  U.  3.  Coast  and 
Geodetic  Survey-( table  7-5),  their  latitudes  are  identically  the  same  and 
their  longitudes  are  the  same  within  0.2°.  On  Figures  7-13  and  7-14,  which 
show  the  digital- to-analog  seismograms  of  the  2  explosions,  the  4  letters 
"a,  b,  c,  and  d"  are  placed  on  key  peaks  or  troughs  that  occur  consistently 
on  all  the  seismograms  shown;  and  these  key  peaks  and  troughs  can  be  used 
to  match  the  seismogram  of  one  explosion  against  another.  Once  they  are 
matched  (over  a  light  table  for  example,)  there  are  literally  tens  or 
scores  of  other  consistencies  that  can  be  noted  between  the  seismograms. 

The  excellent  consistency  of  these  2  seismograms  indicates  that  the 
events  are  from  the  same  area  and  have  essentially  the  same  source  func¬ 
tions.  We  doubt  that  any  earthquake  in  this  area  would  give  a  pattern 
as  nearly  identical  to  these  as  they  are  to  each  other. 

The  seismograms  of  the  other  two  Soviet  probable  underground  explo¬ 
sions  are  also  very  similar  to  these  two;  in  particular,  essentially  all 
the  major  characteristics  are  present.  The  slight  differences  in  character 

—  The  U.  S.  Coast  and  Geodetic  Survey  preliminary  epicenter 
determination  cards  show  all  4  of  these  Soviet  events  as 
"zero  depth"  events. 
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Table  7*5,--  Pertinent  data  incident  to  probable  large  underground  explosions  in  the 
Semipalatinsk  area,  USSR 

Epicentral 

distance 


Date 

(GCT) 

1/ 

Magnitude 

(Richter 

scale) 

to  Dugway 
Utah 

(degrees) 

(km) 

Origin  Time 
(GCT) 
h  m  s 

1/ 

1/ 

Latitude  ~ 
(north) 

Longitude 

(east) 

March  15 

,  1964 

5.6 

90°  00' 
(10,000) 

07 

59 

58.0 

49.7 

78.0 

May  16  , 

1964 

5.6 

89°  46' 

(  9,974) 

06 

00 

58.1 

49.9 

78.3 

July  19, 

1964 

5.5 

89°  48' 
9,977 

05 

59 

58.9 

49.9 

78.1 

Nov.  16, 

1964 

6.0 

90"  00' 
(10,000) 

05 

59 

57.4 

49.7 

78.0 

1/ 


U.S.  Coast  and  Geodetic  Survey  preliminary  epicenter  determination  cards,  which 


indicate  these  events  as  zero-depth  events. 


Figure  7-14.  Digital-to-analog  seismogram  of  probable  underground  explosion  in  the  Semipalatinsk  area, 
USSR,  July  19,  1964;  0.00  =  Q6h  13m  01.4s  GCT. 


be  from  slightly  different  areas. 


PS  converted  waves  from  earthquakes  and  large  underground  explosions 

at  epicentral  distances  of  2,000  to  10,000  kilometers 

In  this  section  wo  will  demonstrate  that  PS  converted  waves  have  been 
observed  from  earthquakes  and  large  underground  explosions  at  epicentral 
distances  of  2,000  to  10,000  kilometers.  In  Going  this,  we  will  also  be 
able  to  explain  certain  events  on  the  seismogram  that  contribute  to  the 
complexity  of  the  coda  of  P.  The  study  will  be  restricted  to  events 
from  the  teleseismic  distances  just  given. 

Available  data 

All  the  seismograms  to  be  discussed  were  recorded  by  the  Benioff 
short-period  matched  seismometers  (magnification  400,000)  at  Dugway, 

Utah.  The  seismograms  were  hand-digitized  and  analyzed  in  a  manner 
identical  to  that  discussed  in  the  previous  sections. 

The  seismograms  of  10  earthquakes  were  studies:  4  from  the  Kurile 
Islands,  2  from  Japan,  one  from  the  Tonga  Islands,  and  3  from  Mexico. 

(See  table  7-6  )  In  addition,  four  Soviet  probable  underground  nuclear 
explosions  from  the  Semipalatinsk  area  were  studied. 

Expected  arrival  time  of  PS  converted  waves 

The  time  of  arrival  of  a  PS  converted  wave  will  generally  be  ex¬ 
pressed  in  terms  of  the  time  interval,  ^>t,  after  the  first  arrival  of 
P.  To  confirm  the  travel  time  of  a  given  body-wave  phase,  however,  the 
Jef f reys-Bullen  tables  of  1958  have  been  used. 
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Table  7-7,  which  was  compiled  from  data  given  in  the  Jeffreys- 
Bullen  tables  of  1958,  gives  thi  time  interval,  At  (in  sec),  after  the 
initial  P,  at  which  various  body  phases  will  arrive  from  epicentral 
distance.  A?  of  iU°  (1,111  km)  to  90“  (lO.uuO  km).  The  phases  studied 
are : 

P  -  the  first  arrival 
pP 
sP 
PcP 

The  tables  include  the  data  for  different  depths  of  focus: 

(1)  zero  depth;  (2)  33-km  depth;  (3)  96-km  depth;  and  (4)  159-knt  depth. 

Figure  7-15  shows  the  earth  crustal  model  of  Berg  and  others  (1960), 
the  Dugway  station,  and  the  approximate  ray  paths  for  the  two  extreme 
ranges  of  epicentral  distances  of  the  events  studied.  The  angles  of  in¬ 
cidence  assumed  at  the  base  of  the  crust  were  computed  with  a  formula  by 
Gutenberg  (Richter,  1958).  Also  shown  on  this  diagram  are  the  P  to  S 
conversions  at  the  25-  and  72-km-depth  boundaries. 

Using  the  assumptions  of  this  model,  and  assuming  further  that  a 
PS  converted  wave  is  generated  at  the  postulated  72-km-depth  boundary  of 
Berg  and  others  (1960),  the  time  interval.  At.  of  the  arrival  of  such  a 
PS  converted  wave  after  the  first  arrival  of  P  would  be  about  8.7  sec  from 
a  shallow  earthquake  in  Mexico  and  about  7.7  sec  from  an  underground  ex¬ 
plosion  in  central  USSR, 

Arrivals  of  shear  (SV)  waves  that  are  observed  at  these  time  inter¬ 
vals  after  P  in  this  study  are  identified  as  PS  converted  waves  genera¬ 
ted  at  the  72-km-depth  boundary. 

Figure  7-16  shows  a  graph  of  the  time  interval.  At,  after  the 
initial  P,  at  which  a  PS  converted  wave  should  arrive  from  an  assumed 
one-layer  earth  crustal  model  for  various  assumed  velocities.  The  assumed 
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Epicenter 

Date 

J 

Latitude 

(degrees) 

J 

Longitude 

(degrees) 

Depth 

of 

focus 

(km) 

Origin 
time 
(GCT) 
h  m  r 

Observed 
t ime  of 
first  P 
Arrival 
(GCT' 
h  m  s 

T  ifee  of 
start  of 
digit  ir  :.*»g 

<ccn 

h  ffi  » 

EpleeM  ra  1 
distance 
(degrees) 
(ks) 

Rithter 
*agni t  ude 

Direct! 

A 

Kurile 

Islands 

Feb.  20,  1964 

44. 6N 

150. OE 

50 

09  53  51.1 

10  04  43. 7 

10  04  43,6 

67.3* 

(7,480) 

5.2 

(C.C.S.) 

-,61s*  -6 

Kurile  May  17,  1963  45, 3N 

Islands 


160. 8E  33 


Oi  06  36  04  17  26,7  04  17  26.6  66,6* 

(7,400) 


5  9  -.617  ♦ 

(f.C.S.) 


Kurile  Feb.  22,  1964  48, 5N 

Islands 


154, 9F,  60  17  50  56.2  18  01  16,6  18  01  16.6  62.6*  5.0  -.604  * 

(6.960)  (C.C.S.) 


Kurile  Jsn.  6,  1963  47. 4N 

Islands  _ 


155. «  33  21  20  56.5  21  31  19,8  21  31  19.8  62,6*  none  -.620  + 

_ _ (6 .960) _ given 


Hokkaido,  Sep:.  24,  1962  42. 8N 

Japan 


145. 3E  23 


14  38  21  14  49  40  5 


14  49  40.2  71.4* 

(7,930) 


none 

given 


..605 


Honshu,  June  3,  1963  34. 2N 

Japan 


138. 7E  43 


07  35  54  07  48  06.3  07  48  06.3  80.9* 

(8,990) 


5.3  -.624  +.V. 

(C.C.S.) 


Tonga  July  4,  1963  26.3*S  177. 7*H  158  10  58  13  11  10  54.8  11  10  54.5  69.6*  6  3/4 

Islands  (9,990)  (p.s)  -.897 

6  3/4-7 
(Ska) 

6) 

_ _ 

Mexico  Oct.  19,  1962  19,8*W  108. 3*W  53  21  21  48  21  26  26.4  21  26  75.3  70.8*  none  -.296  -.8 

(2,310)  given 


Mexico  May  22,  1963  17. 7N  106. ON  33  04  24  29  04  29  35.5  04  29  35.6  23,3*  4.6  -.263  -.9 

(2,590) 


Mexico  June  20,  1964  18.5*N  105. 5M  28  17  12  15  17  17  14.4  17  17  14,0  22.7*  5.5  -.754  ..9 

(7,520)  (C.C.S.) 


*  Tka  Picking  of  in  Individual  phase  1  laced  te  4J  fee  dona  mot  imply  that  the  Oneet  od  that  phaaa  can  8a  identified  oa  (ha  seiemegrea  to  vi * 
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Data  taken  from  U.S.  Coaat  and  Geodetic  Survey  preliminary  epicenter  cards. 
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Figure  7-15.  Ray  paths  of  P  waves  and  PS  converted  waves  penetrating  the  assumed  earth  crustal  model  of 
Berg  and  others  (I960)  at  angles  of  incidence  for  epicentral  distances  ( A )  appropriate  for  events  from 
Mexico  (left)  and  USSR  (right)  studied.  At  *  time  interval  between  PS  converted  wave  and  P. 


Figure  7-16.  Relationship  of  time  interval.  At,  between  PS  converted  wave  and  P  wave  for  an  assumed  one- 
layer  crust  with  different  thicknesses,  h,  and  for  various  angles  of  incidence,  L  ,  at  the  base  of  the 
crust.  Data  are  for  two  sets  of  compressional-velocity  values.  c 


velocities  are  those  included  in  the  earth  crustal  models  of  (1)  Berg 
and  others  (1960)  and  (2)  Pakiser  and  Hill  (1963).  To  obtain  the  approxi¬ 
mate  range  of  values  of  At  involved,  various  thicknesses  of  the  crust 
and  various  angles  of  incidence,  iQ,  at  the  base  of  the  crust  are  assumed. 

In  the  eastern  part  of  the  Basin  and  Range  province,  a  major  and 

persistent  discontinuity  at  a  depth  of.  25  to  28  km  has  been  confirmed  in¬ 

dependently  by  several  teams  of  investigators  (Berg  and  others,  1960; 
Pakiser  and  Hill,  1963;  Ryall  and  Stewart,  1963).  The  principal  use  of 
this  diagram  is  to  show  that  any  PS  converted  wave  generated  at  the 
25-km-depth  boundary  can  be  expected  to  arrive  about  3  to  3  1/2  sec 

after  the  first  P  arrival.  It  should  be  noted  that  for  a  P  to  S  con¬ 

version  at  a  depth  of  25  km,  there  is  a  dependency  of  A  t  on  the  angle 
of  incidence  that  theoretically  should  be  observable  on  seismograms. 

For  example,  the  time  interval  At  for  an  angle  of  incidence  of  about 
50°  is  nearly  1/2  sec  greater  than  the  At  for  an  angle  of  incidence 
of  20°. 

Arrivals  of  shear  (SV)  waves  that  are  observed  at  these  time  inter¬ 
vals  after  P  in  this  study  are  identified  as  PS  converted  waves  genera¬ 
ted  at  the  25-km-depth  boundary. 

Earthquakes 

Figure  7-17  shows  the  vertical  particle-motion  diagram  for  an  earth¬ 
quake  (magnitude  not  given  by  U.S.C.&G.S.)  from  Mexico  on  Oct.  19,  1962, 
with  a  depth  of  focus  of  53  km.  The  first  P  arrival  (Fig. 7-17  A)  has  an 
apparent  angle  of  incidence  at  the  surface  of  about  30°  with  the  vertical, 
which  is  in  accord  with  theory.  The  P  motion  is  essentially  linear  and 
prograde. 
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Figure  7-17.  Vertical  particle-motion  diagrams,  Mexico  earthquake,  Oct.  19,  1962,  time  interval  0.00  to 
31.05  sec  after  first  P  arrival;  0.00  *  21h  26m  <.^.3s  GCT. 


In  Figure  7-17B,  the  PS  converted  wave  apparently  generated  at  the 
25-km-deep  layer  begins  manifesting  itself  at  about  3.5  sec  by  changing 
noticeably  the  slope  of  the  segment  of  the  line  on  the  particle-motion 
diagram;  anc  after  "phasing  in"  for  about  1  sec,  becomes  pronounced  at 
about  4.50  sec  and  shows  excellent  S-like  motion  for  half  a  second.  The 
peak-tc-peak  amplitude  of  the  shear  motion,  as  measured  on  the  vertical- 
particle-motion  diagram,  is  about  one- third  that  of  the  P,  similarly 
measured.  It  should  be  noted  that  the  P  wave  starts  its  large-ampli¬ 
tude  motion  about  1  sec  afte*  the  first  arrival  of  P  energy;  conse¬ 
quently,  the  PS-converted-wave  energy  similarly  starts  its  large  mani¬ 
festation  at  4.50  sec. 

In  Figure  7-17C  "he  PS  conv<  ted  -wk  apparently  ge -y.c.  at  the. 
72-km-depth  boundary  begins  manifesting  itself  at  about  7.95  sec  by  chang¬ 
ing  noticeably  the  slope  of  the  segment  of  the  line  on  the  particle- 
motion  diagram;  and  after  "phasing  in"  for  about  ha^r  a  second,  it  becomes 
pronounced  at  about  8.50  sec  and  shows  excellent  S-like  motion  for  about 
1  sec  thereafter.  The  peak-to-peak  amplitude  of  the  shear  motion,  as 
measured  on  the  vertical  particle-motion  diagram,  is  about  one-third 
that  of  the  P,  similarly  measured. 

From  this  same  earthquake,  the  PS  converted  wave  in  Figure  7-17G  is 
apparently  generated  at  the  25-km-depth  boundary  frotr  the  pP  body  wave  in 
Figure  7-17E.  Also  the  PS  converted  waves  in  Figure  7-17H  and  7-171  are 
apparently  generated  at  the  25-  and  72-km-depth  boundaries,  respectively, 
by  the  sP  body  wave  in  Figure  7-l'7G.  Moreover,  the  PS  converted  wdve  in 
Figure  7-17H  may  conceivably  include  some  shear  energy  generated  at  the 
72-km-depth  boundary  by  the  pP  body  wave  in  Figure  7-17E,  as  this  is  the 
time  at  which  this  energy  should  arrive. 
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Figure  7-18  shows  the  horizontal  particle-motion  diagrams  for  the 
same  time  period  for  the  same  Mexico  earthquake.  During  the  early 
arrivals,  the  motion  is  essentially  linear  and  essentially  along  the 
direction  of  a  great  circle  to  the  epicenter;  later,  however,  the  motion 
is  intermittently  elliptical  and  near-linear.  There  is  a  tendency  for 
some  of  the  PS  converted  wave  arrivals  to  initiate  some  elliptical  mo¬ 
tion,  which  suggests  that  the  converted  wave  is  not  a  simple  SV  w^ve 
but  is  complex. 

Figure  7-19  shows  the  vertical  particle-motion  diagrams  for  an 
earthquake  on  May  17,  1963  from  the  Kurile  Islands  with  a  Richter  magni¬ 
tude  of  5.9  and  a  depth  of  focus  of  33  km.  The  PS  converted  wave  shown 
in  Figure  7-19B  is  apparently  generated  at  the  25-km-depth  boundary.  Its 
peak-to-peak  magnitude  is  about  one- third  that  of  P,  as  measured  on  these 
diagrams. 

No  recognizable  PS  converted  wave  generated  at  the  72-kra-depth 
boundary  was  observed  from  this  earthquake. 

Figure  7-20  shows  the  P-  and  S-like  motion  dot  presentation  and  the 
digital-to-analog  seismogram  of  this  Kurile  Islands  earthquake.  The  ex¬ 
cellent  character  of  the  shear  motion  in  the  PS  converted  wave  generated 
at  the  25-km-depth  boundary  and  with  a  time  interval,  At,  of  about  4  sec 
after  P,  is  shown  on  (1)  the  dot  presentation  (A)  and  (2)  the  radial  hori¬ 
zontal  motion  trace  (B)  and  the  vertical  trace  (C) ,  wherein  the  two 
traces  (B  and  C)  are  180°  out  of  phase. 

Figure  7-21  shows  the  vertical  particle-motion  diagram  for  an  earth¬ 
quake  on  Feb.  20,  1964  from  the  Kurile  Islands  with  a  Richter  magnitude 
of  5.2  and  a  depth  of  focus  of  50  km.  The  PS  converted  waves  in  dia¬ 
grams  C  and  F  are  apparently  generated  at  the  25-  and  72-km-depth 
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Figure  7-18.  Horizontal  particle-motion  diagrams,  Mexico  earthquake,  Oct.  19,  1962,  time  interval  0.00  to 
31.05  sec  after  first  P  arrival;  0.00  *  21h  26m  26.3s  GCT. 
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Figure  7-19.  Vertical  particle-motion  diagrams,  Kurile  Islands  earthquake,  May  17,  1963,  time  interval  0.90 
to  14.75  sec  after  firrc  P  arrival;  0.00  =  04h  17l.  26.6s  GCT. 


Figure  7-20.  (A)  P-  an 

Kurile  Islands  earth 
Q4h  17m  26.6s  GCT. 


Figure  7-21.  Vertical  particle-motion  diagrams,  Kurile  Isa.ids  earthquake,  Feb.  20,  1964,  time  interval 
12.00  to  23.95  sec  after  first  P  arrival;  0.00  »  lOh  04m  43.6s  GOT. 


boundaries,  respectively,  by  the  pP  phase  in  diagram  A.  The  other  PS 
converted  wave  in  diagram  F  is  apparently  generated  at  the  25-km-depch 
boundary  by  the  sP  phase  in  diagram  D. 

Figure  7-22  shows  the  P-  and  S-like-motion  dot  presentation  and  the 
digital-to-analog  seismogram  of  this  Kurile  Islands  earthquake.  The  body 
phases  and  the  PS  converted  waves  just  discussed  are  shown. 

Comparison  with  Hannon’s  theoretical  results 

Figure 7-23C  shows  one  of  the  theoretical  crustal  models  (designated 
CAO  used  by  Hannon  (1964)  in  making  his  synthetic  seismograms  computed 
by  the  Haskell-Thomson  matrix  method.  Although  the  layering  of  this 
model  is  not  strictly  applicable  to  the  eastern  part  of  the  Basin  and 
Range  province,  the  values  of  the  velocities  are  sufficiently  close  to 
them  so  that  his  theoretical  results  can  be  compared  qualitatively  with 
our  observed  results.  We  have  added  to  his  model  the  appropriate  ray 
paths  for  the  incident  P  waves  and  the  PS  converted  waves  generated  at 
the  base  of  *.he  crust.  His  theoretical  results  are  for  two  angles  of 
incidence,  namely  53°  and  33°,  which  are  shown  in  diagram  C.  Hannon’s 
computations  included  the  effect  of  a  layered  half-space. 

Figures  7-23A  and  B  show  Hannon's  (1964)  theoretical  synthetic 
seismograms  from  his  theoretical  model  CAO.  The  top  two  traces  (A) 
are  the  radial-horizontal  and  vertical  traces,  respectively,  for  an 
angle  of  incidence  of  53°;  and  the  bottom  two  traces  (B)  are  the 
corresponding  traces  for  an  angle  of  incidence  of  33°. 

Two  ?S  converted  waves  are  recognizable:  (1)  one,  generated  at 
the  base  of  the  25.43-km-depth  boundary,  which  arrives  with  a  time  inter¬ 
val,  ^v,  of  3.5  sec  after  the  first  P 
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Figure  7-23.  Synthetic  seismograms  of  surface  motion  (Hannon,  1964, 
his  CAO  model,)  for  vertical  and  radial-horizontal  components  for 
angles  of  incidence  of  (A)  50°  and  (B)  33°  at  base  of  crust  for  (C) 
earth  crustal  model  of  Gutenberg,  on  which  are  drawn  P  waves  and  PS 
converted  waves  generated  at  base  of  crust. 


arrival;  and  (2)  another  (with  the  arrow  after  Hannon),  generated  at  the 
32.43-km-depth  boundary  (i.e.,  at  the  base  of  the  crust),  which  arrives 
with  a  time  interval,  .At,  of  about  4  1/2  sec  (actually  4.4  sec  after 
the  first  P  arrival.  These  time  intervals  are  for  the  incident  P  waves 
with  an  angle  of  incidence  of  53°  at  the  base  of  the  crust;  the  time  in¬ 
terval  A  t  is  0.4  tp0.4  sec  smaller  in  each  case  for  the  P  waves  with 
an  angle  of  incidence  of  33®  at  the  bas*»  of  the  crust. 

The  theoretical  amplitudes  of  the  PS  converted  waves  generated  at 
the  25.43-km-depth  boundary  and  at  the  32.43-km-depth  boundary  (base  of 
the  crust),  as  measured  on  the  radial-horizontal-component  trace,  are 
about  10  to  20  percent  of  the  amplitudes  of  the  first  P  arrival,  similar¬ 
ly  measured.  The  theoretical  amplitudes  of  the  PS  converted  waves,  as 
measured  on  the  vertical-component  traces,  are  small. 

In  conclusion,  our  observed  results  give  excellent  experimental 
support  of  Hannon's  (1964)  theoretical  results  for  his  model  CAO.  Our 
work  was  initiated  several  years  ago  and  was  done  quite  independently  of 
that  of  Hannon. 

Underground  explosions  in  Semipalatinsk  area,  USSR 

Figure  7-24  shows  the  vertical  particle-motion  diagrams  for  the  first 
10  sec  after  the  firsc  P  arrival  for  the  underground  explosion  in  the  Semi¬ 
palatinsk  area,  USSR,  on  July  19,  1964.  The  emergence  is  very  steep.  It 
is  possible  that  the  PcP  wave,  which  should  arrive  about  1  re c  after  P, 
may  be  "phased  in"  with  the  P  and  is  giving  an  erroneous  impression  of 
the  steepness  of  P.  The  motion  phases  into  S-like  motion  at  about  2.00 
sec  and  continues  thus  until  3.25  sec;  but  this  value  of  2.00  sec  is  some¬ 
what  ear'ier  than  could  be  expected  for  a  PS  converted  wave  generated  at 
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Figure  7-24.  Vertical  particle-motion  diagrams,  probable  underground 
explosion  in  the  Semipalatinsk  area,  USSR,  July  19,  1964,  time  inter¬ 
val  0.00  to  10.00  sec  after  first  P  arrival;  0.00  =  06h  13m  01.4s 
GCT. 


the  25-km-depth  boundary,  S-like  motion  starting  at  7  3/4  sec  and 
lasting  about  1  1/4  sec  is  believ'  to  be  a  PS  converted  wave  generated 
at  the  72-km-deplh  boundary. 

The  angle  of  incidence  of  the  P  wave  at  the  25-km-depth  boundary  is 
apparently  too  small  to  generate  PS  converted  waves  in  this  case. 

Figure  7-25  shows  the  P-  and  S-like  motion  dot-presentation  of  the 
same  Soviet  explosion.  Note  the  excellent  S-like  motion  at  2  to  3  sec 
and  at  7  3/4  to  9  sec. 

Figure  7-14  shows  the  digital  to-analog  seismogram  for  the  same  Soviet 
explosion.  Note  the  excellent  S-like  motion  at  2  to  3  sec  and  at  7  3/4 
to  9  sec. 

Figure  7-26  shows  the  vertical  particle-motion  diagrams  for  the  first 
9.7  sec  after  the  first  P  arrival  for  the  underground  explosion  in  the 
Semipalatinsk  area,  USSR,  on  March  15,  1964.  The  character  of  these 
particle-motion  diagrams  is  almost  identical  to  those  of  the  ones  just 
shown.  The  S-like  moticn  is  present  at  about  2  to  2  1/2  sec,  and  also 
at  7.3  to  8.2  sec.  In  this  case,  however,  about  half  a  second  should  be 
added  to  these  values  because  the  digitizing  was  started  half  a  second 
late. 

Figure  7-27  shows  the  digital- to-analog  seismogram  of  this  same 
Soviet  explosion. 


Sunanary 

We  have  shown  that  PS  converted  waves  have  been  observed  at  various 

probable 

teleseisraic  distances  from  various  earthquakes  and  from  certain  Soviet  / 
underground  explosions.  From  Mexican  earthquakes--that  is,  at  epicentral 
distances  of  about  2,000  km,  PS  converted  waves  generated  at  both  the 
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Figure  7-25.  (A)  P-  and  S- like-motion  dot  presentation  and  (B  to  E)  digital- to-analog  seismogram  of  probable 

underground  explosion  in  the  Semipalatinsk  area,  USSR,  July  19,  1964,  time  interval  0.00  to  20.00  sec; 

0.00  -  06h  13m  01.4s  GCT. 


Figure  7-26.  Vertical  particle-motion  diagrams,  probable  underground 
explosion  in  the  Semipalatinsk  area,  USSR,  March  15,  1904,  time 
interval  0.00  (approximately)  to  9.70  sec  after  first  F  arrival- 
0.00  -  08h  13m  01.7s  GCT. 
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Figure  7-27 .  Digital-to-analog  seismogram  of  probable  underground  explosion  in  the  Semipalatinsk  area, 
USSR,  March  15,  1964;  0.00  *  08h  13m  01.7s  GCT. 


25-km-depth  boundary  and  the  postulated  72-kw-depth  boundary  have  been 
observed.  At  intermediate  epicentral  distances  of  50°,  PS  converted 
waves  generated  from  both  these  boundaries  are  also  observed.  At  an 
epicentral  distance  of  90°  (10,000  km),  PS  converted  waves  generated 
from  the  72-km-depth  boundary  are  observed;  but  the  PS  converted  waves 
from  the  25-km-depth  boundary  v/ere  not  clearly  observed.  The  observed 
amplitudes  of  the  PS  converted  waves  were  generally  about  one-fifth  to 
one-third  of  that  of  the  parent  P  phases  that  generate  the  PS  converted 
waves.  The  PS  converted  waves  from  earthquakes  at  teleseis>  ic  distances 
were  observed  following  the  phases  P,  pP,  sP,  and  PcP.  Thus,  for  earth¬ 
quakes,  the  PS  converted  waves  add  to  the  complexity  of  the  coda  cf  P  in 
addition  to  the  phases  P,  pP,  sP,  and  PcP.  The  PS  converted  waves  from 
the  Soviet  probable  underground  explosions  were  apparently  generated 
from  the  72- km- depth  boundary  only  and  for  the  P  arrival  only;  however, 
in  this  case  also,  the  PS  converted  waves  add  somewhat  to  the  complexity 
of  the  coda  of  P. 

We  consider  that  our  observed  results  have  given  experimental  support 
to  Hannon's  theoretical  seismograms  using  the  Haskell-Thomson  matrix  formu¬ 
lation  for  ground  motion  for  his  model  CAO;  in  particular,  the  PS  conver¬ 
ted  waves  have  been  observed,  as  predicted  by  him  theoretically. 

Last--but  by  no  means  least  --we  consider  that  our  results  give 
added  support  to  the  existence  of  a  velocity  discontinuity  at  a  depth  of 
about  72  km  beneath  the  eastern  part  of  the  Basin  and  Range  province,  as 
first  postulated  by  Berg  and  others  (1960). 
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CHAPTER  VIII 


SUMMARY  AND  CONCLUSIONS 

S-Iike  motion  appears  on  many  of  the  particle-motion  diagrams  at  epproKi 
mately  the  expected  arrival  time  for  PPS  and  SPS.  Usually,  no  other 
S-like  motion  is  apparent  on  the  particle-motion  diagrams  in  the  first 
ten  seconds  except  at  these  times,  approximately  4.0  and  8.0  sec  for  PPS 
and  SPS,  respectively. 

The  most  consistent  arrival  of  large  amplitude  in  the  early  part  of 
the  seismogram  is  an  event  which  we  have  called  Pm^P  in  Chapter  5.  The 
arrival  time  of  this  event  on  the  seismograms  usually  coincides  very  well 
with  the  theoretical  arrival  time  of  a  reflection  from  the  Mohorovicic 
discontinuity.  The  arrival  is  consistent  in  both  time  and  character,  as 
may  be  seen  from  examining  the  various  digital-to-analog  CalComp  plots  in 
Chapter  5.  This  phase  may  be  a  guided  wave  in  the  crust,  usually  called 
P,  and  not  a  simple  reflection  from  the  base  of  the  crust.  However  the 
cravel-times  of  either  Pj^P  or  P  would  probably  be  very  nearly  the  same. 

The  large  amplitudes  following  P^P  are  attributed  to  multiple  reflections 
of  Pft^P  or  P  in  the  crustal  waveguide. 

Some  evidence  for  an  SPP  transmission  path  is  suggested  by  large 
amplitude  P-motion  at  the  expected  time  of  arrival  of  this  event. 

Hannon  (1964)  showed  that  SV  motion  may  be  accentuated  by  the  presence 
of  thin  low- velocity  layering  beneath  the  recording  site.  This  would 
effectively  raise  the  signal-to-noise  ratio  where  13  convertad  waves  are 
concerned  and  may,  for  example,  partly  explain  large  amplitudes  for  PPS. 

The  characteristics  of  the  S-like  motion  defining  PPS  and  SPS,  and 
other  arrivals,  will  need  to  be  studied  further  in  order  to  define  tne 
differences  in  explosion- generated  shear  and  earthquake-generated  shear, 
and  to  use  these  differences  to  aid  in  distinguishing  between  underground 
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nuclear  explosions  and  earthquakes.  In  particular,  the  signal- to-noise 
ratio  on  the  seismograms  of  underground  nuclear  explosions  as  recorded 
at  a  given  station  can  be  increased  by  "stacking"  the  seismograms  of 
several  blasts,  thus  obtaining  an  average  seismic  signature  which  may 

and 

be  characteristic  of  the  local  blast  site.  Also,  by  stacking  seismograms/ 
using  a  reference  time  measured  with  respect  to  Pn,  shear  arrivals  such 
as  PPS  anu  SP3  which  arrive  after  Pn  at  a  time  independent  of  the  epi- 
central  distance,  will  have  their  signal- to-noise  ratio  increased  over 
what  it  was  for  a  single  recording.  This  can  be  done  for  many  of  the 
LRSM  Geotech  stations  for  regions  over  which  the  crusts.!  section  does  not 
change  very  muc...  Such  a  stacking  procedur  would  cause  other  arrivals, 
whose  arrival  time  as  measured  after  P  Is  a  function  of  the  epicentral 
distance,  to  become  less  co! erent  on  the  particle-motion  d  grams.  The 
signal-to-noise  ratio  of  PPS  and  SP3  should  therefore  be  increased  at 
the  expense  of  distance-dependent  arrivals. 

Average  seismic  signatures  obtained  by  stacking  same-site  explosions 
as  recorded  at  the  same  station  can  then  be  cross-correlated  usi.ig  the 
first  20  to  30  seconds,  say,  of  the  entire  average  signature,  with  a 
single  event  recorded  at  a  later  time,  A  cross-correlation  index  will 
then  aid  in  confirming  the  later  event  as  an  underground  nuclear  explosion. 

The  recordings  of  the  LRSM  Geotech  stations  do  confirm  that  a  large 
amount  of  shear  energy,  at  least  equal  to  the  corapressional  energy,  is 
generated  by  some  underground  nuclear  explosions.  At  the  expected  time 
of  arrival  of  the  body  S-wave  (SSS)  large  amplitudes  are  observed  on 
many  of  the  LRSM  recordings.  A  few  underground  nuclear  explosicns  were 
examined  out  at  this  time  interval,  and  the  particle-motion  diagrams 
showed  shear  motion,  beginning  at  the  correct  time  for  the  SSS  trans- 
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mission  path.  The  large  amount  of  original  shear  also  tends  to  support 
the  identification  of  SPP. 

Seismograms  obtained  from  confirmed  underground  nuclear  explo¬ 
sions  of  Richter  magnitudes  3.9  to  5.8  at  the  Nevada  test  site  (NTS)  and 
4  probable  underground  explosions  in  the  Semipalatinsk  area,  USSR  of 
Richter  magnitude  5.5  or  greater  were  recorded  by  the  Benioff  short- 
period  matched  seismometers  (magnification  400,000)  at  Dugway,  Utah, hand- 
digitized,  and  analyzed  (using  particle-motion  diagrams  and  digital- to- 
analog  seismograms)  during  the  first  15  to  25  sec  after  the  first  arrival 
of  F.  Thf  epicentral  distances  were  about  450  and  9,800  km,  respectively. 
Provided  the  explosions  occurred  in  the  same  area,  a  striking  consistency 
of  the  body-wave  phases  was  observed  from  some  (but  not  universally  from 
all)  underground  explosions  within  that  area.  The  characteristics  of  the 
body-wave  patterns  vary  for  three  different  areas  of  detonation  within 
NTS,  presumably  because  of  different  source  functions  (depending  on  rock 
type,  geologic  environment,  degree  of  coupling,  yield,  etc.)  This 
striking  consistency  of  body-wave  phases  has  been  used  successfully  for 
the  positive  detection  and  identification  of  certain  underground  nuclear 
explosions  detonated  at  NTS  and  recorded  at  Dugway,  Utah  and  to  correctly 
determine,  within  a  few  kilometers,  the  area  of  the  detonation  within  NTS. 

The  following  pnases  were  tentatively  identified:  PjjjP  and/or  PPPPPPP , 

a 

SSPPP,  SSPSS,  PnP  (or  P)  and/or  PPP,  PPS,  SPP,  SPS,  and  P’  (or/multiple 
P  reflection);  the  times  of  arrivals  of  these  phases  agree  moderately  well 
with  the  crustal  model  of  the  area  determined  by  the  refraction  method. 
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Seismograms  obtained  from  (a)  10  earthquakes  of  Richter  magnitudes 
5  to  7  in  the  Kurile  and  Tonga  Islands,  Japan,  and  Mexico  and  (h)  the  same 
probable  explosions  in  the  USSR  mentioned  above,  were  recorded  at  Dugway, 
Utah  at  epicentral  distances  of  (a)  2,300  to  9,700  km  and  (b)  about  9,800 
km,  respectively,  and  similarly  digitized  and  analyzed  during  the  first 
20  or  40  sec  after  the  first  P  arrival.  PS  converted  waves  were  generally 
observed  from  both  the  earthquakes  and  underground  explosions.  These  con¬ 
verted  waves  were  apparently  generated  at  both  the  25-  and/or  72-km-depth 
boundaries;  and  the  time  intervals  ^  t,  after  the  first  arrival  of  P, 
were  about  3  to  3  1/2  sec  and  7  1/2  to  8  1/2  sec,  respectively.  These 
time-interval  values  are  in  excellent  accord  with  the  crustal  model  of 
Berg  and  others  (1960).  Similar  PS  converted  waves  with  approximately 
the  same  At  values  were  generally  observed  following  other  large  phases, 
including  pP,  sP,  and  PcP,  from  the  ea  hquakes.  The  time  interval.  At, 
was  generally  somewhat  less  for  earthquakes  with  greater  epicentral 
distances  because  of  the  smaller  angles  of  incidence  of  the  incident 
P  waves  at  the  base  of  the  crust.  The  observed  maximum  amplitudes  of 
the  PS  converted  waves  as  measured  on  che  seismograms  are  usually  one- 
fifth  to  one-third  of  that  of  the  parent  P  phases  that  generate  the  PS 
converted  waves.  These  results  give  experimental  support  to  the  theoreti¬ 
cal  synthetic  seismograms  for  the  CA0  model  of  Hannon  (1964).  Moreover, 
the  observed  PS  converted  waves  give  added  support  to  the  existence  of  a 

A  *  $ 

velocity/^continuity  at  a  depth  of  about  72  km  beneath  the  eastern  part  of 
the  Basin  and  Range  province,  as  first  postulated  by  Berg  and  others  (1960). 
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An  important  consequence  of  these  results  is  that  (1)  for  the 
Soviet  probable  underground  nuclear  explosions  in  the  Semipalatinsk 
area,  the  PS  converted  waves  add  somewhat  to  the  complexity  of  the 
coda  of  P,  and  (2)  for  earthquakes  at  teleseismic  distances,  the  PS 
converted  waves  add  to  the  complexity  of  the  coda  of  P  in  addition  to 
the  phases  P,  pP,  sP,  and  PcP. 
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1074,  October  1963. 
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Acoustics  and  Seismics  Lab.,  Inst,  of  Science  and  Technology,  Univ.  of 
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Cook  and  J.  K.  Costain,  Program  of  papers  submitted  at  annual  meeting 
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PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 


Computer  Programs 
Theoretical  Programs 


DR- 13  (FORTRAN-  IBM  7040) 

Generalized  P  and/or  S  refraction  time-distance  curves 
for  up  to  seven  layers;  also  layer  thickness  determination 
from  P-wave  refraction  data  using  intercept  times. 

X-l  (FORTRAN) 

Effect  of  free  surface  on  incident  plane  SV  wave. 

(Jeffreys  treatment) 

I  (FORTRAN) 

Computes  Zvolinskii  head-wave  coefficients  for  a  PPP  head  wave. 

5  (FORTRAN) 

Computes  reflection  travel  times  of  a  wave  converted  from  P  to  S 
on  reflection  from  a  single  dipping  interface  for  varying  angles 
of  incidence. 

6  (FORTRAN) 

Computes  reflection  travel  times  for  a  P  wave  reflected  from  a 
single  dipping  interface. 

7  (FORTRAN) 

Computes  arrival  time  for  a  PS  converted  wave  formed  at  the 
base  of  the  nth  crustal  layer. 

9  (FORTRAN) 

Zoeppritz  equations  applied  to  a  multilayered  crust.  Computes 
amplitude  ratios  of  PS  converted  waves.  Effects  of  wave  length 
not  included. 

9A  (FORTRAN) 

Zoeppritz  equations  applied  to  a  multilayered  crust  with  free 
surface  effects  included.  Computes  amplitude  ratios  of  PS 
converted  waves.  Effects  of  wave  length  not  included. 

II  (FORTRAN) 

Computes  travel  times  for  refracted  PPP,  PPS,  SPS,  and  SSS 
head  waves  front  a  single  dipping  interface. 

18  (FORTRAN  -  IBM  7040) 

Energy  ratios,  amplitude  ratios,  and  phase  angles  for  an  incident 
P  wave.  Uses  Zoeppritz  equations  and  Knott's  energy  equation 
applied  to  a  single  interface. 


-141- 


PROGRAM 


PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 


19  (FCRTRAL  -  IBM  7040) 

Energy  ratios,  amplitude  ratios,  and  phase  angles  for  an  incident 
SV  wave.  Uses  Zoeppritz  equations  and  Knott's  energy  equation 
applied  to  a  single  interface. 

21  (FORTRAN  -  IBM  7040) 

Computes  reflected  time-distance  curves  for  up  to  twenty  layers 
with  plane  parallel  interfaces. 


Utility  Programs 


UP-1  (AUTOCODER  -  IBM  1401) 

BCD  (HER)  tape  to  cards  or  print  (IBM  1403  printer). 

UP -2  (AUTOCODER  -  IBM  1401) 

Digital  to  analog  plot  of  HER  data  (IBM  1403  printer,  20  lines/ 
inch) . 

UP-4  (AUTOCODER  -  IBM  1401) 

Binary  (STL)  tape  to  cards  or  print  (IBM  1403  printer). 

UP- 5  (AUTOCODER  -  IBM  1401) 

Digital- to- analog  trace  plot  using  card  input  with  two 
sets  of  data  per  card  (IBM  1403  printer  and  plotting  chain, 

20  lines/inch). 

UP- 6  (AUTOCODER  -IBM  1401) 

Binary  (UED)  tape  to  cards  or  print  (IBM  1403  printer). 

UP -7  (AUTOCODER  -  IBM  1401) 

Digital  to  analog  trace  plot  using  card  input.  One  set  of  data 
per  card.  (IBM  1403  printer  and  plotting  chain,  20  lines/inch). 

DR-15D  (FORTRAN  -  MAP  -  IBM  7040) 

Punches  cards  from  binary  (STL)  extract  tapes  written  with  FORTRAN. 

DR-16B  (FORTRAN  -  MAP  -  IBM  7040) 

Lists  binary  (STL)  tapes. 

DR-16C  (FORTRAN  -  MAT  -  IBM  7040) 

Writes  a  binary  extract  tape  using  FORTRAN  from  binary  (STL) 
tapes. 

DR-18TC  (FORTRAN  -  MAP  -  IBM  7040) 

Punches  cards  from  BCD  (HER)  tapes. 

DR- 22  (FORTRAN  -  MAP  -IBM  7040) 

Punches  cards  from  binary  (UED)  tapes. 
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Data  ysis  Programs 


PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

FROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 


DK-2  (FORTRAN) 

Forms  an  L-Z  product  trace  to  accentuate  shear  and  compress ional 
wave  arrivals, 

DR- 7  (FORTRAN) 

Scans  BCD  (HER)  data  for  shear  motion. 

DR-10B  (FORTRAN  -  IBM  7040) 

Composition  of  vertical  motion  from  east-west  and  nortn-south 
directions  to  ^  vertical  plane  containing  the  azimuth  to  the 
epicenter, 

DR-12B  (FORTRAN  -  IBM  7040) 

Autocorrelation  and  power  density  spectrum. 

DR-15B  (FORTRAN  -  MAP  -  ISM  7040) 

Digital-to- analog  plot;  program  lists  data  using  extract  tapes 
written  with  FORTRAN  from  binary  (STL)  tapes;  uses  plotting 
chain  (10  lines/inch)  of  IBM  1403  printer. 

DR-13C  (FORTRAN  -  MAP  -  IBM  7040) 

Digital- to- analog  plot  -  Program  lists  data  from  binary  (STL) 
tapes;  uses  plotting  chain  (10  lines/inch)  IBM  1403  printer. 

DR-18H  (FORTRAN  -  MAP  -  IBM  7040) 

Digital-to -analog  plot  -  Program  lists  data  from  BCD  (HER)  tapes 
uses  plotting  chain  (20  lines/inch)  of  IBM  1403  printer) 

DR-18F  (FORTRAN  -  MAP  -  IBM  7040)  digital  to  analog  plot. 

Program  lists  data  from  BCD  (HER)  tapes  using  a  faster  plot 
subroutine  but  having  less  resolution  (uses  plotting  chain  (10 
lines/inch)  of  1403  printer) 

DR-18P  (FORTRAN  -  IBM  7040) 

Digital-to-analog  plot  using  card  input  (uses  IBM  1403  printer, 
10  lines/inch). 

DR-18FC  (FORTRAN  -  IBM  7040) 

Alphanumeric  particle-motion  plot  and  listing  from  card  input 
(uses  IBM  1403  printer,  20  lines/inch). 

DR-18FP  (FORTRAN  -  MAP  -  IBM  7040) 

Alphanumeric  particle-motion  plot-listing  for  BCD  (HER)  tapes 
(used  IBM  1403  printer,  20  lines/inch). 

DR-18FP1  (FORTRAN  -  MAP  -IBM  7040) 

Alphanumeric  particle-motion  plot  for  BCD  (HER)  tapes  using 
faster  tape  reading  subroutine  SKIPPD  also  gives  data  listing 
(uses  IBM  1403  printer,  20  lines/inch). 
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PROGRAM 

PROGRAM 

PROGRAM 


PROGRAM 

PROGRAM 

PROGRAM 

PROGRAM 


DR-18FP2  (FORTRAN  -  IBM  7040) 

Same  as  DR-18FP1,  but  with  card  input.  Normalizes  the  digitized 
input  data  for  the  alphanumeric  particle  motion  plot.  Uses 
IBM  1403  printer,  20  lines/ inch. 

DR-18RHS  (FORTRAN  -  IBM  7040) 

L-Z  product  trace  using  IBM  7040  U.  of  U.  Computer  Center 
FORTRAN  plot  (uses  IBM  1403  printer,  20  lines/inch). 

DR-24B  (FORTRAN  -  IBM  7040) 

Computes  epicenter- to-station  distance  using  Richter's  method. 
DJG-12  (FORTRAN  -  IBM  7040) 

Prints  out  digital  to  analog  trace  plot  of  east-west,  north- 
south,  vertical  and  radial  components,  the  latter  composed  from 
east-west  and  north-south  seismograms.  Differentiates  the 
seismogram  traces,  and  plots  the  differentiated  traces.  Shows 
by  another  plot  whether  differentiated  motion  is  compression  or 
shear.  The  program  then  resolves  the  horizontal  vector  of 
motion  into  a  component  parallel  to  the  azimuth  and  a  compo¬ 
nent  perpendicular  to  the  azimuth,  and  plots  the  resolved 
components . 

CC-1  (FORTRAN  -  IBM  7040) 

CalComp  particle-motion  plot  for  use  with  Phillips  Fetroleum  Co. 
plot  facilities  at  Idaho  Falls,  Idaho. 

CC-2  (FORTRAN  -  IBM  7040) 

CalComp  particle-motion  plot  for  card  input  (HER)  -  uses  IBM  1627 
plotter  at  U.  of  U.  Computer  Center. 

CC-3  (FORTRAN  -  IBM  7040) 

CalComp  particle-motion  plot  for  card  input  (UED)  -  uses  IBM  1627 
plotter. 

CC-5  (FORTRAN  -  IBM  7040) 

CalComp  digitial  to  analog  trace  plot  of  seismograms  using  card 
input  (HER)  and  IBM  1627  plotter. 
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EXHIBIT  A 


AMPLITUDE,  ENERGY,  AND  PHASE  ANGLES  OF  PLANE 
SV  WAVES  AND  THEIR  APPLICATION  TO  EARTH 
CRUSTAL  STUDIES 


by 

John  K.  Costain,  Kenneth  L.  Cook,  and  S.  T.  Algermissen 


(See  attached  reprint  from  the  Bulletin  of  the  Seismological 
Society  of  America,  v.  53,  no.  5,  p.  1039-1074,  October,  1963) 
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EXHIBIT  B 


TABLES  OF  ENERGY  RATIOS,  AMPLITUDE  RATIOS,  AND 
PHASE  ANGLES  FOR  PLANE  SV  WAVES* 

by 

John  K.  Costain  and  Kenneth  L,  Cook 

ABSTRACT 

This  report  presents  tabulations  of  energy  ratios  from  Knott's 
energy  equation,  amplitude  ratios  i;om  the  Zoeppritz  equations,  and  phase 
angles  computed  from  the  complex  amplitude  ratios  for  a  plane  SV  wave  in¬ 
cident  on  a  plane  elastic  discontinuity.  All  computations  were  programmed 
in  FORTRAN  for  an  IBM  7040  digital  computer.  For  a  Poisson's  ratio  of 
0.250,  incident  angles  were  varied  from  0°  to  88°  in  increments  of  2°, 
except  near  phase  changes  of  180°  in  the  real  domain  and  near  the 
critical  angles  where  the  ratios  were  calculated  in  increments  of  0.25°. 
For  a  Poisson's  ratio  of  0.400,  incident  angles  were  varied  from  0°  to 
85°  in  increments  of  5°.  Both  real  and  imaginary  coefficients  were  con¬ 
sidered  in  the  calculations.  The  varying  parameters  were  the  velocity 
ratio  and  density  ratio  across  each  interface,  and  the  angle  of  inci¬ 
dence.  Compressional  velocity  ratios  (transmitted/incident)  of  0.700, 
0.800,  and  0.900,  and  density  ratios  (transmitted/incident)  of  0.700, 
0.800,  0.900,  and  1.000  were  used.  The  data  are  presented  in  24  tables. 


i 


*  Contribution  No.  57,  Dept,  of  Geophysics,  Univ.  of  Utah,  Salt  Lake  City, 
Utah;  to  be  published  as  VESIAC  Special  Report  No.  4410-96-X;  .'.06  p. 
incl.  illus.,  tables,  refs;  in  press. 
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,  AND  1  MIASM  AND  MAS  Of.  PLANK  SI 
APPLICATION  TO  HAITI  II  (MU'STAL 


S  IM  M )  I  MS 


By  .Ion %  K.  Costain,  Kenneth  L.  Cook,  on  S,  T.  Ai.ueumis.skn 


ABSTRACT 

Knorgy  rut  it  m  from  Hunt  I  ’s  merit  v  e<|Uiil  ioii  anti  amplitude  ratios  and  phase  angles  from  the 
/.neppritz  equations  were  eah’tiiutetl  for  a  plane  .S’l  wave  incident  on  a  plane  elnstii-  diseontin 
uity.  All  computations  were  programmed  in  fortran  for  an  IBM  1(520  digital  eomptiler.  In¬ 
cident  angles  were  varied  from  0’  to  !K>®  in  increments  of  two  degrees  rxeepl  near  the  critical 
angles,  where  the  ratios  were  calculated  in  increments  of  (1.25  degree  generally.  Both  real  anti 
imaginary  coefficients  were  considered  in  the  calculations.  The  varying  parameters  were  the 
velocity  ratio  and  density  ratio  across  each  interface  and  the  angle  <>f  incidence.  Poisson's 
ratio  was  kept  constant  at  11.25.  Compressional  velocity  ratios  (1  r.tnsmined, /incident  i  of  0.7, 
0.8,  and  0.0,  and  density  ratios  ftransniilled/incident)  of  0  7,0.8,  0.0,  and  1.0  were  used.  The 
data  are  presented  in  an  integrated  album  which  consists  of  a  total  of  144  curves  (48  curves  for 
each  of  the  3  characteristics  studied).  The  theoretical  amplitudes  of  PS  converted  waves,  which 
are  seismic  Imdy  waves  resulting  from  flic  conversion  of  an  incident  parent  P  wave  at  a  bound¬ 
ary  within  (lie  earth's  crust  to  a  refracted  vertically  polarized  S t'  wave,  were  computed  (using 
a  Fortran  program  on  an  IBM  1020  digital  computer)  for  several  multi -layered  hypothetical 
and  actual  crustal  models.  Preliminary  results  indicate  that  the  amplitudes  of  successive  PS 
converted  waves  arriving  from  successively  deeper  layers  will  continually  increase  provided 
that  I  he  velocity  ratio*  (upper  iiiedium/lower  medium)  continually  <krrr.a*t  with  depth.  The 
effects  of  the  wave  period  have  not  yet  tieen  considered.  Particle  motion  diagrams  obtained 
from  seismograms  of  two  underground  nuclear  explosions  show  some  shear  motion  in  the  first 
ten  seconds  which  is  provisionally  interpreted  to  he  front  PS  converted  waves. 

Inthomvtiox 

Tim  general  properties  of  plane  AM'  (vertically  polarized  shear)  seismic  hotly 
waves  as  they  impinge  upon  a  plane  elastie  discontinuity  and  are  subsequently 
reflected  and  transmitted,  have  long  been  known.  Recent  studies  indicate  the 
desirability,  however,  of  having  available  a  more  complete  album  of  curves  and 
tables,  than  published  heretofore,  of  specific  characteristics  of  NT  wax  °s  to  assist 
in  the  solution  of  manv  current  problems  in  seismology.  The  purpose  of  this  paper 
is  to  present,  at*  integrated  album  of  cnrx'es  for  the  energies,  amplitudes,  and  phase 
angles  of  waves  derived  from  plane  NT  waxes  which  are  incident  upon  a  plane 
elastie  discontinuity,  and  to  <bo\v  how  such  information  can  lx*  applied  to  studies 
of  earth  crustal  structure,  as  an  example.  It  is  Itelieved  that  the  album  of  eurx’es 
will  have  widespread  use  in  other  types  of  seistrie  problems  in  which  a  quantitative 
exahiation  of  the  characteristics  is  desirable. 

Pkkvkm's  VVokk 

Theoretical  studies  of  the  partition  of  energy  incident  at  an  elastic  interface  and 
the  amplitudes  of  the  resulting  refracted  and  reflected  waves  relative  to  the  incident 
xvax’e  have  Ijeen  made  by  many  inxestigators.  These  include  Knott  (ISAM),  Zoeppritz 
(19151),  Jeffreys  (192b),  Muskat  (1933),  Shriller  (1933),  Slichter  and  Oahricl  (1933), 
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Mnskat  and  Meres  (1010),  t  hiTenlierg  (lOllu),  Medan  (107.%,  l>),  sale  (i077). 
Schwind  (1000),  Richards  (1000,  1001),  Sleinhart  and  others  ( 1001),  Haskell  (100:2), 
and  McCamy  and  others  (1002). 

The  published  data  pertaining  to  incident  .ST  waves  oniy  are  summarized  in 
table  1.  Mtlmugh  the  table  does  not  given  an  exhaustive  list  of  data  and  curves 
that  have  been  published,  it  is  believed  to  summarize  most  of  the  recent  data  and 
curves  available  for  incident  .S’ I'  waves. 

(IkNEHAI,  MkTIIOJ;  or"  Com  CITATION  AM)  At’ITKACY  OF  1)\T\ 

The  determination  of  the  relative  energies  and  amplitudes  of  waves  derived  from 
a  plane  .SI'  wave  incident  on  a  plane  interface  involves  the  solution  of  four  simul¬ 
taneous  linear  equations  (Knott,  1800;  Zocppritz,  1010)  with  real  or  imaginary 
itn  UteH-llls.  The  eqthilitrhs  w.ti  solved  in  the  tial  and  cirtuphw  doinaius  by  the 
(’rout  method  (Crotit,  1041)  which  was  programmed  in  foktkan  i  and  foktkan  ii 
for  use  on  an  IBM  1(520  digital  computer.  The  variable  field  length  capability  of  the 
IliLt)  FoTii  iiA.\  o  SAslt  in  pro  ided  a  means  oi  checking  ihc  accuracy  of  the  solutions, 
most,  of  which  weix:  obtained  using  a  float ing-point  precish  n  of  eight  decimal 
places.  Those  solutions  which  were  obtained  in  foktkan  ii  using  a  floating  point 
precision  of  sixteen  decimal  places  differed  from  those  obtained  using  a  precision  of 
eight  decimal  places  only  in  the  seventh  decimal  place.  A  conservation  of  energy 
cheek  on  all  solutions  was  satisfied  in  every  set  of  computations  to  at  least  the  sixth 
decimal  place  in  foktkan  i,  and  to  the  fourteenth  decimal  place  when  the  higher 
precision  was  used  in  foktkan  ii. 

PaRAME'i  .-.as  1  NVEsTlUATElJ 
Velocity  and  Density  Ratios 

The  angle  of  incidence  of  plane  .ST  waves  on  a  plane  interface  separating  two 
semi -infinite  media  of  different  velocity  and  density  was  varied  in  increments  of 
two  degrees  from  zero  to  ninety  degrees.  In  the  neighborhood  of  the  critical  angles, 
the  increment  ranged  from  0.01  to  0.27  degree  hut  generally  was  held  at  0.27  degree. 
The  values  of  the  density  ratios  (density,  pi,  of  medium  with  transmitted  waves 
density,  p2,  of  medium  with  incident  wave )  wen*  taken  as  0.700,  0.800.  0,000  and 
1.000.  For  each  of  these  density  contrasts,  the  values  of  the  velocity  ratios  (longi 
tudinai  velocity,  >',  of  medium  with  transmitted  waves  longitudinal  velocity, 
of  medium  with  incident  wave)  were  taken  as  0.700,  0.800  and  O.OtX).  In  all  the 
calculations,  Poisson’s  ratio,  a,  was  assigned  the  value  0.2.70  for  media  on  each 
side  of  the  interface.  The  subscripts  2  and  1  refer  to  properties  of  the  medium  with 
incident  waves  and  medium  with  transmitted  waves,  respectively  (see,  for  example, 
figure  2). 

In  addition  to  the  systematic  variations  of  velocity  and  density,  various  hypo¬ 
thetical  and  actual  crustal  sections  were  examined.  The  density  and  velocity  ratios 
involved  in  these  sections  are  shown  on  figures  10  and  20. 

Critical  Angles 

Figure  1  shows  the  reflected  and  transmitted  waves  derived  from  a  plane  ,SF 
wav«  incident  on  a  plane  interface. 


According  to  Snell's  law 


CLANK  A'l  WAV  Kb 
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sin  6  sin  b  sin  r  sind  bin  r  sin  / 

v  ~  ~v  r~  =  r"  ”  v'  ~  z 

l  or  Poisson’s  ratio  equal  to  0.2.40,  V  V  =  Y  Z  —  \/o,  and 

r  =  sin  '(  v  0  sin  8  ? 

l  or  the  sine  of  the  angle  c  (angle  of  reflected  P  wave!  to  lit*  less  Mian  or  equal  to 
one,  8  should  not  exceed  sin  1  (  1  \/8)  =  0*1.204  degrees.  This  is  the  Inst  critical 
angle.  Fora  velocity  ratio  of  Y  I 

e  =  sin  '[(  T  1'}  sin  d\  =  sin'K.r  (')•(('  rising  —  bill  '{( v^*i)  -  <  Y  O  sin  0] 


v  =  longitudinal  velocity 

IN  MEDIUM  I 
7  -  SHEAR  VELOCITY 
IN  MEDIUM  I 

U  *  LONGITUDINAL  VELOCITY 
IN  MEDIUM  2 
V  •  SHEAR  VELOCITY 
IN  MEDIUM  2 


Fie.  I.  Iteflcclcd  and  transmitted  waves  derived  from  a  plane  SF  wave  incident  on  a  plane 

interface  at  an  angle  6. 


For  the  sine  of  the  angle  e  (angle  of  transmitted  P  wave)  to  lie  less  than  or  equal 
to  one,  the  angle  9  should  not  exceed  sin  '|1  [ (  \/8 )  •  (  Y  (  )]'.  For  velocity  ratios 
Y  ('  equal  to  0.700,  0.800  and  0.900,  this  corresponds  to  values  of  the  incident 
angle.  9,  of  .Vi.oOT,  40.194  and  .‘>9.904  degrees,  respectively.  These  are  the  second 
critical  angles.  Physically,  this  means  the  energy'  of  the  reflected  P  wave  becomes 
zero  when  the  angle  of  ineidence  of  the  .VI'  wave  reaches  Ho. 204  degrees,  and  the 
energy  of  the  transmitted  P  wave  becomes  zero  when  the  angle  of  ineidence  of  the 
.ST  wave  reaches  .*>*>..*>1)7,  40,194  and  89.904  degrees  for  velocity  ratios  Y  I  of 
0.700,  0.800  and  0.900,  respectively. 

Amplitude  Patios 

The  Zoeppritz  equations  for  a  plane  .ST  wave  of  amplitude  R  incident  on  a 
plane  interface  are  (Richter,  19.*>8,  pp.  070  078): 


( R  -p  I) )  sin  h  +  ('  cos  a  —  /•.  cos  r  —  F  sin  /  =  0 
i  R  —  i) )  cos  h  T  ('  sin  a  +  /•,’  sin  c  —  h '  cos  /  =  0 


i 


f 


W.WKs 
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'■  li  +  D)  cos  2b  -  C 

■  ill  —  I))  sin  2b  -j- 


sin  2c  ~  FK 


(c)  ""  +  KK  (it)  ’ 

('  (,,,s  +  A’A'  cos  2/  4-  FK 


cos  2/  =  0 


sin  2 !/  =  0 


where  the  anpi’es  are  as  detinetl  in  figure  !.  The  angle  a  is  the  angle  of  incidence  of 
a  compressional  wave  of  velocity  (  which  would  generate  hair  waves  with  the  same 
reference  angles  as  shewn  in  figure  1.  K  is  the  density  contrast,  equal  to  pt  p> , 
where  Pi  and  p->  are  the  densities  of  the  upper  and  lower  media,  respectively. 

The  equations  were  solved  for  C,  l),  E  and  F,  the  amplitudes  of  the  leflected  F 
reflect'  N,  transmitted  F  and  transmitted  .S  waves,  irsjiectively.  A  similar  set  of 
eqnatioi  with  C,  /),  E  and  F  as  the  unknown  amplitudes  exists  for  a  plane  F  wave 
of  amplitude  .1,  incident  on  a  plane  interface.  The  equations  are  not  given  hen*  hut 
are  given  by  Hichter  ( lb. 18,  pp.  U70  U7d).  For  the  analysis  of  PS  converted  waves, 
it  was  necessary  to  consider  the  conversion  of  a  F  wave  incident  on  a  plane  interface. 
A  FonruAN  program  with  the  basic  structure  of  that  used  for  the  incident  «ST  wave 
was  employed.  When  only  a  single  interface  was  involved,  tee  amplitude  of  the 
incident  waves,  .4  :>r  H.  was  set  oq  al  to  unity. 

The  (juantities  f\  I),  E  and  F  am  complex  when  the  angle  of  incidence  of  an  AT 
wave  is  greater  than  2.V2G4  degrees  for  the  parameters  used  in  this  study.  Both  the 
real  and  imaginary  parts  of  the  complex  amplitudes  were  printed  out  to  facilitate 
the  determination  of  the  phase  shifts  associated  with  the  complex  domain.  The 
quantity  having  physical  significance  is  the  modulus  of  the  complex  amplitude, 
which  is  the  square  root  of  the  sum  of  the  squares  of  the  re.  1  and  imaginary  parts 
of  the  complex  amplitude. 

Figures  2b  through  17b  show  the  amplitude  ratios  as  computed  from  the  Zoeppritz 
v  .".ations  for  an  AT  wave  incident  on  a  single  interface.  The  dashed  portions  of  the 
curves  refer  to  angles  of  incidence  for  which  some  of  the  coefficients  of  that  particular 
reflected  or  refracted  wave  amplitude  are  imaginary,  that  is,  the  angle  of  reflection 
or  refraction  is  imaginary.  The  wave  whose  angle  of  reflection  or  refraction  is 
imaginary  does  not  exist,  and  no  physical  significance  should  1*  attached  to  the 
modulii  of  the  i  nplitudes  (or  energies)  within  these  regions. 

Inspection  of  the  amplitude  ratio  curves  on  figures  2b  through  17b  shows  the 
amplitudes'  :o  tie  relatively  insensitive  to  changes  in  the  density  ratio.  Changes  in 
the  velocity  ratio,  however,  cause  marked  changes  in  the  amplitudes.  For  the  most 
part,  the  greater  the  velocity  contrast  (lower  velocity  ratio  of  Y/l"),  the  greater  will 
fie  the  amplitude  of  th°  reflected  or  refracted  dor:  ed  wave.  At  and  near  the  critical 
angles  the  amplitudes  change  rapidly,  sometimes  passing  through  zero  and  resulting 
in  phase  changes  of  180  degree's.  Curve  .‘3  on  figure  10b,  for  example,  shows  a  phase 
change  of  ISO  degrees  near  the  first  critical  .jigle.  The  critical  angles  are  easily 
identified  with  the  sharp  peaks  and  troughs  on  the  craves. 

Small  changes  in  Poisson’s  ratio  hare  little  effect  on  the  amplitude  ratios.  For 
Poisson’s  ratio  equal  *o  0.400  tor  both  media,  the  greatest  changes  were  noticed  for 
transmitted  F  waves,  whose  amplitudes  were  decreased.  More  variation  occurs  if 
the  upper  (transmitted  wave)  medium  has  Poisson  s  ratio  equal  to  0.4(H),  and  the 
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lower  medium  is  assigned  a  value  of  0.2.10.  However,  ' .  -ause  of  the  relative  in¬ 
sensitivity  of  the  amplitude  ratios  to  changes  in  Poisson’s  a*io  foi  conversions  of 
interest  in  this  study  (incident  /*  to  transmitted  P  and  transmitted  S;  and  incident 
S  to  transmitted  P  and  transmitted  «S),  Poisson’s  ratio  was  assigned  the  value  0.2.K) 
for  all  calculations  presented  heroin. 

Energy  Ratios 

The  distribution  of  energy'  among  the  four  waves  derived  from  an  incident  SP 
wave  can  l>e  obtained  from  Knott’s  energy  equation  for  an  incident  .S’  V  wave,  which 
is  (Macciwane,  1932,  p.  173,  Equation  7.30): 

_  (cot  c)C'-  D  ’  pi  (cot  e ) E  "  pi  (cot  f)F>m  . , 

1  ~  (cut  d'.B*  +  It-  +  pt  (cot  d)B'-  +  ps  (cot  </)S'2 

where  the  angles  c,  d,  etc.  are  defined  as  in  figure  1,  and  pi  and  p>  are  the  densities 

of  the  upper  and  lower  media,  respectively.  The  ((uantities  B  ,(' ,  D  ,  E  and  F  arc 

the  amplitudes  of  the  displacement  potentials  of  the  incident  »SP  wave  and  the 
reflected  P,  reflected  .S’,  transmitted  /'and  transmitted  .S  waves,  respectively.  These 
quantities  are  replaced  by  BY,  (T,  DY,  EY  and  FZt  respectively,  where  /i,  (\  I), 
E  and  E  are  the  amplitudes  of  the  particle  displacements  as  defined  in  the  Zoepprit/, 
equations.  Equation  1  then  becomes 

(cot  e)C3{T?  t/V"  pi  (cot  e)E~Y"  ,  pi  (cot  f)F‘Z‘ 

1  ~  (cot ~d)fr-Y-  +  Fp2  +  (h  (cot  d)B-V-  +  p,  (cot  f i)&V* 

or 

_  (cot  c)CJC-  2  Pi  (cot  e)E‘Yi  Pi  lent  f)F’Z ’  ^ 

(cot  d)  P2  PJ  (cot  d)P2  Pi  (cot  d)  \  2 

The  terms  on  the  right  side  of  etjuation  2  represent  the  energies  of  the  reflected  P, 
reflected  „S,  transmitted  P  and  transmitted  S  waves,  respectively,  with  respect  to 
the  energy  of  an  incident  .ST  wave  of  amplitude  B  (B  set  equal  to  unity).  The 
energy  ratios  are  obtained  by  determining  the  amplitudes  C,  D,  E  and  E  from  the 
Zoeppritz  eonations,  inserting  these  values  in  equation  2,  and  then  evaluating  each 
term  in  the  equation.  These  individual  terms  were  then  printed  out  as  the  ratios 
of  the  energy  of  the  reflected/',  reflected  .S,  transmitted  P  and  transmitted S  waves, 
respectively,  to  the  energy  of  an  incident  SPr  wave. 

Equation  2  serves  as  a  useful  conservation  of  energy'  check  on  the  accuracy  of  the 
amplitudes  determined  from  the  Zoeppritz  equations,  and  energy'  ratios  as  computed 
from  this  equation.  In  all  eases,  the  sum  of  the  terms  on  the  right  liand  side  of 
equation  2  differed  from  unity  only  in  the  seventh  decimal  place  at  worst,  when  a 
floating  point  precision  of  e’Tht  decimal  places  was  used  in  the  computations. 

A  similar  conservation  of  energy'  equation  exists  for  an  incident  plane  P  wave  on 
a  nlane  interface  (Maechvane,  1932,  p.  167,  Equation  7.62).  Where  conversion  of 
incident  P  to  transmitted  P  or  S  was  desired,  as  in  the  multilayer  case,  the  con¬ 
servation  of  energy  etjuation  for  a  P  wave  was  similarly  applied  to  cheek  the 
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accuracy  of  the  amplitudes  a*  determined  from  the  Zoeppritz  equations  for  an 
incident  P  wave. 

The  same  generalizations  I  old  for  the  energy  ratio  curves  shown  on  figures  2a 
through  17a  as  for  the  amplitude  ratios,  (lenerally,  increasing  the  velocity  contrast 


Fio.  4.  Plane  SP  wave  incident  on  a  plane  interface  f<*r  a  density  ratio  ftransmitted/in- 
cident)  of  0.900  and  velocity  ratios  ft  ranainii  ted/incident)  of  0,7tX),  0.800  and  0.000.  (a)  Ratio 
of  energy  of  reflected  P  wave  to  energy  of  incident  51  wave.  lb)  Ratio  of  amplitude  of  reflected 
P  wave  to  amplitude  of  incident  51'  wave,  fc.;  Phase  angle  of  reflected  P  wave. 
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will  increase  the  energy  of  the  derived  wave  with  respect  to  the  incident  wave. 
1'he  Hashed  portions  of  the  curves  correspond  to  those  incident  angles  for  which 
certain  of  the  angles  of  reflection  or  refraction  of  a  particular  derived  wave  are 
imaginary.  No  physical  significance  is  to  be  attached  to  these  regions.  Xo  energy 
is  carried  away  from  the  interface  by  these  waves  at  these  angles  of  incidence. 


20  25  K  U  «0  O  »C  ft  M  H  »  K  N  W 
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Pm.  5,  Plane  AT  wave  incident  on  a  plane  interface  for  a  density  ratio  (transmit  ted/in- 
eident)  of  1.000  and  velocity  ratios  {transmitted/incident)  of  0.700,  U.S00  and  0.900.  (a)  Ratio 
of  energy  of  reflected  P  wave  to  energy  of  incident  SP  wave,  (b)  Ratio  of  amplitude  of  reflected 
/'  wave  to  amplitude  of  incident  AT  wave,  (c)  Phase  angle  of  reflected  /'  wave, 
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Fi«,  6.  Plane  ST  wave  incident  on  a  plane  interface  for  a  density  ratio  (transmitted  in¬ 
cident)  of  0.700  and  velocity  ratios  (transmitted/incident)  of  0,700.  0.800  and  0.900,  fa)  Ratio 
of  energy  of  reflected  .St'  wave  to  energy  of  incident  M  wave.  ib)  Ratio  of  amplitude  of  re¬ 
flected  St*  wave  to  amplitude  of  incident  .St  wave,  (c)  Phase  angle  of  reflected  .St  wa\c. 


. . . . . 


Flu.  7.  Plane  SF  wave  incident  on  a  plane  interface  for  a  density  ratio  (transmitted/in¬ 
cident)  of  0.800  and  velocity  ratios  ttransniitted/incident)  of  0.700,  O.SX)  and  0.90C.  fa)  Ratio 
or  energy  of  reflected  »S V  wave  to  energy'  of  incident  SI'  wave,  (b)  Ratio  of  amplitude  of  re 
fleeted  5F  wave  vo  amplitude  of  incident  SV  wave,  (c)  Plmie  angle  of  reflected  51'  wave. 


Flo.  S.  Plane  Slr  wave  incident  on  a  plane  interface  for  a  density  ratio  {transmit ted/in  «  j 

cident)  of  0.000  and  velocity  ratios  (transmitted/incident)  of  0.700.0,800  ami  0.900.  (a)  Ratio 
oi  energy  if  reflected  .ST  wave  to  energy  if  incident  SI'  wave,  (b)  Ratio  of  amplitude  of  re¬ 
flected  .S' I "  wave  to  amplitude  of  incident  SF  wave,  (c)  Phase  angle  of  reflected  SV  wave. 
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Flo.  it.  Plane  ST  wave  incident  on  a  plane  interface  for  a  density  raiio  (transmitted/in 
eident)  of  1.000  and  velocity  ratios  (transmitted/incident  of  0.700,  0,800  and  0.900.  (a)  Hatio 
of  energy  of  reflected  SV  wave  to  energy  oNnciJent  SF  wave,  (b)  Ratio  of  amplitude  of  re¬ 
flected  »ST  wave  to  amplitude  of  incident  .ST  wave,  fc)  Phase  angle  of  reflected  SI’  wave. 
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At  the  critical  angle  lieyoud  which  the  angle  of  reflection  or  refraction  of  a  particular 
derived  wave  becomes  imaginary,  the  energy  of  this  wave  becomes  zero.  In  applying 
the  conservation  of  energy  equation  after  35.204  degrees,  for  the  parameters  used 
in  this  study,  the  terms  in  energy  equation  2  are  complex.  Fhe  imaginary  parts  of 
all  of  these  terms  must  add  up  to  z->ro  and  the  real  parts  to  unity,  in  all  eases.  The 
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stun  of  the  modulii  of  the  complex  energies  of  those  derived  waves  whose  angles  < 
reflection  or  refraction  are  not  imaginary  will  l>e  unity.  The  moduli i  of  those  waves 
whose  angles  of  reflection  or  refraction  are  imaginary  are  not  included  in  the 
summation. 
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Fig.  11.  Plane  SI  wave  incident  on  a  plane  interface  for  a  density  ratio  (transmit led/in¬ 
cident)  of  0.800  and  velocity  ratios  (t ransmi tted/ineideut .  of  0.700,  0  8*1  and  0.000  fa)  Ratio 
of  energy  of  transmitted  /'  wave  to  energy  of  incident  .ST  wave,  (b)  Ratio  of  amplitude  of 
transmitted  /'  wave  to  amplitude  of  incident  SI'  wave,  (c)  Phase  angle  of  transmitted  /'  wave. 
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Phase  Axuufs 

Hev<  Hu  I  tin-  iirst  critical  angle,  the  amplitudes  ef  the  derived  waves  are  complex, 
and  all  wave  amplitudes  will  generally  undergo  phase  ehanges.  The  complex  am 
plitudes  computed  from  the  Zoeppritz  equat  ions  were  printed  out  in  the  form  u  +  il>. 
The  phase  angle,  Q,  of  the  derived  wave  is  given  by 


Flo  12-  Fianc  SF  wave  incident  on  a  plane  interface  fur  a  density  ratio  (transmit led/ill 
cidenl)  of  0.900  and  velocity  ratios  (transmitted  ncidenl )  of  tt.TtX),  tt.StX)  and  tt.iXXt.  (a)  Ratio 
of  energy  of  transmitted  P  wave  to  eti  -rgy  of  incident  SF  wave,  (b)  Ratio  of  amplitude  of 
1  ransrr:- ted  P  wave  to  amplitude  of  incident  SY  wave,  (e.i  Fliase  angl.-  of  transmitted  P  wave. 
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i  b 

4>  =  tan 

a 

The  quantities  a  and  h  identify  the  quadrant  and  thus  the  phase  angle  of  the  derived 
wave. 
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Fits.  13.  Plane  SY  wave  incident  on  a  plane  interface  for  a  density  ratio  (transmitted/ 
incident)  of  1,000  and  velocity  ratios  (transmitted  'incident )  of  0,700,0,^00  and 0.000.  ia>  Ratio 
of  energy  of  transmitted  F  wave  to  energy  of  incid-nt  S1T  wave.  <b)  Ratio  of  amplitude  of 
transmitted  F  wave  to  amplitude  of  incident  »S|'  wave,  (c)  Phase  angle  of  transmitted  P  wave. 
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Fig.  14.  Plane  SF  wave  incident  on  a  plane  interface  for  a  density  ratio  i transmitted/ 
incident)  of  0.700  and  velocity  ratios  (transmitted,  incident )  of  0.700,  O.S00  and  0.900.  fa) 
Ratio  of  encigv  ot  transmitted  SF  wave  to  energy  of  incident  SF  wave,  fh)  Ratio  of  ampli¬ 
tude  ol  transmitted  SF  wave  to  amplitude  of  incident  SF  wave,  (c)  Phase  angle  of  trans¬ 
mitted  SI  wave. 
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Fio.  15.  Plano  SI’  wave  incident  on  a  plane  interface  for  a  density  ratio  (transmitted/ 
incident)  of  0.800  and  velocity  ratios  (transmitted/inrident )  of  0.700,  O.StX)  ami  0.900.  (a) 
Ratio  of  energy  of  transmitted  SI’  wave  l»  energy  of  incident  SV’  wave.  ( 1 » )  Ra'ioof  amplitude 
of  transmitted  <SF  wave  1 1>  amplitude  of  incident  SP  wave.  <c)  Phase  angle  of  transmitted 
SI’  wave. 
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ANGLE  Of  IHClOfSCE,  S  (DC  Of  EES) 

F ,t;  jo.  plane  SP  wave  incident  on  a  plane  interface  for  a  density  ratio  (transmitted/ 
incident)  of  0.000  and  velocity  ratios  (transmitted /incident)  of  0.700  .  0.800  and  O'.MK).  (a) 
Ratio  of  energy’  of  transmitted  ST  wave  to  energy  of  incident  .VI  wave.  (I>)  Ratio  of  an'plitude 
of  transmitted  .ST  wave  to  amplitude  of  incident  .SI'  wave.  >e)  Phase  angle  **f  transmitted 
•SI’  wave. 
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TIk  significance  of  u  phase  change  in  a  derived  wave  is  to  change  the  arrival  time 
off  peak  or  trough  of  the  wave  by  an  amount  SI  =  (<j>7’i  (2ir),  where  T  is  the  pe¬ 
riod  of  the  wave.  If  the  phase  angle,  is  negative,  the  peak  or  trough  is  delayed 
an  amount  St\  if  t lie*  phase  angle  is  positive,  it  is  advanced  in  time  this  amount, 
i  or  example,  a  transmitted  .s'l  wave  of  period  one  second  with  a  phase  angle  of  MX) 
degrees  will  arrive  approximately  0.27  second  early  because  of  the  phase  shift.  Prac¬ 
tical  application  of  phase  angles  would  include  the  correction  of  observed  travel 


Fm.  18.  Polar  coordinate  diagram  of  amplitudes  and  phase  angles  of  a  transmitted  /'  wave 
derived  from  an  incident  5F  wave  for  a  velocity  ratio  (transmitted/incident)  of  0. IKK)  and  a 
density  i  at  to  (Iransmittod/incident )  of  0.700.  Numbers  along  the  curve  refer  to  the  angle  of 
incidence,  0,  of  the  incident  SF  wave  at  a  plane  interface.  The  curve  is  dashed  after  de 
groea,  at  which  point  the  transmitted  H  wave  vanishes. 

times  by  the  appropriate  time  lead  or  lag  determined  from  the  phase  angle,  in 
order  to  locate  more  precisely  the  depths  to  the  interfaces  which  generate  the  de¬ 
rived  waves  (see  e.g.,  Richards.  1901). 

Figures  2c  through  17c  give  the  phase  angles  associated  with  the  complex  ampli¬ 
tudes  computed  front  the  Zoeppritz  equations.  Phase  changes  of  180  degrees  in  the 
real  domain  are  also  noted  on  the  curves  where  they  occur.  It  will  he  noted  that  the 
phase  changes  near  the  critical  angles  are  often  large  and  rapid  (figure  10c),  The 
dashed  portions  of  the  curves  have  no  physical  significance;  they  are  associated 
with  incident  angles  ol  derived  waves  whose  angles  of  reflection  or  refraction  are 
imaginary.  At  grazing  incidence  (9  =  i0  degrees),  for  all  velocity  ratios,  the  phase 
angle  of  a  reflected  .S'l  wave  reaches  180  degrees,  as  shown  on  figures  0c  through  0c. 


AMPLITUDE  OF  TRANSMITTED  P  (P.I/AMPLlTUDE  OF  INCIDENT  P  <P*) 
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Figure  IS  i>  uii  alternate  method  «t  displaying  the  data  shown,  for  example,  on 
figure  10h  and  U»c,  curve  In  particular,  the  polar  coordinate  diagram  shows  how 
the  amplitude  of  a  wave  must  pass  through  zero  \\  hen  a  phase  change  of  ISO  degrees 
takes  place.  This  is,  of  couise,  true  in  either  the  complex  domain  or  the  real  domain. 
In  figure  IS,  the1  real  part  of  the  complex  amplitude  is  plotted  along  the  abscissa, 
and  the  imaginary  part  along  the  ordinate.  The  phase  angle  is  then  read  directly 
from  tin  graph.  Number*  plwtltd  along  th<  curve  irftT thr  rmgfr  of  mridri'ifl ,  ft 
In  the  real  domain,  the  amplitudes  pass  through  zero  with  a  phase  change  of  180 
degrees,  as  also  shown,  for  .ample,  on  figure  Oh  and  IF*,  curves  1 , 2  and  d. 


Flo.  19.  Hypothetical  crustal  section  in  which  the  velocity  ratio  (f ntBsmiUed/ineideni ) 
decreases  with  depth.  The  amplitudes  of  successive  PS  converted  waves  successively  increase 
with  increasing  depth  to  the  interface.  Free  surface  effects  are  not  included.  Arrival  of  PSs 
converted  wave  is  shown  in  insert . 

Application's  of  Data 

The  thickness  of  a  single  crustal  layer  ean  theoretically  be  determined  by  observ¬ 
ing,  with  suitable  instrumentation  at  a  single  station,  the  difference  in  time  of 
arrival  of  the  first  P  wave  and  the  arrival  of  an  .S  wave  (a  PS  converted  wave) 
derived  from  this  P  wave  at  the  base  of  the  layer  (  Vndreev,  1957 ;  Cook  and  others, 
1902).  For  a  multilayered  crust,  PS  converted  wave  arrivals  corresponding  to  each 
interface  will  theoretically  exist.  If  the  compressional  and  shear  velocities  in  each 
of  the  different  layers  are  known,  the  thickness  of  each  layer  ean  be  computed  from 
the  differences  in  time  of  arrivals  of  the  parent  P  wave  and  of  the  individual  PS 
converted  waves. 

Using  seismograms  from  both  earthquakes  and  blasts,  Russian  workers  (Andreev, 
1957,  Andreev  and  Shebalin,  197)7,  Kuz’mina,  I9.’>9)  have  observed  that  the  ampli¬ 
tudes  of  the  various  successive  PS  converted  wave  arrivals  increase  witn  succes¬ 
sively  increasing  depth  to  the  interface  where  the  conversion  took  place.  Bulin 


AMPLITUDE  OF  PS  CONVERT EO  MRVE  (PS.) /AMPLITUDE  OF  INCIDENT  P  (P*  | 
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(lUtiOi,  a-  well  as  Schwiud  and  utliers  (11M50,  y-  11823),  find  that  this  is  not  always 
true*. 

\  preliminary  investigation  of  the  theoretical  amplitudes  to  Ik*  expected  from  > 
converted  waves  shows  that  if  the  velocity  ration  (.upper  medium  lower  medium) 
continually  decrease  with  depth,  then  the  amplitudes  of  successive  PS  converted 
waves  arriving  from  successively  deeper  layers  mil  continually  increase  with  depth 
Figure  IP  shows  this  situation  fora  hypothetical  five-layered  crust.  The  curves  were 
obtained  by  repeatedly  applying  the  Zoeppritz  equations  at  each  interface  to  obtain 


on  Crustal  section  in  the  eastern  part  of  the  Basin  and  Range  province,  according  o 
and  others  (1‘ttiOl.  The  amplitudes  of  successive  PS  converted  waves  do  not  successively 
ne reuse  with  increasing  depth  to  the  interface.  Free  surface  effects  are  not  included. 


the  amplitudes  of  the  various  converted  waves1  The  results  show  that  the  various 
PS  converted  wave  arrivals  should  have  the  largest  amplitude  when  the  angle  of 
incidence  (6)  of  the  parent  P  wave  at  the  base  of  the  crustal  section  is  approxi¬ 
mates  BO  degrees,  as  indicated  in  figure  IP.5  Figure  20  shows  that  for  the  crustal 
section  as  interpreted  by  Berg  and  others  (1900.  figure  10)  in  l  tali,  where  the 
velocity  ratios  do  not  continually  decrease  with  depth,  the  amplitudes  of  the  PS, 
converted  waves  are  theoretically  larger  than  the  amplitudes  of  the  PS,  waves  which 
ir  turn  arc  larger  than  the  amplitudes  of  the  PS,  waves.  Figures  IP  and  20  do  not 


1  In  traveling  to  the  surface  from  the  bottom  layer  (layer  4),  the  wave  designated  PS,  Cotv 
,  inws  as' a  P  wave  through  layers  3  and  2.  and  converts  loan  ,Xi  wave  in  passing  from  lay  er 
,,  laver  1  (surface  laver)-  the  wave  designated  PS,  converts  from  a /  to  an  M  wave  at 
! he  boundary  lietweer  layers  4  and  3  and  continues  through  layers  3,  2,  and  1  as  an  SI 

"  a=V  A  dd  i  t  i  on  a  Fa"  ud  IP  "  a  re  'in  progress  to  ascertain  the  range  of  parameters  for  which  fins 
generalization  may  hold. 
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include*  the*  effects  of  a  free  surface  on  the  /’  and  PS  arrivals,  'Hie  general  conclusions 
would  remain  the  same,  however,  with  the  amplitudes  of  both  the  P  and  PS  ar  ivals 
being  approximately  doubled, 

Figures  21  ami  22.  which  include  the  effects  of  a  free  surface,  show  the  ratios  of  the 


Fin.  21.  Ratio  of  the  horizontal  component  of  the  ground  notion,  ((  *„).  due  to  t he  arrival 
of  PS.  converted  waves  to  the  horizontal  component  of  the  nnlransformed  transmitted  /'  wave, 
(T  ),  for  the  crustal  section  in  the  eastern  part  of  the  Basin  and  Uanse  province,  according 
to  Berg  and  others  (1900),  as  shown  in  figure  20.  Free  surface  effects  are  included. 


Fm.  22.  Ratio  of  the  vertical  conmonent  of  the  ground  motion,  fft'O,  due  to  the  arrival  of 
I‘Sr  converted  waves  to  the  vertical  component  of  the  uiitratisformed  transmitted  H  wave, 
(W  J,  for  the  crustal  section  in  the  eastern  part  of  the  Basin  and  Range  province,  according  to 
Berg  and  others  (19110),  as  shown  in  figure  2tt.  Free  surface  effects  are  included. 
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magnitudes  of  the  horizontal  and  vertical  components  of  the  ground  motion  due  to 
the  arrivals  of  the  various  /\S  converted  waves,  to  the  horizontal  and  vertical  com¬ 
ponents,  respectively,  of  the  untraiisforined  transmitted  P  wave,  for  the  three-layer 
crustal  section  of  Berg  and  others  (1900,  figure  10).  As  Shchter  (19,1:1,  p.  240, 


KiO.  23.  Particle-mot i»n  diagram  of  t lie  i.ik;an  nuclear  explosion  recorded  al  a  distance  of  * 

714  km  in  northern  New  Mexico  with  short -period  BeniofT  seismometers.  Interval  shown  is 
from  (i.44  to  8  59  sec  after  the  firsl  motion,  liigitizing  interval  is  0.05  sec. 


figure  5,  and  p.  248)  points  out,  the  horizontal  ground  motion  of  the  PS  converted 
wave  with  respect  to  the  horizontal  ground  motion  due  to  the  transmitted  un- 
transformed  parent  P  wave  may  be  prominent  at  all  epiceutral  distances,  as  also 
shown  for  PSt  in  figure  21 .  For  the  crustal  section  assumed,  the  resultant  magnitude 
of  these  PS  converted  waves,  however,  is  greatest  for  angles  of  incidence  of  ap- 
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proximately  60  degrees,  corresponding  to  earthquakes  at  an  epic-ntral  distance  of 
approximately  18  degrees.  Bulin  (1060,  p.  519)  noted  that  PS  converted  waves  of 
maximum  amplitude  were  observed  on  seismograms  of  local,  deep-focus  earthquakes 
for  which  th"*  average  value  of  the  angle  of  incidence,  6,  at  the  hase  of  the  crust  was 
56  degrees. 
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Fig.  24.  Particle  motion  diagram  of  the  gnome  nuclear  explosion  recorded  at  Vernal. 
L- 1 all  (vnut).  Interval  shown  is  from  0.4(1  to  7.90  sec  after  the  first  mot  ion.  Digitizing  interval 
is  0.05  sec. 

Vasil Vv  (1956)  discusses  various  PS  converted  waves  recorded  i.s  a  result  of 
explosion  seismic  work  by  the  Geophysical  Institute  of  the  Academy  of  Sciences  of 
the  U.S.S.K.  conducted  in  1954  in  the  Tatar  A.S.S.K.  Vasil’ev  found  that  the  wave 
PmmSi  has  slightly  lower  amplitude  than  the  wave  which  travelled  along  the 
entire  path  as  a  longitudinal  wave.  The  indices  1 , 2,  etc.  refer  to  each  layer.  He  found 
the  amplitude  of  the  converted  wave  PmtSKl  to  be  “at  least  twice  the  value  of  the 
amplitude  of  the  corresponding  longitudinal  wave.”  He  noted  that  these  large 
amplitudes  are  “to  a  certain  extent”  in  contradiction  with  calculated  amplitude 
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values  and  attributed  this  to  a  “partial  screening  and  weakening”  of  the  pun*  longi 
tudinai  waves. 

Selected  copies  of  three-component  seismograms  (made  available  by  Dr.  Carl  F. 
Romney)  from  the  underground  nuclear  explosion  logax,  near  Mercury,  Nevada, 


Fin.  25.  Hend  wave  coefficients  for  compressional  waves  for  density  ratios  (upper  medium/ 
lower  medium)  of  0.50.  0.70  and  0.90  for  velocitv  ratios  (up|)er  medium/lower  medium)  from 
0.05  to  0.95.  as  computed  by  the  method  of  Zvolinskii  (1957,  1958).  Two  coefficient*  given  by 
Heelan  (1953b,  p.  893)  are  shown  for  comparison. 

which  were  recorded  at  temporary  field  stations  out  to  epicentral  distances  of  2000 
km.  were  enlarged,  and  selected  intervals  were  digitized  for  the  study  of  PS  con¬ 
verted  waves.  For  the  analysis  of  particle  motion,  the  records  were  enlarged  to  a 
scale  of  1  mm  =  0.05  sec.  These  enlargements  were  then  digitized  at  1  mm  inter¬ 
vals.  The  resulting  data  were  used  to  plot  particle-motion  diagrams  in  both  (1)  the 
vertical  plane  in  a  radial  direction  from  the  epicenter  and  (2)  the  vertical  plant 
transverse  to  this  direction.  On  a  seismogram  obtained  in  northeastern  Xew  Mexico 
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from  the  loh  an  explosion  at  an  epieentral  distance  of  714  km,  the  first  ten  seconds 
of  arrivals  show  /'  motion  except  for  an  interval  of  approximately  O.G  see  (sec  figure 
2d)  which  begins  7.5  sec  after  the  first  arrival.  This  interval  shows  shear  motion 
(predominantly  .ST).  Various  crustal  models  used  by  ether  workers  (Steinhart 
and  Wool  lard,  1959,  p.  347;  Werth  and  others,  10(i2,  p.  1589)  in  northern  Arizona 
and  New  Mexico  place  the  Mohorovific  discontinuity  at  depths  ranging  from 
34  to  57  km.  I  sing  velocities  found  by  these  investigators  for  the  crust  tieneuth 
northern  Arizona  and  New  Mexico,  the  time  interval  of  7.5  sec  between  the  first 
arrival  and  the  shear  arrival  could  coincide  with  the  arrival  of  a  PS  converted 
wave  generated  at  the  base  of  the  crust  at  a  depth  of  approximately  45  km. 

Selected  events  from  three-component-  recordings  by  the  fleotech  field  recording 
teams  of  the  underground  nuclear  explosion  gnome,  made  available  by  Dr.  C  arl  F. 
Romney,  were  digitized  at  intervals  of  0.05  sec.  Particle-motion  diagrams  were 
prepared  from  the  data.  Station  vnut  (Vernal,  Utah)  shows  a  single  interval  of 
shear  motion,  which  is  provisionally  interpreted  es  the  arrival  of  a  PS  converted 
wave,  in  the  first  10  seconds  (7.50  to  7.85  sec)  from  the  gnome  explosion.  (See  fig¬ 
ure  24.) 

PS  converted  waves  generated  from  head  waves  will  have  amplitude  coefficients 
as  given  by  Heelan  (1953b)  and  Zvolinskii  (1957,  1958).  The  Zvchuskii  method  x 
determining  P- wave  head  coefficients  was  programmed  in  koutuan,  and  head  wave 
coefficients  for  density  ratios  of  0.700,  0.800,  and  0.900  were  calculated  for  velocity 
ratio*  j tsuuprinaikulul  velocity  iu  upper  medium 'Vm  ipwhsrumd  vehrity  in  Luwer 
medium )  ranging  from  0.50  to  0.95.  These  data  are  shown  in  figure  25.  They  suggest 
that  there  might  lx*  sonic  inherent  difficulty  in  correlating  crustal  sections  as 
indicated  from  PS  converted  waves,  with  those*  interpreted  from  conventional 
seismic  refraction  travel-time  curves.  For  example,  figure  25  shows  that  the  lower 
the  velocity  contrast  (higher  the  numerical  value  of  the  velocity  ratio),  the  larger 
w4U  4*  tW  AHLfittuidr  .4  ll<r  l*-*d  warr  Ifcwwvw,  ki*h  vC+wity  Mtiua  art  Irw 
efficient  for  the  conversion  of  P  to  »S’F  waves.  Conceivably,  an  interface  across  which 
the  velocity  contrast  is  low  might  generate  a  largo  amplitude  head  wave,  hut  result 
in  a  PS  converted  wave  of  amplitude  too  small  to  lx*  detected.  It  might,  therefore, 
l>e  difficult  to  expect  the  same  information  about  depths  to,  and  uum’ner  of,  inter¬ 
faces  from  both  methods. 
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